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Foreword 

The purpose of these notes is to present as exhaustively and pedagogically as possible 
some recent mathematical results bearing on the Bose-Einstein condensation phenomenon 
observed in ultra-cold atomic gases. One of the numerous theoretical problems posed by 
these experiments is the understanding of the link between effective models, describing 
the physics with a remarkable precision, and the first principles of quantum mechanics. 
The process leading from the fundamental to the effective theories is often called a mean- 
field limit , and it has motivated a very large number of investigations in theoretical and 
mathematical physics. In this course, we shall focus on one of the methods allowing to 
deal with mean-field limits, which is based on de Finetti theorems. The emergence of the 
mean-field models will be interpreted as a fundamental consequence of the structure of 
physical states under consideration. 

This text will touch on subjects from mathematical analysis, probability theory, con¬ 
densed matter physics, ultra-cold atoms physics, quantum statistical mechanics and quan¬ 
tum information. The emphasis will be on the author’s speciality, namely the analytic 
aspects of the derivation of mean-field equilibrium states. The presentation will thus have 
a pronounced mathematical style, but readers should keep in mind the connection be¬ 
tween the questions adressed here and cold atoms physics, in particular as regards the 
experiments leading to the observation of Bose-Einstein condensates in the mid 90’s. 

A few words about the experiments. Bose-Einstein condensation (BEC) is at the 
heart of a rapidly expanding research field since the mid 90’s. The extreme versatility of 
cold atoms experiments allows for a direct investigation of numerous questions of funda¬ 
mental physics. For more thorough developments in this direction, I refer the reader to 
the literature, in particular to [U 05]071 021 112511152111581 172105] and references therein. 
French readers will find very accessible discussions in [0511571142j . 

The first experimental realizations of BEC took place at the MIT and at Boulder, in 
the groups of W. Ketterle on the one hand and E. Cornell-C. Wieman on the other hand. 
The 2001 Nobel prize in physics was awarded to Cornell-Wieman-Ketterle for this re¬ 
markable achievement. The possibilities opened up by these experiments for investigating 
macroscopic quantum phenomena constitute a cornerstone of contemporary physics. 

A Bose-Einstein condensate is made of a large number of particles (alkali atoms usually) 
occupying the same quantum state. BEC thus requires that said particles be bosons, i.e. 
that they do not satisfy Pauli’s exclusion principle which prevents mutliple occupancy of 
a single quantum state. 

This macroscopic occupancy of a unique low-energy quantum state appears only at very 
low temperatures. There exists a critical temperature T c for the existence of a condensate, 
and macroscopic occupancy occurs only for temperatures T < T c . The existence of such 
a critical temperature was theoretically infered in works of Bose and Einstein [191 63] in 
the 1920’s. Important objections were however formulated: 

(1) The critical temperature T c is extremely low, unrealistically so as it seemed in 
the 1920’s. 

(2) At such a temperature, all known materials should form a solid state, and not be 
gaseous as assumed in Bose and Einstein’s papers. 
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(3) The argument of Bose and Einstein applies to an ideal gas, neglecting interactions 
between particles, which is a serious drawback. 

The first objection could be bypassed only in the 1990’s with the advent of powerful 
techniques such as laser coolingj and evaporative cooling. These allowed to reach tem¬ 
peratures in the micro-Kelvin range in quantum gases trapped by magneto-optics means. 
As for the second objection, the answer has to do with the diluteness of the samples: 
three-particles collisions necessary to initiate the formation of molecules, and ultimately 
of a solid phase, are extremely rare in the experiments. One thus has the possibility of 
observing a metastable gaseous phase during a sufficiently long time for a condensate to 
form. 

The third objection is of a more theoretical nature. Most of the material discussed in 
these notes is part of a research program (of many authors, see references) whose goal is 
to remove that objection. We will thus have the opportunity to discuss it a length in the 
sequel. 

Many agreeing observations have confirmed the experimental realization of BEC: imag¬ 
ing of the atoms’ distribution in momentum/energy space, interference of condensates, 
superfluidity in trapped gases ... The importance thus acquired by the mathematical 
models used in the description of this phenomenon has motivated a vast literature de¬ 
voted to their derivation and analysis. 

Some mathematical questions raised by experiments. In the presence of BEC, 
the gas under consideration can be described by a single wave-functions if : M. d >—> C, 
corresponding to the quantum state in which all particles reside. A system of N quantum 
particles should normally be described by a N -particle wave-function : K . dN *—> C. 
One thus has to understand how and why can this huge simplification be justified, that is 
how the collective behavior of the N particles emerges. The investigation of the precision 
of this approximation, whose practical and theoretical consequences are fundamental, is a 
task of extreme importance for theoretical and mathematical physicists. 

One may ask the following questions: 

(1) Can one describe the ground state (that is the equlibrium state at zero tempera¬ 
ture) of an interacting Bose system with a single wave-function if ? 

(2) Start from a single wave-function and let the system evolve along the natural 
dynamics (A r -body Schrodinger flow). Is the single wave-function description pre¬ 
served by the dynamics ? 

(3) Can one rigorously prove the existence of a critical temperature T c under which the 
finite temperature equilibrium states may be described by a single wave-function ? 

These questions are three aspects of the third objection mentioned in the preceding 
paragraph. We thus recall that the point is to understand the BEC phenomenon in the 
presence of interactions. The case of an ideal gas is essentially trivial, at least for questions 
1 and 2. 

One should keep in mind that, in the spirit of statistical mechanics, we aim at jus¬ 
tifying the single wave-function description asymptotically in the limit of large particle 
numbers, under appropriate assumptions on the model under consideration. Ideally, the 
assumptions should reduce to those ensuring that both the iV-body model one starts from 

■*■1997 Nobel prize in physics: S. Chu-W. Phillips-C. Cohen-Tannoudji. 
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and the 1-body model one arrives at are mathematically well-defined. Let us note that, 
for interacting quantum particles, the former is always linear, while the latter is always 
non-linear. 

Many recent results presented in the sequel are generalizations to the quantum case 
of better known results of classical statistical mechanics. Related questions indeed occur 
also in this simpler context. For pedagogical reasons, some notions on mean-field limits 
in classical mechanics will thus be recalled in these notes. 

This course deals with question number 1, and we will be naturally lead to develop tools 
of intrinsic mathematical interest. One may use essentially two approaches: 

• The first one exploits particular properties of certain important physical models. 
It thus applies differently to different models, and under often rather restrictive 
assumptions, in particular as regards the shape of the inter-particle interactions. 
A non-exhaustive list of references using such ideas is [ 173 . 172 . 167 . 168 . 179] for 
the classical case, and [15] 1133] I176[ [86 , 178. 125 . 123 ] for the quantum case 

• The second one is the object of these lectures. It exploits properties of the set 
of admissible states, that is of V-body wave-functions Tat. That we consider 
bosonic particles implies a fundamental symmetry property for these functions. 
This approach has the merit of being much more general than the first one, and 
in many cases to get pretty close to the “ideally minimal assumptions” for the 
validity of the mean-field approximation mentioned before. 

A possible interpretation is to see the mean-field limit as a parameter regime 
where correlations between particles become negligible. We will use very strongly 
the key notion of bosonic symmetry. There will be many opportunities to discuss 
the literature in details, but let us mention immediately [1461 I3T1198: 99. Tdl . T70J 
and [TO] (71, 1154 : 11591 11131 115] for applications of these ideas in classical and 
quantum mechanics respectively. 

The distinction between the two philosophies is of course somewhat artificial since one will 
often benefit from borrowing ideas to both, see e.g. [ 149] for a recent example. 

To keep these notes reasonably short, question 2 will not be treated at all, although 
a vast literature exists, see e.g. [92] [79] [801 HM m M ESI 0 [691 ESI d65l E2 dm [[SI] 
and references therein, as well as the lecture notes [821111]. We note that the quantum de 
Finetti theorems that we will discuss in the sequel have recently proved useful in dealing 
with question 2, see El El El [361 ESI - The use of classical de Finetti theorems in a dynamic 
framework is older [1891119(11 fi91j 

As for question 3, this is a famous open problem in mathematical physics. Very little is 
known at a satisfying level of mathematical rigor, but see however [nzunj. We will touch 
on questions raised by taking temperature into account only very briefly in Appendix iBl 
in a greatly simplified framework. This subject is further studied in the paper [ 114] , 

Plan of the course. 

These notes are organized as follows: 

• A long introduction, Chapter [IJ recalls the basic formalism we shall need to give a 
precise formulation to the problems we are interested in. We will start with classical 
mechanics and then move on to quantum aspects. The question of justifying 
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the mean-field approximation for equilibrium states of a given Hamiltonian will 
be formulated in both contexts. The proof stategy that will occupy the bulk 
of the notes will be described in a purely formal manner, so as to introduce as 
quickly as possible the de Finetti theorems that will be our main tools. Similarities 
between the classical and quantum frameworks are very strong. Differences show 
up essentially when discussing the proofs of the fundamental de Finetti theorems. 

• Chapter [2] is essentially independent from the rest of the notes. It contains the 
analysis of classical systems: proof of the classical de Finetti theorem (also called 
Hewitt-Savage theorem), application to equilibrium states of a classical Hamilton¬ 
ian. The proof of the Hewitt-Savage theorem we will present, due to Diaconis and 
Freedman, is purely classical and does not generalize to the quantum case. 

• We start adressing quantum aspects in Chapter [3l Two versions (strong and 
weak) of the quantum de Finetti theorem are given without proofs, along with 
their direct applications to “relatively simple” bosonic systems in the mean-field 
regime. Section 13.41 contains a discussion of the various versions of the quantum 
de Finetti theorem and outlines the proof strategy that we shall follow. 

• Chapters [Hand [5] contain the two main steps of the proof of the quantum de Finetti 
theorem we choosed to present: respectively “explicit construction and estimates 
in finite dimension” and “generalization to infinite dimensions via localization in 
Fock space”. The proof should not be seen as a black box: not only the final result 
but also the intermediary constructions will be of use in the sequel. 

• Equiped with the results of the two previous chapters, we will be able to give in 
Chapter [6] the justification of the mean-field approximation for the ground state of 
an essentially generic bosonic system. Contrarily to the case treated in Chapter [3l 
the quantum de Finetti theorem will not be sufficient in this case, and we will have 
to make use of some of the ingredients introduced in Chapter El 

• The mean-field limit is not the only physically relevant one. In Chapter [7] we will 
study a dilute regime in which the range of the interactions goes to 0 when N —>• oo. 
In this case one obtains in the limit Gross-Pitaevskii (or non-linear Schrodinger) 
functionals with local non-linearities. We will present a strategy for the derivation 
of such objects based on the tools of Chapter [H 

The main body of the text is supplemented with two appendices containing each an 
unpublished note of Mathieu Lewin and the author. 

• Appendix [A] shows how, in some particular cases, one may use the classical de 
Finetti theorem to deal with quantum problems. This strategy is less natural (and 
less efficient) than that presented in Chapters [3] and El but it has a conceptual 
interest. 

• Appendix IBldeviates from the main line of the course since the Hilbert spaces under 
consideration will be finite dimensional. In this context, combining a large tem¬ 
perature limit with a mean-field limit, one may obtain a theorem of senri-classical 
nature which gives examples of de Finetti measures not encountered previously. 
This will be the occasion to mention Berezin-Lieb inequalities and their link with 
the considerations of Chapter EJ 

Acknowledgements. The motivation to write the french version of these notes came 
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1. Introduction: Problems and Formalism 

Here we describe the mathematical objects that will be studied throughout the course. 
Our main object of interest will be many-body quantum mechanics, but the analogy with 
some questions of classical mechanics is instructive enough for us to also describe that 
formalism. Questions of units and dimensionality are systematically ignored to simplify 
notation. 

1.1. Statistical mechanics and mean-field approximation. For pedagogical reasons 
we shall recall some notions on mean-field limits in classical mechanics before going to the 
quantum aspects, relevant to the BEC phenomenon. This paragraph aims at fixing nota¬ 
tion and reviewing some basic concepts of statistical mechanics. We will limit ourselves 
to the description of the equilibrium states of a classical system. Dynamic aspects are 
voluntarily ignored, and the reader is refered to [82j for a review on these subjects. 

Phase space. The state of a classical particle is entirely determined by its position x and 
its speed v (or equivalently its momentum p ). For a particle living in a domain D C R d 
we thus work in the phase space D x R d , the set of possible positions and momenta. For 
a iV-particle system we work in £l N x W IN . 

Pure states. We call a pure state one where the positions and momenta of all particles are 
known exactly. Equilibrium states at zero temperature for example are pure: in classical 
mechanics, uncertainty on the state of a system is only due to “thermal noise”. 

For a IV-particle system, a pure state corresponds to a point 

(X;P) = (xi,.. .,x N ;pi,... ,p N ) € ft N x R dN 

in phase-space, where the pair (xf,pi) gives the position and momentum of particle num¬ 
ber i. Having in mind the introduction of mixed states in the sequel, we will identify a 
pure state with a superposition of Dirac masses 

Px ; p = ^ S Xa] p a - (1-1) 

crgEjv 

This identification takes into account the fact that real particles are indistinguishable. 
One can actually not attribute the pair ( xp,pi ) of position/momentum to any one of the 
N particles in particular, whence the sum over the permutation group in (HHD- Here 
and in the sequel our notation is 

X(j (^cr(l) 5 • • • ? ^cr(N) ) 

(^cr(l) •> • • • ?Pcr(AI))* (1‘^) 

Saying that the system is in the state px-.p means that one of the particles has position 
and momentum ( xp,pi),i = 1...N, but one cannot specify which because of indistin- 
guishability. 

Mixed states. At non-zero temperature, that is when some thermal noise is present, one 
cannot determine with certainty the state of the system. One in fact looks for a statistical 
superposition of pure states, which corresponds to specifying the probability that the 
system is in a certain pure state. One then speaks of mixed states, which are the convex 
cominations of pure states, seen as Dirac masses as in Equation (11.11) . The set of convex 
combinations of pure states of course corresponds to the set of symmetric probability 
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measures on phase-space. A general IV-particles mixed state is thus a probability measure 
p, N € V s (£l N x R dN ) satisfying 

dfx N (X-,P) = dn N (X a \P a ) (1.3) 

for all permutations a £ S^r. One interprets fi N (X;P) as the probability density that 
particle i has position X{ and momentum pi, i = 1... N. Pure states of the form D are 
a particular type of mixed states where the statistical uncertainty is reduced to zero, up 
to indistinguishability. 


Free energy. The energy of a classical system is specified by the choice of a Hamiltonian, 
a function over phase-space. In non-relativistic mechanics, the kinetic energy of a particle 
of momentum p is always m\p\ 2 /2. Taking m = 1 to simplify notation, we will consider 
an energy of the form 

H N (X-,P):=^2 l -^- + ^2 v (x j ) + A £ w(xi- Xj ) (1.4) 

3 =1 3 = 1 l<i<j<N 

where V is an external potential (e.g. electrostatic) felt by all particles and w a pair- 
interaction potential that we will assume symmetric, 

w(—x) = w(x). 


The real parameter A sets the strength of interparticle interactions. We could of course 
add three-particles, four-particles etc ... interactions, but this is seldom required by the 
physics, and when it is, there is no additional conceptual difficulty. 

The energy of a mixed state p N € V S (Q N x R dN ) is then given by 

£\p N }-= f H N {X-P)dp N {X-P) (1.5) 

Jfl N x R dN 

and by symmetry of the Hamiltonian this reduces to H^{X\P) for a pure state of the 
form ED- At zero temperature, equilibrium states are found by minimizing the energy 
functional (11,51) : 

E(N) = inf {s\ja n \,h n £ V a (Sl N x M d7V )} (1.6) 

and the infimum (ground state energy) is of course equal to the minimum of the Hamil¬ 
tonian H^. It is attained by a pure state ED where {X ; P) is a minimum point for Hn 
(in particular P = (0,..., 0)). 

In presence of thermal noise, one must take the entropy 

S\f* N ]-=~ f dp N (X;P) log(p N (X; P)) (1.7) 

Jn N xR dN 

into account. This is a measure of the degree of uncertainty on the state of the system. 
Note for example that pure states have the lowest possible entropy: = —oo if p N 

if of the form ED. At temperature T, one finds the equilibrium state by minimizing the 
free-energy functional 




[ H N (X-P)dp N (X-P)+T [ dn N {X;P)log(p N {X;P)) (1.8) 
Jn N x«. dN Jn N xM. dN 
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which amounts to saying that the more probable states must find a balance between having 
a low energy and having a large entropy. We shall denote 

F{N) = inf {x[n N ],V N G V s (tt N x R dJV )} (1.9) 

without specifying the temperature dependence. A minimizer must be a sufficiently regular 
probability so that (minus) the entropy is finite 


Momentum minimization. In the absence of a prescribed relation between the position and 
the momentum distribution of a classical state, the minimization in momentum variables 
of the above functionals is in fact trivial. A state minimizing (11.51) is always of the form 


Hn — dp =o ® ^ x=x ° 

ceEjv 

where X° is a minimum point for Hn(X; 0,... ,0), i.e. particles are all at rest. We are 
thus reduced to looking for the minimum points of Hn(X; 0,... , 0) as a function of X. 

The minimization of (ESI) leads to a Gaussian in momentum variables mutliplied by a 
Gibbs state in position variable^! 

/Tv = 77 ^ exp ® z“ exp (-^iv(AT;0,... ,0) 

Momentum variables thus no longer intervene in the minimization of the functionals de¬ 
termining equilibrium states and they will be completely ignored in the sequel. We will 
keep the preceding notation for the minimization in position variables: 



N 

H n (X) = y ^V(xj) + A w(xi-Xj) 

j= 1 1<2<7 <7V 

£[v N ] = I H N (X)dfi N (X ) 

Jn N 

?[?*]=[ H N {X)dn N {X) + T [ d» N (X) log(n N (X)) (1.10) 

Jn N Jn N 

where G V S (Q N ) is a symmetric probability measure in position variables only. 


Marginals, reduced densities. Given a IV-particles mixed state, it is very useful to consider 
its marginals, or reduced densities, obtained by integrating out some variables: 

IJ-n\x i,...,x n )= / fj,(x 1 ,...,x n ,x' n+1 ,...,x' N )dx' n+1 ...dx' N eT s (n n ). ( 1 . 11 ) 

(n) 

The n—th reduced density fi y N is interpreted as the probability density for having one 
particle at x\, one particle at X2, etc... one particle at x n . In view of the symmetry of 
the choice of which N — n variables over which one integrates in Definition (11.111) is 
irrelevant. 


2 The partition functions Zp and Zn normalize the state in L 1 . 
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A first use of these marginals is a rewriting of the energy using only the first two 
marginal^: 

£[H N ] = n J V{x)d^\x) + ^~ JJ w ( x ~ 

= 11 (^V(x) + ^V{y) + \ N ( N 9 ^ w{x~y)) d^\x,y), ( 1 . 12 ) 

where we used the symmetry of the Hamiltonian. 

Mean-field approximation. Solving the above minimization minimization problems is a 
very difficult task when the particle number N gets large. In order to obtain more tractable 
theories from which useful information can be extracted, one often has to rely on approx¬ 
imations. The simplest and most well-known of these is the mean-field approximation. 
One can introduce it in several ways, the goal being to obtain a self-consistent one-body 
problem starting from the IV-body problem. 

Here we follow the “molecular chaos” point of view on mean-field theory, originating in 
the works of Maxwell and Boltzmann. The approximation consists in assuming that all 
particles are independent and identically distributed (iid). We thus take an ansatz of the 
form 

N 

/Tv( x 'i, .. .,x N ) = p® N (x 1 ,.. .,x N ) = (1.13) 

where p S is a one-body probability density describing the typical behavior of one 

of the iid particles under consideration. 

The mean-held energy and free-energy functionals are obtained by inserting this ansatz 
in m or (USD. The mean-held energy functional is thus 

S MF [p] = N~ 1 £[p 0N ] = f V(x)dp(x) + X N [[ w(x - y)dp{x)dp{y). (1.14) 

Jq a JJuxq 

We shall denote E MF its inhmum amongst probability measures. In a similar manner, the 
mean-held free energy functional is given by 

= N -ljr[ p ®N} 

= f V(x)dp(x) + X N ff w(x - y)dp(x)dp(y) + T [ plogp (1.15) 

and its inhmum shall be denoted F mf . The term “mean-held” is motivated by the fact 
that (I1.14P corresponds to having an interaction between the particles’ density p and the 
self-consistant potential 

p*w = / w(. — y)dp(y) 

Ja 

whose gradient is the so-called mean-held. 


^More generally, an energy depending only on a n-body potential may be rewritten by using the n-th 
marginal only. 
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1.2. Quantum mechanics and Bose-Einstein condensation. After this classical in- 
temezzo we can now introduce the quantum objects that are the main theme of this 
course. We shall be content with a rapid review of the basic principles of quantum me¬ 
chanics. Other “mathematician-friendly” presentations may be found in [12^1H' 88 ] . Our 
discussion of the relevant concepts is in places voluntarily simplified. 

Wave-functions and quantum kinetic energy. One of the basic postulates of quantum me¬ 
chanics is the identification of pure states of a system with normalized vectors of a complex 
Hilbert space Sj. For particles living in the configuration space the relevant Hilbert 
space is T 2 (M rf ), the space of complex square-integrable functions over M. d . 

Given a particle in the state if € L 2 (R rf ), one identifies \if \ 2 with a probability density: 
\if(x)\ 2 gives the probability for the particle to be located at x. We thus impose the 
normalization 

[ M 2 = i- 

J R d 

We thus see that, even in the case of a pure state, one cannot specify with certainty the 
position of the particle. More precisely, one cannot simultaneously specify its position and 
its momentum. This uncertainty principle is the direct consequence of another fundamen¬ 
tal postulate: \if\ 2 gives the momentum-space probability density of the particle, where if 
is the Fourier transform of if. 

In quantum mechanics, the (non-relativistic) kinetic energy of a particle is thus given 
by 

[ ^~\^(p)\ 2 dp = [ \\S/if{x)\ 2 dx. (1.16) 

That the position and momentum of a particle cannot be simultaneously specified is a 
consequence of the fact that it is impossible for both \if\ 2 and \if\ 2 to converge to a 
Dirac mass. A popular way of quantifying this is Heisenberg’s uncertainty principle: for 
all xq € 

( f |V^(x)| 2 dx) ( f \x — xo\ 2 \if(x)\ 2 dx ) > C. 

\J R<* / \JR d J 

Indeed, the more precisely the particle’s position is known, the smaller the second term of 
the left-hand side (for a certain xo). The first term of the left-hand side must then be very 
large, which, in view of (11.161) rules out the possibility for the momentum distribution to 
be concentrated around a single p$ € M. d . 

For many applications however (see inn for a discussion of this point), this inequality is 
not sufficient. A better way of quantifying the uncertainty principle is given by Sobolev’s 
inequality (here in its 3D version): 

[ |vvf>c(7 IVfY /3 . 

J R 3 \J R 3 / 

If the position of the particle is known with precision, \if\ 2 must approach a Dirac mass, 
and the right-hand side of the above inequality blows up. So do the integrals (11.161) . with 
the same interpretation as previously. 
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Bosons and Fermions. For a system of N quantum particles living in M. d , the appropriate 
Hilbert space is L 2 (R dN ) ~ 0 jV L 2 (M d ). A pure state is thus a certain T E L 2 (M. dN ) 
where ^(xi,..., xn)\ 2 is interpreted as the probability density for having particle 1 at x\, 
particle N &t xjy. As in classical mechanics, the indistinguishability of the particles 
imposes 

l'F(A)! 2 = |'J / (A cr )| 2 (1.17) 

for any permutation d G S]\r. This condition is necessary for indistinguishable patricles, 
but it is in fact not sufficient for quantum particles. To introduce the correct notion, we 
denote U a the unitary operator interchanging particles according to a E Sjv- 

U a u\ <g>... ® un = iV(i) <8> ■ ■ ■ <8> w CT ( at) 

for all ui,...,un E L 2 ( R d ), extended by linearity to L 2 (R dN ) ~ 0 iV L 2 (R d ) (one may 
construct a basis with vectors of the form u\ <8> ■ ■ • <83 un). For T E L 2 (W iN ) to describe 
indistinguishable particles we have to require that 

(T, AT) L 2( Rt i2V) = (Uer^, AUcr^) L 2^dN^ (1-18) 

for any bounded operator A acting on L 2 (W lN ). Details would lead us to far, but suffice it 
to say that condition (11,181) corresponds to asking that any measure (corresponding to an 
observable A) on the system must be independent of the particles’ labeling. In classical 
mechanics, all possible measurements correspond to bounded functions on phase-space 
and thus (11.31) alone guarantees the invariance of observations under particle exchanges. 
In quantum mechanics, observables include all bounded operators on the ambient Hilbert 
space, and one must thus impose the stronger condition (11.181) . 

An important consequence^ of the symmetry condition (11.181) is that a system of in¬ 
distinguishable quantum particles must satisfy one of the following conditions, stronger 
than (11.171) : either 

T(A) = T(X ct ) (1.19) 

for all X € M. dN and a E S^r, or 

T(X) = £(a)^(X a ) (1.20) 

for all X E and a E Sat, where e(<r) is the sign of the permutation a. One refers 
to particles described by a wave-function satisfying (11.191) (respectively (11.201) 1 as bosons 
(respectively fermions). These two types of fundamental particles have a very different 
behavior, one speaks of bosonic and fermionic statistics. For example, fermions obey the 
Pauli exclusion principle which stipulates that two fermions cannot simultaneously occupy 
the same quantum state. One can already see that (11.201) imposes 

T(xi, ... ,Xi,...,Xj,.. .,x N ) = -T(xi, ...,Xj,...,Xi,.. .,x N ) 

and thus (formally) 

T(xi, ... ,Xi,... ,Xi,.. .,x N ) = 0 

which implies that it is impossible for two fermions to occupy the same position x,;. One 
may consult [ 124) for a discussion of Lieb-Thirring inequalities, which are one of the most 
important consequences of Pauli’s principle. 

4 This is an easy but non-trivial exercise. 
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Concretely, when studying a quantum system made of only one type of particles, one 
has to restrict admissible pure states to those being either of bosonic of fermionic type. 
One thus works 

• for bosons, in L 2 (R dN ) ~ T 2 (R d ), the space of symmetric square-integrable 

wave-functions, identified with the symmetric tensor product of N copies of L 2 (R d ). 

• for fermions, in L 2 s (R dN ) ~ (^) as L 2 (W l ), the space of anti-symmetric square- 
integrable wave-functions, identified with the anti-symmetric tensor product of N 
copies of L 2 (R d ). 

As the name indicates, Bose-Einstein condensation can occur only in a bosonic system, 
and this course shall thus focus on the first case. 

Density matrices, mixed states. We will always identify a pure state T € L 2 (R dN ) with 
the corresponding density matrix, i.e. the orthogonal projector onto \k, denoted |\k) (\k|. 
Just as in classical mechanics, mixed states of a system are by definition the statistical 
superpositions of pure states, that is the convex combinations of orthogonal projectors. 
Using the spectral theorem, it is clear that the set of mixed states conciides with that of 
positive self-adjoint operators having trace 1: 

s(L 2 (R dN )) = |r e &\L 2 (R dN )),T = r*,r > o,Tfr = 1 } ( 1 . 21 ) 

where © 1 (iJ) is the Schatten-von Neumann class [161 . 174 . 1831 of trace-class operators on a 
Hilbert space S). To obtain the bosonic (fermionic) mixed states one considers respectively 

s(L 2 s/as (R dN )) = {re 6\L 2 s/as (R dN )) 1 r = r*,r> o,ivr = 1 }. 

Note that in the vocabulary of density matrices, bosonic symmetry consists in imposing 

U a T = T (1.22) 

whereas fermionic symmetry corresponds to 

U a T = e{a) T. 

One sometimes considers the weaker symmetry notion 

U a TU* = T (1.23) 

which is for example satisfied by the (unphysical) superposition of a bosonic and a 
fermionic state. 

Energy functionals. The quantum energy functional corresponding to the classical non- 
relativistic Hamiltonian (11.41) is obtained by the substitution 

p o —zV, (1-24) 


consistent with the identification (11.161) for the kinetic energy of a quantum particle. The 
quantized Hamiltonian is a (unbounded) operator on L 2 (R dN ): 
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where —A j = (—iVj) 2 corresponds to the Laplacian in the variable Xj £ R d . The corre¬ 
sponding energy for a pure state 'k € L 2 (R dN ) is 

£[^] = (^,H N ^) L 2^ R dNy (1.26) 

By linearity this generalizes to 

£[T] = Tr L2{wdN) [H N T}. (1.27) 

in the case of a mixed state T € S(L 2 (R dN )). At zero temperature, the equilibrium state 
of the system is found by minimizing the above energy functional. In view of the linearity 
of (11.271) as a function of T and by the spectral theorem, it is clear that one may restrict 
the minimization to pure states: 

E s/as (N) = inf {£[T],T £ S(L 2 s/as (R dN ))} 

= inf {£[*],* € L 2 s/as (R dN ), II^IIl^R'W) = l} • 

(1.28) 

Here we use again the index s (respectively as) to denoted the bosonic (respectively 
fermionic) energy. In the absence of an index, we mean that the minimization is performed 
without symmetry constraint. However, because of the symmetry of the Hamiltonian, one 
can in this case restrict the minimization to mixed states satisfying (11.231) . or to wave- 
functions satisfying (11.171) . 

In the presence of thermal noise, one must add to the energy a term including the 
Von-Neumann entropy 

S[T] = -Tr L 2 (M dAr)[riogr] = -^Ojlogaj (1.29) 

j 

where the aj’s are the eigenvalues (real and positive) of T, whose existence is guaranteed 
by the spectral theorem. Similarly to the classical entropy dm Von-Neumann’s entropy 
is minimized (<S[r] = 0) by pure states, i.e. orthongonal projectors, which of course only 
have one non-zero eigenvalue, equal to 1. 

The free-energy functional at temperature T is then 

jqr] = £[T] - TS[T] (1.30) 

and minimizers are in general mixed states. 

Alternative forms for the kinetic energy: magnetic fields and relativistic effects. In the 
classical case, we have only introduced the simplest possible form of kinetic energy. The 
reason is that for the minimization problems that shall concern us, the kinetic energy plays 
no role. Other choices for the relation between p and the kinetic energ}0 are nevertheless 
possible. In quantum mechanics the choice of this relation is crucial even at equilibrium be¬ 
cause of the non-trivial minimization in momentum variables. Recalling (11.161) and (|1.24l) . 
we also see that in the quantum context, different choices will lead to different functional 
spaces in which to set the problem. 

Apart from the non-relativistic quantum energy already introduced, at least two gener¬ 
alizations are physically interesting: 

^Other dispersion relations. 
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• When a magnetic field B : i —> M interacts with the particles, one replaces 

p —iV + A (1.31) 

where A : M. d i— > M. d is the vector potential , satisfying 

B = curl A. 

Of course B determines A only up to a gradient. The choice of a particular A is 
called a gauge choice. The kinetic energy operator, taking the Lorentz force into 
account, becomes in this case 

(p + A) 2 = (-zV + Af = - (V - iA) 2 , 

called a magnetic Laplacian. 

This formalism is also appropriate for particles in a rotating frame: calling X 3 the 
rotation axis one must then take A = £l(—x 2 , x\ , 0) with fl the rotation frequency. 
This corresponds to taking the Coriolis force into account. In this case one must 
also replace the potential V(x) by V(x) — fl 2 |x | 2 to account for the centrifugal 
force. 

• When one wishes to take relativistic effects into account, the dispersion relation 

becomes _ 

Kinetic energy = cv'p 2 + rn 2 c 2 — me 2 

with m the mass and c the speed of light in vacuum. Choosing units so that c = 1 
and recalling (11.2411 . one is lead to consider the kinetic energy operator 

\Jp 1 + m 2 — m = \/—A + m 2 — m (1-32) 

which is easily defined in Fourier variables for example. In the non-relativitic limit 
where \p\ <C m = me we formally recover the operator —A to leading order. A 
caricature of (11.321) is sometimes used, corresponding to the “extreme relativistic” 
case |p| me where one takes 

V / P^ = V Z A (1-33) 

as kinetic energy operator. 

• One can of course combine the two generalizations to consider relativistic particles 
in a magnetic field, using the operators 

\J (—iV + A) 2 + m 2 — m and | — zV + A\ 

based on the relativistic dispersion relation and the correspondance (11.311) . 

Reduced density matrices. It will be useful to define in the quantum context a concept 
similar to the reduced densities of a classical state. Given T £ S(L 2 (W iN )), one defines 
its n-th reduced density matrix by taking a partial trace on the last N — n particles: 

r< n) =T 'tu+i^nT, (1.34) 

which precisely means that for any bounded operator A n acting on L 2 (M rfn ), 

TV L2(Mdn) [A n rW] := Tr L2{RdN) [A n ® 1®^] 

where 1 is the identity on L 2 (W l ). The above definition is easily generalized to any 
Hilbert space, but in the case of L 2 , the reader is maybe more familiar with the following 
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equivalent definition: We identify T with its kernel, i.e. the function T(X;Y) such that 
for all T G L 2 (R dN ) 

TT= / r (X-,Y)V(Y)dY. 

J R dN 

One can then also identify with its kernel 


— / r(xi ,... , x n , z n -j -1, • • •, Z]\i ,y i,..., y n , z n -\- 1, • • •, zjv)d^n+i • • • dzjy. 

In the case where T has a symmetry, bosonic or fermionic, the choice of the variables over 
which we take the partial trace is arbitrary. Let us note that even if one starts from a 
pure state, the reduced density matrices are in general mixed states. 

Just as in (11.121) one may rewrite the energy (11.271) in the form 


£[F] — YTr L 2( R d) 
= Tr L 2 ( R 2 d) 


(-±a + e) 


+ —— Tr i 2 (R 2d ) [m(x - y)r (2) ] 


y + V(x)~ ±A y + V(y)) + \^^j-^-w(x - y)^j 


pi 2 ) 

(1.35) 


Mean-field approximation. Even more than so in classical mechanics, actually solving the 
minimization problem (11.281) for large N is way too costly, and it is necessary to rely on 
approximations. This course aims at studying the simplest of those, which consists in 
imitating (11.131) by taking an ansatz for iid particles, 

T(xi, ... ,x N ) = u® N (x i,... ,x N ) = u(x i) . ..u(x N ) (1.36) 

for a certain u € L 2 (M d ). Inserting this in the energy functional (11.351) we obtain the 
Hartree functional 

£ n [u] = N~ 1 £[u® n ] 

= f (h^u ] 2 + V\u\A + ff \u{x)\ 2 w(x-y)\u{y )\ 2 (1.37) 

JR d J 1 J J R d xR d 

with the corresponding minimization problem 

e H = inf {f H [«], hllL 2 (R d) = l} • (1.38) 

Note that we have transformed a linear problem for the fV-body wave-function (since the 
energy functional is quadratic, the variational equation is linear) into a cubic problem for 
the one-body wave-function u (the energy functional is quartic and hence the variational 
equation will be cubic). 

Ansatz (|1.36l) is a symmetric wave-function, appropriate for bosons. It corresponds to 
looking for the ground state in the form of a Bose-Einstein condensate where all particles 
are in the state u. Because of Pauli’s principle, fermions can in fact never be completely 
uncorrelated in the sense of (11.361) . The mean-held ansatz for fermions rather consists in 
taking 

^(aq ,...,x N ) = det (uj(x k )) 
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with orthonormal functions u±,, un (orbitals of the system). This ansatz (Slater deter¬ 
minant) leads to the Hartree-Fock functional that we shall not discuss in these notes (a 
presentation in the same spirit can be found in [166]) 

Note that for classical particles, the mean-field ansatz (11.131) allows taking a mixed 
state p. To describe a bosonic system, one always takes a pure state u in the ansatz ([1.36D . 
which calls for the following remarks: 

• If one takes a general 7 € S(L 2 (R d )), the IV-body state defined as 

T = 7 ® N (1.39) 

has the bosonic svmmetrv (11.221) if and onlv if 'v is a Dure state (see e.e. EU), 
7 = |u)(it|, in which case T = |ri® Ar )(rt® Ar | and we are back to Ansatz ([1.361) . 

• In the case of the energy functional (11.351) . the minimization problems with and 
without bosonic symmetry imposed coincide (see e.g. EMI Chapter 3]). One can 
then minimize with no constraint and find the bosonic energy. This remains true 
when the kinetic energy is (jl.321) or (11.331) but is notoriously wrong in the presence 
of a magnetic field. 

• For some energy functionals, for example in presence of a magnetic or rotation 
field, the minimum without bosonic symmetry is strictly lower than the minimum 
with symmetry, cf |lT5j . The ansatz ([1.391) is then appropriate to approximate 
the problem without bosonic symmetry (in view of the Hamiltonian’s symmetry, 
one may always assume the weaker symmetry ([1.231) 1 and one then obtains a 
generalized Hartree functional for mixed one-body states 7 € S(L 2 (M . d )): 

£n\l] = Tr L 2 (Rd) [(-^A + V) 7 ] + X N 2 Tr L 2 (R2d) [w(x - y)7 ® 2 ] . (1.40) 


The next two sections introduce, respectively in the classical and quantum setting, the 
question that will occupy us in the rest of the course: 

Can one justify, in a certain limit, the validity of the mean-field ansatze ([1.131) 
and ([1.361) to describe equilibrium states of a system of indistinguishable particles ? 


1.3. Mean-field approximation and the classical de Finetti theorem. 

Is it legitimate to use the simplification (11.131) to determine the equilibrium states of a 
classical system ? Experiments show this is the case when the number of particles is large, 
which we mathematically implement by considering the limit N —>• 00 of the problem at 
hand. 

A simple framework in which one can justify the validity of the mean-field approximation 
is the so-called mean-fied limit , where one assumes that all terms in the energy (11,51) have 
comparable weights. In view of (11.121) . this is fulfilled if A scales as IV” 1 , for example 


A 


1 

IV- 1’ 


(1.41) 


in which case one may expect the ground state energy per particle N~ 1 E(N ) to have a 
well-defined limit. The particular choice (ll.41|) helps in simplifying certain expressions, 
but the following considerations apply as soon as A is of order N^ 1 . We should insist on the 
simplification we introduced in looking at the N —>• 00 limit under the assumption (ll.4ip . 
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This regime is far from covering all the physically relevant cases. It however already leads 
to an already very non-trivial and instructive problem. 

The goal of this section is to dicuss the mean-field approximation for the zero- 
temperature equilibrium states of a classical system. We are thus looking at 


E(N) = inf 



H]\rdn N , n N € V s (£l N ) 


(1.42) 


with 

N 1 

Hn(X) = ^V ( Xj ) + N _ i ^2 w(xi-Xj). (1.43) 

j =1 < AT 

We will sketch a formal proof of the validity of the mean-field approximation at the level 
of the ground state energy, i.e. we shall argue that 


= £ MF = inf {£ MF Ip\,P € V(Sl)} 


(1.44) 


where the mean-held functional is obtained as in (11.141) . taking (11.411) into account: 

£ MF [p\ = ( Vdp+ \ [[ w(x -y)dp(x)dp{y). (1.45) 

J£1 ^ J JPixPl 


Estimate (11.441) for the ground state energy is a first step, and the proof scheme in fact 
yields information on the equilibrium states themselves. To simplify the presentation, we 
shall postpone the discussion of this aspect, as well as the study of the positive temperature 
case, to Chapter [21 


Formally passing to the limit. Here we wish to make the “algebraic” structure of the 
problem apparent. The justification of the manipulations we will perform requires analysis 
assumptions that we shall discuss in the sequel of the course, but that we ignore for the 
moment. 

We start by using (I1.12|) and (|1.41l) to write 

T = \ inf { f[ H 2(x,y)diu,^\x,y), n N eV s (n N )\ 

where H 2 is the two-body Hamiltonian defined as in (11 .43|) and fx y N is the second marginal 
of the symmetric probability pi N . Since the energy depends only on the second marginal, 
one may see the problem we are interested in as a constrained optimization problem for 
two-body symmetric probability measures: 

a ^ 2 ( x,y )^ (2} ( x,y ) ,M(2) ep ^ ) } 

with 

vP = {^ (2) € V s (n 2 ) \3pi N € Vs(n N ),^ = ^ } } 

the set of two-body probability measures one can obtain as the marginals of a IV-body 
state. Assuming that the limit exists (first formal argument) we thus obtain 

Jim \ Jim inf { [[ H 2 (x, y)dp 2) (x, y), p 2) € Vp \ 

n^-oo iV 2 Af->oo IJJnxn J 


(1.46) 




















22 


NICOLAS ROUGERIE 


and it is very tempting to exchange limit and infimum in this expression (second formal 
argument) in order to obtain 


lim 

TV—>■ oo 


E(N) 

N 



H 2 (x, y)dp ly2 \x 1 y), G lim V 




TV—^oo 


( 2 ) 

N 


We have here observed that the energy functional appearing in (jl.4611 is actually indepen¬ 
dent of N. All the ^-dependence lies in the constrained variational space we consider. 
This suggests that a natural limit problem consists in minimizing the same functional but 

on the limit of the variational space, as written above. 

( 2 ) ( 2 ) 

To give a meaning to the limit of V y N , we observe that the sets V N form a decreasing 
sequence 


,( 2 ) 

' N+l 


,( 2 ) 

N 


as is easily shown by noting that if p G T^yv+i’ then for a certain fJ-N+i G 




( 2 ) 


( 2 ) 

Pn+i 



and of course p^li 


€ V S (Q N ). One can thus legitimately identify 



lim = 

N^-oo JV 


n v< n 

N>2 


(1.47) 


and the natural limit problem is then (modulo the justification of the above formal ma¬ 
nipulations) 

lim = Jinf { [ I H 2 (x, y)dp {2 \x, y), p {2) G v£> \ . (1.48) 

at^oo N 2 [JJaxn J 

This is a variational problem over the set of two-body states that one can obtain as reduced 
densities of IV-body states, for all N. 


A structure theorem. We now explain that actually 


E^ = E 


MF 


( 2 ) 

as a consequence of a fundamental result on the structure of the space Voo . 

Let us take a closer look at this space. Of course it contains the product states of the 
form p® p, p G V(yt) since p® 2 is the second marginal of p® N for all N . By convexity V ^ 
also contains all the convex combinations of product states: 

{/ p® 2 dP(p),PeV(V(n))\ cv£> (1.49) 

[Jpevm I 


with V(V(£l)) the set of probability measures over V(Q). In view of (11.121) and (|1.45l) we 
will have justified the mean-field approximation (11.441) if one can show that the infimum 
in (11. 48 p is attained for p^ = p® 2 for a certain p G V(£l). 

The structure result that allows us to reach this conclusion is the observation that there 
is in fact equality in (11.491) : 


[ P ® 2 dP(p),PeV(V(n)) \ =vgl 

Jp&v(n) I 


(1.50) 
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Indeed, in view of the the linearity of the energy functional as a function of one may 
write 

E 00 = Uni\[ ff H 2 (x,y)dp 92 (x,y)dP(p),P €V(Vm 
z { Jp&v(p.) J Jnxn 

= mf If £ MF [p]dP(p),P€V(Vm 
[Jpev( n) 

_ ^tMF 


since it is clear that the infimum over P £ V(V(£l)) is attained for P = 5 q mf, a Dirac 
mass at £ MF , a minimizer of £ MF . 

We thus see that, if we accept some formal manipulations (that we will justify in Chap¬ 
ter [2]) to arrive at (11.48p . the validity of the mean-field approximation follows by using very 
little of the properties of the Hamiltonian, but a lot the structure of symmetric IV-body 
states, in the form of the equality (11.501) . 

The latter is a consequence of the Hewitt-Savage, or classical de Finetti, theorem [ 501 
ED [93], recalled in Section ELD and proved in Section 12.21 Let us give a few details for the 
familiarized reader. We have to show the reverse inclusion in (11.49p . We thus pick some 
pS 2 ’ satisfying 


k <2) = PS 


for a certain sequence £ V S (Q N ). Recalling (11.471) we may assume that 


(iV) 

Pn+i = Pn 


and we thus have to deal with a sequence (hierarchy) of consistant IV-body states. There 
then exists (Kolmogorov’s extension theorem) a symmetric probability measure over the 
sequences of Q, p, £ 'P S (H N ) such that 


Pn = P 


(TV) 


where the IV-th marginal is defined as in (11.111) . The Hewitt-Savage theorem [93] then 
ensures the existence of a unique probability measure P £ V(V(£l)) such that 

Pn= [ P® N dP{p). 

JpeV(n) 

We obtain the desired result by taking the second marginal. 

Chapter [2] contains the details of the above proof scheme. We will specify the adequate 
assumptions to put all these considerations on a rigorours basis. It is however worth to 
note immediately that this proof (inspired from [ 146 . 951 [991 [ST , 101] ) used no structure 
property of the Hamiltonian (e.g. the interactions can be attractive or repulsive or a 
mixture of both) but only compactness and regularity properties. 


1.4. Bose-Einstein condensation and the quantum de Finetti theorem. 

We now sketch a strategy for justifying the mean-field approximation at the level of 
the ground state energy of a large bosonic system. The method is similar to that for the 
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classical case presented above. We consider a mean-field regime by setting A = (N — 1) 1 
and thus work with the Hamiltonian 

N 1 

h n = ^-{Vj + iA{x j )) 2 + V{xj) + —— w{xi-Xj) (1.51) 

j=l 1 l<i<j<N 


acting on L 2 (R dN ). Compared to the previous situation, we have added a vector potential 
A to start emphasizing the generality of the approach. It can for example correspond to 
an external magnetic field B = curl A. 

The starting point is the bosonic ground state energy defined as previously 

E(N) = inf {Ti-^tv) [H n T n ] , V N € S(L 2 (R dN ))} 

= mi[(^ N ,H N ^ N ),^ N eL 2 s (R dN ),\\^ N \\ L2{RdN) = l} (1.52) 


where we recall that we can minimize over pure or mixed states indifferently. 

Our goal is to show that for large N the bosonic energy per particle can be calculated 
using Hartree’s functional 

(1.53) 

where 


r E(N) 

lim - = eiT 

N—> oo N 


e H = inf [u],u € L 2 (R d ), ||w||= 1 } 

£k[u\ = f (|(V + iA)u \ 2 + V\u\ 2 ) + \ ff 

J R d * J jR d x\ 


u(x)\ 2 w(x 


y)\u(y)\ 2 dxdy. 


(1.54) 


Since Hartree’s energy is obtained by inserting an anzatz Tat = u® N in the IV-body 
energy functional, the asymptotic result (11.531) is already a strong indication in favor of 
the existence of BEC in the ground state of bosonic system, in the mean-field regime. We 
will come back later to the consequences of (11.531) for minimizers. As in the classical case, 
the mean-field regime is a very simplified, but already very instructive, model. We will 
present in Chapter [7] an anlysis of other physically relevant regimes. 


Formally passing to the limit. The first step to obtain (11.531) is as in the classical case to 
reduce formally to a simplified limit problem. We start by rewriting the energy using (|1.35l) 
and the assumption A = (N — l) -1 


E(N) 

N 



{T r L 2 (R 2d ) 


f 2 t ( 2) ,r^€74 2) } 


(1.55) 


where 

V {2) = {r( 2 ) € S(L 2 s (R 2d ))\3T N € S(L 2 s (R dN )),T& = Tr 3 ^v[rjv]} 

is the set of ” A r -representable” two-body density matrices that are the partial trace of a 
IV-body state. 

Characterizing the set V y N is a famous problem in quantum mechanics, see e.g. |43U124j . 
and the formulation (I1.55P is thus not particularly useful at fixed N. However, as in the 
classical case, the representability problem can be given a satisfactory answer in the limit 
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N —>■ oo, and we shall rely on this fact. We start by noticing that, taking partial traces, 
we easily obtain 

?%+! c vf (1.56) 

so that (I1.55P can be seen as a variational problem set in a variational set that gets more 
and more constrained when N increases. 

Assuming again that one can exchange infimum and limit (which is of course purely 
formal) we obtain 


TV— >• oo N 


inf jTr L 2 ( R 2 d } \h 2 r< 2) ] ,r® € (1.57) 


with 


^- l~K 2) 

N>2 


(1.58) 


the set of two-body density matrices that are ./V-representable for all N. As previously we 
have used the fact that the energy functional does not depend on N to pass to the limit 
formally. 

( 2 ) 

A structure theorem. It so happens that the structure of the set Vdo is entirely known 
and implies the equality 

Eoo = e H , (1.59) 

which concludes the proof of (|1.53l) . up to the justification of the formal manipulations 
we have just performed. 

The structure property leading to (11.591) is a quantum version of the de Finetti-Hewitt- 
Savage theorem mentioned in the prevous section. Pick a € V& . We then have a 
sequence Tat € 5(L 2 (M dAr )) such that for all N 

r( 2 ) = r^ } . 

Without loss of generality one can assume that this sequence is consistent in the sense 
that 

^ n+i = Triv+i[riv+i] = Tat. 

The quantum de Finetti theorem of Stprmer-Hudson-Moody 1192i (94| then guarantees the 
existence of a probability measure P € V(SL 2 (R d )) on the unit sphere of L 2 (M d ) such 
that 

r N = f \u® N )(u® N \dP(u) 

Ju£SL 2 (R d ) 

and thus 

= [ \u® 2 )(u® 2 \dP(u). 

Jui 

We can then conclude that 


lu<ESL 2 (JBL d ) 


Eca = inf • 


= inf ■ 


1 


Tr LHR 2 d) [H 2 \u® 2 )(u® 2 \}dP(u),P <E P(SL 2 (R d )) 


2 JueSL 2 (9. d ) 

[ £u[u]dP(u), P € P(SL 2 (R d )) 

JueSL 2 (R d ) 













26 


NICOLAS ROUGERIE 


where the last equality holds because it is clearly optimal to pick P = <5 Uh with uh a 
minimizer of Hartree’s functional. 

We thus see that the validity of (11.531) is (at least formally) a consequence of the 
structure of the set of bosonic states and only marginally depends on properties of the 
Hamiltonian (11.511) . The justification (under some assumptions) of the formal manipula¬ 
tions we just performed to arrive at (11.571) (or a variant), as well as the proof the quantum 
de Finetti theorem (along with generalizations and corollaries) are the main purpose of 
these notes. 


Bose-Einstein condensation. Let us anticipate a bit on the implications of (11.531) . We 
shall see in the sequel that results of the sort 

^ [ \u® n )(u® n \dP{u) (1.60) 

Ju&SL 2 {R d ) 

for all fixed re € N when N —>• oo, follow very naturally in good cases. Here is 
(the density matrix of) a minimizer of the IV-body energy and P a probability measure 
concentrated on minimizers of £h- The convergence will take place (still in good cases) in 
trace-class norm. 

When there is a unique (up to a constant phase) minimizer uh for £h we thus get 

|re^ n )(re^ ?l | when IV —»• oo, (1.61) 

which proves BEC at the level of the ground state. Indeed, BEC means by definition 
(see (125) and references therein, in particular (151) 1 the existence of an eigenvalue of 
ordeiO 1 in the limit IV —>• oo in the spectrum of T)J . This is clearly implied by (11.611) . 
which is in fact stronger. 

One can certainly wonder whether stronger results than (11.611) may be obtained. One 
could for example think of an approximation in norm like 


Tat — u 


(g> N 

H 


—> 0 . 

L 2 (R dN ) 


Let us mention immediately that results of this kind are wrong in general. The good 
condensation notion is formulated using density matrices, as the following two remarks 
show: 


• Let us think of a of the form (<8> s stands for the symmetric tensor product) 

<g> s p 

where p is orthogonal to ua . Such a state is “almost condensed” since all particles 
but one are in the state ren- However, in the usual L 2 (M. dN ) sense, is orthogonal 
to The two states thus cannot be close in norm, although their re-body density 
matrices for re <C N will be close. 

• Following the same line of ideas, it is natural to look for corrections to the IV-body 
minimizer under the form 

'LlV — ‘PO u n + 1% Pi + P2 T ■ ■ ■ T pn 

with po E C and pk € L^(M. dk ) for k = 1... N. It so happens that the above ansatz 
is correct if the sequence {pk)k=o,... N is chosen to minimize a certain effective 


6 


Let us recall that all our density matrices are normalized to have trace 1 in this course. 
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Hamiltonian on Fock space. The occurence of non-condensed terms implying the 
<Pk for k > 1 contributes to leading order to the norm of 'kw, but not to the reduced 
density matrices, which rigorously confirms that the correct notion of condensation 
is necessarily based on reduced density matrices. This remark only scratches the 
surface of a beautiful subject that we will not discuss here: Bogoliubov’s theory. 
We refer the interested reader to HUJ 113(11113111187111761155111171115(11155] for recent 
mathematical results in this direction. 

A remark concerning symmetry. We have been focusing on the bosonic problem because of 
our main physical motivations. The fermionic problem is (at least at present) not covered 
by such considerations, but one might be interested in the problem without symmetry 
constraint mentioned previously. 

In the case where A = 0 in (11.511) . the problems with and without bosonic symmetry 
coincide, but it is not the case in general. One may nevertheless use the same method 
as above to study the problem without symmetry constraint. Indeed, because of the 
symmetry of the Hamiltonian, one may without loss of generality impose the weaker 
symmetry notion (11.231) in the minimization problem. The set of two-body density matrices 
appearing in the limit is also covered by the Stprmer-Hudson-Moody theorem, and one 
then has to minimize an energy amongst two-body density matrices of the form 

T (2) = f 7 ® 2 dP(a) (1.62) 

J7G<S(L 2 (R d )) 

where P is now a probability measure on mixed one-particle states, that is on positive self- 
adjoint operators over L 2 (R d ) having trace 1. One deduces in this case the convergence of 
the ground state energy per particle to the minimum of the generalized Hartree functional 
(cf (H3DJ) 

£h[ 7] = Tr A 2 (R d ) [(“( V + iA ) 2 + V ) l\ + A jV 9 Tr L2{R 2 d) [w(x - y)7® 2 ] 
and the minimum is in general different from the minimum of Hartree’s functional (| 1.541) 

(see e.g. urn 

We already recalled that y ® 2 can have bosonic symmetry only if 7 is a pure state 
7 = |«)(«|. it is thus clear why the problem without symmetry can in general lead to a 
strictly lower energy. In fact the minimizer of (11.401) when there is no magnetic held is 
always attained at a pure state, in coherence with the different observations we already 
made on symmetry issues. We finally note that taking bosonic symmetry into account 
in the mean-held limit is done completely naturally using the different versions of the 
quantum de Finetti theorem. 
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2. Equilibrium statistical mechanics 

In this chapter we shall prove the classical de Finetti theorem informally discussed above 
and present applications to problems in classical mechanics. 

To simplify the exposition and with a view to applicaions, we will consider particles 
living in a domain Q cM. d which could be W 1 itself. We denote Q N and the cartesian 
product of N copies of and the set of sequences of Q respectively. The space of probability 
measures over a set A will always be denoted T’(A). In places we may make the simplifying 
assumption that 0 is compact, in which case V(Ll) is compact for the weak convergence 
of measures. 

2.1. The Hewitt-Savage Theorem. 

We mention only in passing the first works on what is now known as de Finetti’s 
theorem [50] mum Eg. In this course we shall start from [93] where the classical de 
Finetti theorem is proved in its most general form. Our point of view on the de Finetti 
theorem is here resolutely analytic in that we shall only deal with sequences of probability 
measures. More probabilistic versions of the theorem exist. We refer the reader to [2]|96] 
and references therein for developments in this direction. 

Informally, the Hewitt-Savage theorem [93] says that every symmetric probability mea¬ 
sure over Ll N approaches a convex combination of product probability measures when N 
gets large. A symmetric probability measures p N has to satisfy 

AHvC^i x ... x An) = /Tv(At(i),..., A a ^) (2-1) 

for every borelian domains A\,... ,An C II and every permutation of N indices a £ 

We denote V S (Q N ) the set of such probability measures. A product measure built on 
p € V(yt) is of the form 

p® N (A 1 ,...,A N ) = p(A 1 )...p(A N ) (2.2) 

and is of course symmetric. We are thus looking for a result looking like 

/Tv ~ [ P m dP tlN {p) when N ->• oo (2.3) 

JpeV(n) 

where Pp, N € V(V(Ll)) is a probability measure over probability measures. 

A first possibility for giving a meaning to m consists in taking immediately N = oo, 
that is, consider a probability with infinitely many variables p € V(Q N ) instead of p N , 
which lives on Cl N . This is the natural meaning one should give to a “classical state 
of a system with infinitely many particles”. We assume a notion of symmetry inherited 
from (| 2.1 1) : 

V>(Ai,A 2 , • • •) = MArli), Ar(2), • • •) (2-4) 

for every sequence of borelian domains (Afc)fc g M C and every permutation of infin¬ 
itely many indices a € Eqo. We denote P S (Q^) the set of probability measures over 
satisfying m- Hewitt and Savage [93] proved the following: 

Theorem 2.1 (Hewitt-Savage 1955). 

Let p € V s (Lt^) satisfy (12.41) . Let p^ be its n-th marginal, 

P {n) (A u ... , A n ) = p(A 1 ,...,A n ,Q,...,n,...). 


(2.5) 
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There exists a unique probability measure P ^ G V(V(£l)) such that 



( 2 . 6 ) 


In statistical mechanics, this theorem is applied as a weak version of the formal ap¬ 
proximation 112.311 in the following manner. We start from a classical iV-particle state 
whose marginals 

...,A n ) = ...,A n , Q N ~ n ) (2.7) 

converg^ as measures in V(Q n ), up to a (non-relabeled) subsequence: 

^ p {n) € v{^ n) ). ( 2 . 8 ) 

This means that 

^\A n )^^\A n ) 

for all borelian subset A n of fl n , and thus 

[ fndfX^ [ f n dH (n) 

Jft Jq 

for all bounded continuous functions decaying at infinity f n from to M. Modulo a 
diagonal extraction argument one may always assume that the convergence (|2.8I) is along 
the same subsequence for any n € N. Testing (12.81) with a Borelian A n = A m x $V n ~ n for 
m < n one obtains the consistency relation 

^( n )^ = pS- m \ for all m < n (2-9) 


which implies that (p^ n ')neN does describe a system with infinitely many particles. One 
may then see (Kolmogorov’s extension theorem) that there exists p € (P(fl N ) such that 
is its n-th marginal (whence the notation). This measure satisfies (12.41) and we can 
thus apply Theorem 12.II to obtain 



when N —>• oo 


( 2 . 10 ) 


where G V(V{£1)). This is a weak but rigorous version of (12.31) . In words the ra¬ 
th marginal of a classical IV-particle state may be approximated by a convex 
combination of product states when N gets large and n is fixed . Note that the 
measure P ^ appearing in (12.101) does not depend on n. 

Of course, one must first be able to use a compactness argument in order to obtain (12.81) . 
This is possible if Q is compact (then (P(f l n ) is compact for the convergence in the sense of 
measures (ESD). More generally, if the physical problem we are interested in has a confining 
mechanism, one may show that the marginals of equilibrium states of the system are tight 
and deduce (|2.8D . 


As for the proof of Theorem 12.11 there are several possible approaches. We mention 
first that the uniqueness part is a rather simple consequence of a density argument in 
C'b('P(n)), the bounded continuous functions over V(£l) due to Pierre-Louis Lions [ 14C1 . 


^We denote this convergence 


to distinguish it from a norm convergence. 
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Proof of Theorem, \ 2. 11 Uniqueness. We easily verify that monomials of the form 

C b (V(n)) 3 M k j(p) := j <j>(x 1 ,...,x k )dp 9k (x 1 ,...,x k ),k€N,<j>€C b (n k ) (2.11) 

generate a dense sub-algebra of C b (fP(Q)), the space of bounded continuous functions over 
V(£l), see Section 1.7.3 in [82] . The main point is that for all p G V{Ui) 

^k+i,cf> <S> • 

That this sub-algbra is dense is a consequence of the Stone-Weierstrass theorem. 

We thus need only check that if there exists two measures and P ^ satisfying (12.51) . 
then 

[ MkMdP^p) = [ MkMdP'^p) 

Jpev(n) Jp&v(n) 

for all k G N and (j) G C b (V(Q k )). But this last equation simply means that 


/ (/ 

Jp&v(P) \Jn k 


4>(x 1 ,...,x k )dp® k (xi,...,x k ) dP /J (p) 


I (/ 


4>(x 1 ,...,x k )dp (Z)k (x 1 ,...,x k ) ) dP'(p) 


lp&V(n) \Jn k 

which is obvious since both expressions are equal to 


/ (f(xi,...,x k )dpi { ' k \x 1 ,...,x k ) 

Jn k 


by assumption. 


□ 


For the existence of the measure, which is the most remarkable point, we mention three 
approaches 

• The original proof of Hewitt-Savage is geometrical: the set of symmetric probability 

measures over is of course convex. The Choquet-Krein-Milman theorem says 
that any point of a convex set is a convex combination of the extremal points. It 
thus suffices to show that the extremal points of coincide with product 

measures, corresponding to sequence of marginals of the form (p® n )neN- This 
approach is not constructive, and the proof that the sequences (/9® n ) nS M are the 
extremal points of is by contradiction. 

• An entirely constructive approach is due to Diaconis-Freedman m ■ In this prob¬ 
abilistic argument, Theorem 12. II is a corollary of an approximation result at finite 
N, giving a quantitative version of (12.101) . 

• Lions developed a new approach for the needs of mean-field game theory [liOj . 
This is a dual point of view where one starts from “weakly dependent continuous 
functions” with many variables. A summary of this is in the lecture notes [82 , 
Section 1.7.3], and a thorough presentation in [ 147 , Chapter I], Developments and 
generalizations may be found in [91|. 

It so happens that the proof of the Hewitt-Savage theorem following Lions’ point of view 
is for a large part a rediscovery of the method of Diaconis and Freedman. In the sequel 
we will follow a blend of the two appraoches. See [147] for a more complete discussion. 
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2.2. The Diaconis-Freedman theorem. 

As we just announced, it is in fact possible to give a quantitative version of (J2.10I) 
which implies Theorem 12.11 Apart from its intrinsic interest this result naturally leads to 
a constructive proof of the measure existence. The approximation will be quantified in 
the natural norm for probability measures over a set S, the total variation norm: 



sup 

06 C b (S) 


/ (jxlfJL 

Jn 


( 2 . 12 ) 


which coincides with the L 1 for absolutely continuous measures. We shall prove the 
following result, taken from m- 


Theorem 2.2 (Diaconis-Freedman). 

Let fi N € V s (Ll N ) be a symmetric probability measure. There exists P tlN € V(V(Ll)) such 
that, setting 

Pn-= [ P® N dP» N {p) (2.13) 

■JpeT(n) 


we have 




n(n — 1) 

< 2 —-- 

TV N 


(2.14) 


Proof. We slightly abuse notation by writing p, N (Z)dZ instead of dfi N (Z) for integrals in 
(z ±,..., zn) = Z e Ll N . It is anyway already instructive enough to consider the case of an 
absolutely continuous measure. 

By symmetry of p, N we have, for all X = (x \,..., xn) € f l N , 

» N (X) = [ Hn(Z) Y m~ 1 Sx=z a dZ (2.15) 

JaN a^ N 

where Z a is the IV-tuple (z a m ,..., ^ CT ( n ))- We define 




MZ) Y N~ N 5 x =z,dZ , 

Terjv 


(2.16) 


where Tjy is the set of all application^ from {1,..., N} to itself and is defined is the 
same manner as X a . The precise meaning of (|2.16l) is 


/ 

Jn N 


<t>(X)di* N (X) = 


^ n ~ n L 


(j){Z^)dfi N (Z) 


Terjv 


for every regular function <f>. 
Noting that 


Y N N dx=Z 1 

7er N 


/ N y N 


(2.17) 


'Compared to Sjv, r jv thus allows repeated indices. 
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one can put (12.1611 under the form (12.1311 by taking 

p n N (p) = I 5 p=P zH N ( z )dZ, p z (x) :=£n-% 

J ° N i= i 


N 


(2.18) 


Note in passing that P Pn charges only empirical measures (of the form pz above). We 
now estimate the difference between the marginals of p, N and fi N . Diaconis and Freedman 
proceed as follows: of course 

/ / \ (n) / \ («)' 


(n) ~ (n) 

/% ~ /% = 




V 


E ( N -) 1 dx=z a 


but 


E n ~ NSx 
N 

(n) 


= Z~, 


I A 1 A r{Z)dZ, 


:=Z CT 


E ( m )~ l6x = 

ytrGSjv 

is the probability law for drawing n balls at random from an urn containing N ball^O, 
without replacement, whereas 

^ (n) 

E N~ n 5x=z 7 

is the probability law for drawing n balls at random from an urn containing N, with 
replacement. Intuitively it is clear that when n is small compared to N, the fact that we 
replace the balls or not after each drawing does not significantly influence the result. It is 
not difficult to obtain quantitative bounds leading to (|2.14j) . see for example [751 . 

Another way of obtaining (|2.14j) . which seems to originate in [87j, is as follows: in view 
of (12.161) and (12.171) . we have 






(xi,.. .,x N ) dZ. 


We then expand the tensor product and compute the contribution of terms where all 

(n) 

indices are different. By symmetry of i± N we obtain 


(n) _ N(N — 1) ... (IV — TO + 1) ( n ) 


P'N ~ 


N n 


I^N p Vn 


(2.19) 


where fn is a positive measure on Q n (all terms obtained by expanding (I2.17|) are positive). 
We thus have 

N(N — 1)... (N — n + 1)' 


(n) ~{n) I -i 

Mat — Mn — 1 _ 


N n 


( n ) 

V N ~ v n 


( 2 . 20 ) 


and since both terms in the left-hand side are probabilities, we deduce 

N(N — 1)... (IV — n + 1)' 


/ 

J 


du n = 1 — 


N n 


^Labeled zi,.. ., zn- 
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Moreover, since the first term of the right-hand side of (12.201) is positive and the second 
one negative we obtain by the triangle inequality 


d 


(n) ~ (n) 

/% - 


N(N - 1)... (N - n + 1) 


= 2 1 - 


It is then easy to see that 

N(N - 1)... (N - n + 1) 

m 


n 

3 =1 


N-j + 1 
N 


nb 

3 = 1 


j ~ 1 

N 


> 1 - 


n — 1 


iV 


> 1 - 


n(n — 1) 

iv : 


which proves (12.141) with C = 2. A better constant C = 1 can be obtained, see [57 , 751 . □ 

The following remark can be useful in applications: 


Remark 2.3 (First marginals of the Diaconis-Freedman measure). 
We have 


Aw (*) = Mw(z) 


and 


I*n\xi, x 2) = A N + ^w( x i)<^i=z 

as direct consequences of Definition (12.171) . Indeed, using symmetry, 

AT . 

-.(F/'.vO _ AT— 1 — ..i 1 )/ 


( 2 . 21 ) 

( 2 . 22 ) 


A;v ( x ) = N / ^ N {Z)d Zj=x dZ = n [ N \x) 


3 =1 




AwC^i^s) = AT 


-2 


(|r^"‘) (|r 


$Zj =X2 I 




M TV ( =#i $Zj =X2 


= ‘ v - 2 E 

iV 

+ iV ' 2 E/ Ma r(Z)S Zi=xl 6 Zi=X2 dZ 

= l.® 2 ) + 

Higher-order marginals can be obtained by similar but heavier computations. 


(2.23) 

□ 


As a corollary of the preceding theorem we obtain a simple proof of the existence part 
in the Hewitt-Savage theorem: 
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Proof of Theorem, \ 2. 11 Existence. We start with the case where Pi is compact. We apply 
Theorem 12.21 to obtaining 

M (n) -/ P m dP N (p ) 

JpeP{n) 

with Pjy € V{fP{Pl)) the measure defined in (12.1811 . When Pi is compact, V{Pl) and V{V{Pl)) 
also are. We thus may (up to a subsequence) assume that 

Pn^Pp P(vm 


TV 


- N 


(2.24) 


in the sense of measures and there only remains to pass to the limit IV -A oo at fixed n 
in (frail . 

When Pi is not compact, we follow an idea of Lions [140] (see also [82] )• We want to 
ensure that the measure Px obtained by applying the Diaconis-Freedman theorem to ^ N ' 1 
converges. For this it suffices to test against a monomial of the form (12.1111 : 

[ M k ^p)dP N (p) = f M ki * ( I d ^\ z ) 

Mvm Jz&n N \ N ~r[ ) 

= 1 [ 0(xi,...,x fc )n ( Tf ) dp, (N \ z ) 

= [ (f{z 1 ,...,z k )d^ k \z) + 0(N- 1 ) (2.25) 

Jzen k 

by a computation similar to that yielding (12.1911 . The limit (12.2511 thus exists for any 
monomial and by density of monomials for any bounded continuous function over 

V(Pl). We then deduce 

Pn * P 

and we can conclude as previously. □ 


Some remarks before going to the applications of Theorems 12.11 and 12.21 

Remark 2.4 (On the Diaconis-Freedman-Lions construction). 

(1) We first note that the measure defined by (12.1811 is that Lions uses in his approach 
of the Hewitt-Savage theorem. Using empirical measures is the canonical way to 
construct a measure over V{V{Pl)) given one over V s {Pl N ). Passing to the limit 
as in ([2.2511 can replace the explicit estimate (12.1411 if one is only interested in the 
proof of Theorem 12.11 This construction and the combinatorial trick (12.1911 are 
also used e.g. in [83, . 148. 9T| . 

(2) Having an explicit estimate of the form (12. 1411 at hand is very satisfying and can 
prove useful in applications. One might wonder whether the obtained convergence 
rate is optimal. Perhaps suprisingly it is. One might have expected a useful 
estimate for n <S N, but it happens that y/n -C N is optimal, see examples in m- 

(3) The formulae ([2.2311 are useful in practice (see | T70 j for an application). It is quite 
satisfying that [iff = jiff and that faff can be reconstructed using only fi Z> 
and ■ 
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(4) As a drawback of its generality, the previous construction actually behaves very 
badly in many cases. Note that (12.181) charges only empirical measures, which 
all have infinite entropy. This causes problems when employing Theorem 12.21 to 
the study of a functional with temperature. Moreover, in a situation with strong 
repulsive interactions, one typically applies the construction to a measure with 
zero probability of having two particles at the same place, p, y N (x,x) = 0. In this 
case fi,ff(x,x) = N~ x fj,^\x) is non-zero and the energy of fi N will be infinite if 
the interaction potential has a singularity at the origin. 

(5) It would be very interesting to have a construction leading to an estimate of the 
form (12.141) . bypassing the aforementioned inconveniences. For example, is it pos¬ 
sible to guarantee that jl N € L 1 (Q N ) if p, N € L 1 (Q N ) ? One might also demand 
that the construction leave product measures p® N invariant, which is not at all 
the case for (12.131) . 

(6) If fl is replaced by a finite set, say II = {l,...,d}, one may obtain an error 
proportional to dn/N instead of n 2 /N by using the original proof of Diaconis- 
Freedman. One may thus replace (12.141) by 


(n) 

/% 


~ (n) 

f*N 


c 

< — mm 
TV N 


(dn , n 2 ) 


(2.26) 


We will not use this point anywhere in the sequel. 


□ 


Finally, let us make a separate remark on a possible generalization of the approach 
above: 


Remark 2.5 (Weakly dependent symmetric functions of many variables.). 

In applications (see next section), one is lead to apply de Finetti-like theorems in the 
following manner: given a sequence (u/v)veN of symmetric functions of N variables, we 
study the quantity 

f u N (X)dn N (X) 

Jn N 

for a symmetric probability measure p, N . The results previously discussed imply that if 

in) 

this happens to depend only on a marginal fx y N for fixed n, a natual limit object appears 
when N —>• oo. In particular, if 

u n(X)=( N \ 0(®ii,...,Xi B ), (2.27) 

' l<ii<...<i n <N 

we have 

f u N {X)dfjL N {X) = [ cj)(xi,..., x n )dn$ [ ( [ <t> dp®™) dPfj,(p) (2.28) 

Jn N ID" Jpev(n) \Jn n / 

for a certain probability measure € P(P(fI)). One might ask whether this kind of 
results is true for a class of functions un depending on N in more subtle a manner. The 
natural assumptions seems to be that un depends weakly on its N variables , in the sense 
introduced by Lions }l40j and recalled in | i82l Section 1.7.3]. Without entering the details, 
one may easily see that to such a sequence there corresponds (modulo extraction of a 
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subsequence) a continuous function over probabilities U € C(V(£l)): one may show that, 
along a subsequence, 


[ u N {X)dv N {X ) 
Jn N 


L 


pev( n) 


U(p)dP fl (p). 


(2.29) 


Think of the case where ujy depends only on the empirical measure: 


un(x i, • • -,x N ) = F 



(2.30) 


with a sufficiently regular function F. Such a function depends weakly on its N variables 
in the sense of Lions, but cannot be written in the form (12.2711 . More generally, this kind 
of considerations can be applied under assumptions of the kind 

(j 

\V Xj u N (xi,...,x N )\ < — Vj = 1... N, V(xi,... ,x N ) € Q N . 

One may then wonder whether the convergence rate in (12.2911 can be quantified. Results 
in this direction may be found in [91] . □ 


2.3. Mean-field limit for a classical free-energy functional. 

In this section we apply Theorem 12. II to the study of a free-energy functional at positive 
temperature, following [1461131119911101] . We consider a domain 12 C M. d and the functional 

F N [n\= [ H N (X)p{X)dX + T I /x(X)log n(X)dX (2.31) 

J xm N Jn N 

defined for probability measures /x e Here the temperature T will be fixed in the 

limit IV —>• oo and the Hamiltonian is in mean-held scaling: 

N 1 

H N {X) = Y J V{x j ) + w — [ Y, ( 2 - 32 ) 

j=1 1 

We denote by V a lower semi-continuous potential. To be in a compact setting we shall 
assume that either 12 is bounded or 


V(x) —>• oo when |x| —>• oo. (2.33) 

The interaction potential w will be bounded below and lower semi-continuous. To be 
concrete, one may think of w € L°°, or a repulsive Coulomb potential: 


w(x) = 

^ if <7 "I 

(2.34) 

. . , 0 11 Ul — O 

X ^ " 

w(x) = 

— log |x| if d = 2 

(2.35) 

w(x) = 

— |x| if d = 1, 

(2.36) 


ubiquitous in applications. We will always take w even, 

w(x) = w(—x), 

and if the domain is not bounded we will assume 


w(x — y) + V (x) + V(y) —>• oo when |x| -A oo or |y| —>• oo 
w(x — y) + V(x) + V ( y ) is lower semi-continuous. 


(2.37) 
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We are interested in the limit of the Gibbs measure minimizing f!2.31li in "P s (f ) N ): 

V N (X) = exp (\ dX 


'N 


and to the corresponding free-energy 


F n = inf 7i v [n] = FnIVn} = ~TlogZ N . 
U G7W) 


(2.38) 

(2.39) 


The free-energy functional is rewritten 

X N [fi\ = N [ V{x)dn [l \x) + ^ I f w(x - y)dp, {2 \x,y)+ T f plog 
Jn 1 JJnxn Jn N 


P 


_ N 

~ t 

using the marginals 


{w(x - y) + V(x) + V(y)) d^ 2 \x, y) + T / /rlog/x. (2.40) 
nxo Jn N 


/z (n) (xi, ... ,x n ) = / dfj.(xi,...,x N ). 

^ ^n+lv) 


(2.41) 

(2.42) 


Inserting an ansatz of the form 

v = p® N , P €V(n) 

in (12.3111 we obtain the mean-field functional 

F MF [p] := N-'XnIp®^ 

= f I f{x)dp{x) + \ ff w(x — y)dp(x)dp(y) + T [ plogp (2.43) 
JQ ^ J JQxQ J 


with minimum F mf and minimizer (not necessarily unique) £ MF amongst probability 
measures. Our goal is to justify the mean-field approximation by proving the following 
theorem: 


Theorem 2.6 (Mean-field limit at fixed temperature). 

We have 

—— —> F mf when N —>• oo. 

N 

Moreover, up to a subsequence, we have for every n € N 


(n) 

p n 


/ 

J p&AP 


p® n dP{p). 


as measures where P is a probability measurer over A4 MF , the set of minimizers of . 


(2.44) 


(2.45) 


rMF 


In particular, if has a unique minimizer we obtain, for the whole sequence, 


(n) 

Pn 


(^ MF ) lg)n 


Proof. We follow [ 4461 131, 98] [99] - An upper bound to the free energy is eeasily obtained 
using test functions of the form p® N and we deduce 

Fn 


N 


< F 


MF 


(2.46) 
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The corresponding lower bound requires more work. We start by extracting a subsequence 
along which 

^ /^ (n) (2-47) 

for all n £ N, with fi £ V S (Q N ). This is done as explained in Section 12.11 using either 
the fact that f l is compact or Assumption (12.371) . which implies that the sequence is tight. 
Indeed, using simple energy upper and lower bounds one easily sees that for all e > 0 
there exists R e such that 

(B( 0, R e ) x B{ 0, R e )) > 1 - e 

for all N. Tightness of all the other marginals follows. 

By lower semi-continuity we immediately have 


hminfi ff (w{x - y) + V(x) + V(y)) d^\x, y) 

N^oo 2 JJ nxn 

-\[[ ( w{x-y) + V(x)+V(y))dn {2] (x,y ). (2.48) 

^ J ifixSI 


For the entropy term we use the sub-additivity property (this is a consequence of Jensen’s 
inequality, see or the previously cited references) 


ln N 


/Tv lo gAbv > 


N 

n 


In™ 


/*!?' /»$?* + 


. l#J A N -"^ j) iog„r”^) 


where [_ • J denotes the integer part. Jensen’s inequality implies that for probability mea¬ 
sures n and z/, the relative entropy of /x with respect to u is positive: 



v— log — > 

V V 


= 0. 


We deduce that for all vq £ 'P(fl) 





Hk —'log n *-"^ log ( -W7 




®(N-n LfJ) 


+ 


L 


'-"LnJ 


R'N 


N 1(J & 


> 


^N- n 


N 

n 



logt'o- 


Choosing zx 0 £ 'P(fl) of the form zx 0 = coexp(—ciV) it is not difficult to see that the last 
integral is bounded below independently of N and we thus obtain that for all n £ N 

Urn inf [ n N log n N > - [ /x (n) log/x (n) (2.49) 

N xoo N J^n n Jqtl 

by lower semi-continuity of (minus) the entropy. 
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Gathering (I2.48f) and (12.4911 we obtain a lower bound in terms of a functional of /x: 

{w(x - y) + V(x) + V(y)) d/x (2) (x, y) 

+ T lim sup — f ^ n ' > log /x^. (2.50) 

n-¥oo Tl J<; 2" 


lim inf ^n^n] > Hu] ■= \ 

N^OO N 2 ././Q X Q 


The second term in the right-hand side is called (minus) the mean entropy of /x € lP(fl N ). 
We will next apply the Hewitt-Savage theorem to /x. The first term in T is obviously 
affine as a function of /x, which is perfect to apply (12.61) . One might however worry at the 
sight of the second term which rather looks convex. In fact a simple argument of [163] 
shows that this mean entropy is affine. It is a part of the statement that the lim sup is in 
fact a limit: 


Lemma 2.7 (Properties of the mean entropy). 

The functional 


1 / 

sup — /X 

n—Hxo n Jq 


(n) log /x (n) = lim - f /x^ log /x (n) 


IHOO Tl jQn 


(2.51) 


is affine over (P(fl N ). 


Proof. We start by proving that the sup equals the limit. Denote the former by S and 
pick j such that 


- [ log /*&') = sup - f n {n) log /x (n) — £j = S — 

j Jflj n-too Tl Jfin 


Next let k > j and write (Euclidean division) 

k = mj + r with 0 < r < j. 

Then, by the subbaditivity property already mentioned 

- [ ^ fe) l0g/X W >— f nWWnW) 
k Jnk k 


£n. 


/x UJ log /x u ' ) + j [ /x! r ) log/x 
M k Jnr 


(r) 


k 


k 


where Cj only depends on j. We then pass to the limit k —> oo at fixed j. Noting that 
then m —>• k/j we obtain 


lim inf — I /x^ log > S — £j. 

k^oo k J^k 

Passing finally to the limit j —>• oo along an appropriate subsequence we can make £j —>• 0, 
which proves the equality in (12.511) . 
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To prove that the functional is affine, let p 1; g 2 € 'P(fl N ) be given. We use the convexity 
of x eA x log x and the monotonicity of x eA log x to obtain 


If Mi 0 !°g n) + l[ /4 n)l °g/4 n) ^ / 

Z JO" 4 J^n J^n 


1 


m!' 0 + ^2 tJ ) lo g ( + X/4 ? 


1 


.0) 


.( n ) 


.(n) 


(n) 


>l[ Ml' 0 log Ml' 0 + l I V 2 ] log l4 

£ Jsi n ^ JQ™ 

- ^ ([ M n> + / d" 1 ) 

z Vif2™ J!) n / 

= r ) f P\ l) log Mi'° + 1/ /4 n) lo gM2 n) - lo g( 2 )- 
^ ^in™ 

Dividing by n and passing to the limit we deduce that the functional is indeed affine. □ 


Thus P[p] is affine. There remains to use (12.61) . which gives a probability P tl G V{V(p.)) 
such that 

hnnnfljWtgv]^ [ P[p®°°]dP„(p) 

N ^°° M Jpavipi) 

= [ .F MF [p]dP M (p) > F mf . 

Jp&V(P) 

Here we denote p®°° the probability over which has p® n for n-th marginal for all 

n and we have used the fact that P ^ has integral 1. This concludes the proof of (|2.44l) 
and (|2.45l) follows easily. It is indeed clear in view of the previous inequalities that P ^ 
must be concentrated on the set of minimizers of the mean-field free-energy functional. □ 


2.4. Quantitative estimates in the mean-field/small temperature limit. 

Here we give an example where the Diaconis-Freedman construction is useful to supple¬ 
ment the use of the Hewitt-Savage theorem. As mentioned in Remark l2.4l this construction 
behaves rather badly with respect to the entropy, but there is a fair number of interesting 
problems where it makes sense to consider a small temperature in the limit N —?> oo, in 
which case entropy plays a small role. 

A central example is that of log-gases. It is well-known that the distribution of eigen¬ 
values of certain random matrices ensembles is given by the Gibbs measure of a classical 
gas with logarithmic interactions. Moreover, it so happens that the relevant limit for 
large matrices is a mean-field regime with temperature of the order of IV -1 . Consider the 
following Hamiltonian (assumptions on V are as previously, taking w = — log | . |): 

N 1 

H N {X) = Y J V(x j )- w — [ Y, '°S\xi-Xj\ (2-52) 

j =1 1 <i<j<N 

where X = (aq,... ,xn) € W iN . The associated Gibbs measure 

P N (X) = ex P (~PNH n {X)) dX (2.53) 

corresponds (modulo a /3-dependent change of scale) to the distribution of of the eigenval¬ 
ues of a random matrix in the following cases: 
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• d = l,/3 = 1,2,4 and V(x) = ■^j-. We respectively obtain the gaussian real 
symmetric, complex hermitian and quarternionic self-dual matrices. 

• d = 2,/3 = 2 and V(x) = We then obtain the gaussian matrices without 
symmetry conditon, the so-called Ginibre ensemble [78] , 

In these notes we give no further precisions on the random matrix aspect, and will simply 
take the previous facts as a sufficient motivation to study the N —>• oo limit of the mea¬ 
sures (12.531) with /3 fixed. This corresponds (compare with (12.381) 1 to taking T = [3~ 1 N~ l , 
i.e. a very small temperature. For an introduction to random matrices and log-gases we 
refer to [8, 73], 1145] . For precise studies of the measures (12.531) following different methods 
than what we shall do here, one may consult e.g. pH Q21 EH [201 ED EQl S3 EH HZ3 HT2]. 

In the case d = 2, measures of the form (12.531) also have a natural application to the 
study of certain quantum wave-functions appearing in the study of the fractional quantum 
Hall effect (see [ 1681 169. 1701 1171 ] and references therein). Here too it is sensible to 
consider (3 as being fixed. 

The singularity at the origin of the logarithm poses difficulties in the proof, as indi¬ 
cated in Remark 12.41 but one may bypass them easily, contrarily to those linked to the 
entropy. The following method is not limited to log-gases and can be re-employed in 
various contexts. 


Again, (12.531) minimizes a free-energy functional 

Wv[R] = [ H N {X)ti(X)dX + -±- [ fi(X)log v(X)dX (2.54) 

Jxea N pA Jqjv 

whose minimum we denote lyy. The natural limit object is this time an energy functional 
with no entropy term: 

£ MF [p] ■= [ Vdp-\[[ log\x-y\dp(x)dp(y), (2.55) 

JR d * J J R d xR d 


obtained by inserting the ansatz p® N in (12.541) and neglecting the entropy term, which is 
manifestly of lower order for fixed /3. We denote E MF and g MF respectively the minimum 
energy and the minimizer (unique in this case by strict convexity of the functional). It is 
well-known (see the previous references as well as [ 101] ) that 


N~'F n 


~\og Z n ^E mf 


when N —>• oo 


(2.56) 


and 


We shall prove (12.571) and give 
from [ 170 ]: 


(2.57) 

a quantitative version of (12.561) . taking inspiration 


Theorem 2.8 (Free-energy estimate for a log-gas). 

For all f3 € M, we have 

E mf — CN -1 (/3 -1 + log IV + l) <~\og Z N < E MF ECP^N- 1 . (2.58) 

Fine estimates for the partition function of the log-gase seem to have become avail¬ 
able only recently [TTMITH81HH7HT78HT72] . Amongst other things, the previous references 
indicate that the correction to E MF is exactly of order N log N. 
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Proof. We shall not elaborate on the upper bound, which is easily obtained by taking 
the usual product ansatz and estimating the entropy. For the lower bound we shall use 
Theorem 12.21 We first need a crude bound on the entropy term: by positivity of the 
relative entropy (via Jensen’s inequality) 


/ log — > 0 for all ii,iy £ V{Pl 

J R dN u 

one may write, using the probability measure 

u N = (c 0 exp(-F(|x|)))® 7V 

the lower bound 


N\ 


IR dN 


P-N log P N > 


L 


fi N logu N = -N I Vd/iff — N log cq ■ 


(2.59) 


To obtain a lower bound to the energy we first need to regularize the interaction potential: 
Let a > 0 be a small parameter to be optimized over later and 

-loga + \ (l - if \z\ < a. 

v — log \z\ if \z\ > a. 

Clearly — log Q \z\ < — log \z\ is regular at the origin. Moreover, we have 


-logaM = 


(2.60) 


d i . , 
-~r^g a \z\ = 
da 


_i Hr 

a a 3 


si \z\ < a 
si \z I > a. 


(2.61) 


Using the lower bound — log Q \z\ < — log \z\ to obtain 


[ H N (X)ii(X)dX >N f Vdp$ -^ [[ 

JXgQ n J z J JR d xR d 

we are now in a position where we can apply Theorem 12.21 
explicit formulae (12.2311 : 


log a \x - y\d^\x,y) 

More precisely we use the 


/ H N (X)ii N (X)dX > 

Jxen N 

N [ Vdji$ - ^ [f logo \ x ~ y\dp^{x,y) + CTog a (0). (2.62) 

JR d z \ lsl ~ l) J JWL d xR d 

Combining (12.59]) , (12.1311 and recalling that the temperature T is equal to we 

obtain 

> [ £T\p\dP„M 
Jp&v( R d ) 

+ CN- 1 (log Q (0) - /r 1 ) > Ea F - CN~ l log (a) - C(PN)- 1 (2.63) 

where is the minimum (amongst probability measures) of the modified functional 

£™ F [p\-= [ U(1 ~ P~ l N~ l )dp- ^ [[ log a \x-y\dp(x)dp(y). 

jR d 2(iV - 1) J J K d xR d 
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Exploiting the associated variational equation, (12.611) and the Feynman-Hellmann princi¬ 
ple, it is not difficult to see that for small enough a 

\ E MF _ £MF| < Ca d + CN -1 

and we thus conclude 

N- ] F n [h n ] > E mf - CN- 1 log(a) - CfiN- 1 - CIV" 1 - Ca d - C , 
which gives the desired lower bound after optimizing over a (take a = N~ l / d ). □ 

Remark 2.9 (Possible extensions). 

(1) One can also prove quantitative versions of (12.571) following essentially the above 
method of proof. We shall not elaborate on this point for which we refer to 170l . 

(2) Another case we could deal with along the same lines is that of unitary, orthogonal 

and symplectic gaussian random matrices ensembles introduced by Dyson [601 [61 . 
[62] , In this case M. d is replaced by the unit circle, /3 = 1, 2,4, V = 0 in (12.521) . In 
this case, dividing the interaction by N — 1 is irrelevant. □ 
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3. The quantum de Finetti theorem and Hartree’s theory 

Now we enter the heart of the matter, i.e. mean-field limits for large bosonic sys¬ 
tems. We present first the derivation of the ground state of Hartree’s theory for confined 
particles, e.g. living in a bounded domain. In such a case, the result is a rather straight¬ 
forward consequence of the quantum de Finetti theorem proved by Stprmer and Hudson- 
Moody [1921 194] , The latter describes all the strong limits (in the sense of the trace-class 
norm) of the reduced density matrices of a large bosonic system. 

We will then move on to the more complex case of non-conhned systems. In this 
chapter we will assume that the interaction potential has no bound state (it could be 
purely repulsive for example), the general case being dealt with later (Chapter [6]). In 
the absence of bound states it is sufficient to have at our disposal a de Finetti theorem 
describing all the weak limits (in the sense of the weak-* topology on S 1 ) of the reduced 
density matrices. 

The weak de Finetti theorem (introduced in [ 113] ) implies the strong de Finetti theorem 
and in fact the two results can be deduced from an even more general theorem appearing 
in [1~92[ [9il . In these notes I chose not to follow this approach, but rather that of 0 
m\ which is more constructive. This will be discussed in details in Section which 
announces the plan of the next chapters. 

3.1. Setting the stage. 

To simplify the discusstion, we will focus on the case of non relativistic quantum par¬ 
ticles, in the absence of a magnetic field. The Hamiltonian will thus have the general 
form 

N 1 

r n = Y1 T i + w ( Xi ~ ( 3 - 1 ) 

j=l 1 <i<j<N 

acting on the Hilbert space fji) = S), i.e. the symmetric tensor product of N copies of 

fj where Sj denotes the space L 2 (fi) for HcR'*. The operator T is a Schrodinger operator 

T = —A + V (3.2) 

with V : H H > R and T) acts on the j-th variable: 

Tjibi <g)... <g) ipN = Vh <8> • • • <8> Tjipj <g) ... <g) 

We assume that T is self-adjoint and bounded below, and that the interaction potential 
w : M i-a R is bounded relatively to T (as operators): for some 0 < /3_,/3_|_ < 1 

- (3-(Ti + T 2 )-C< w(x i - x 2 ) < /3+(Ti + T 2 ) + C. (3.3) 

We also take w symmetric 

w(—x) = w(x), 

and decaying at infinity 

w € L P (H) + L°°(n), max(l, d/2) < p < oo —>• 0, w(x) —> 0 when \x\ -A oo. (3.4) 

We will always make an abuse of notation by writing w for the multiplication operator 
by w(x\ — x 2 ) on L 2 (ki). 
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Remark 3.1 (Checking assumptions on w). 

One may check that assumptions (13.41) imply operator bounds of the form (13.31) . using 
standard functional inequalities. We quickly sketch the method here for completness. 
Suppose that w € L p (M. d ), such that the conjugate Holder exponent q is smaller than half 
the critical exponent for the Sobolev injection i4 1 (M rf ) ^A 

2d 


Q = 


P 


P 


p- 1<7’ P =d- 2' 


(3.5) 


We decompose |u;| as 


M = MlH<R + 

Then, pick / € C)? 0 (M 2rf ) and write 

|2 




< 


R d xR d 


w(x 1 ~ X 2 )\f(xi,X 2 )\ dx\dx 2 
w(x 1 - ^2)1|«,|>_r| \f{xi,x 2 )\ 2 dxidx 2 + R 


R d xR d 


\f(xi 1 x 2 )\ 2 dx l dx 2 . 


Next, using Holder’s and Sobolev’s inequalities: 


JL 


R d xR d 


\w{xi - a: 2 )l| w |>i?| \f(xi, x 2 )\ 2 dxidx 2 

' f ( [ \ w ( Xl ~ x i)\ Pl H>R dx 2 ) 

J R4 \J R4 / 


1/9 


\f{xi,x 2 )\ 2q dx 2 \ dx 1 

Id J 

< l\wt\ w \> R \\ [ Wf(.,x 2 )\\ 2 Hl dx 2 

JR d 

with p -1 + q~ 1 = 1. All in all we thus have 

\{f,uif)\<\\wt\ w \ >R \\ // \\7 Xl f(x 1 ,x 2 )\ 2 dx 1 dx 2 + R \f(xi, x 2 )\ 2 dx l dx 2 . 

J J R^ JjRd xR d 

Since w € L p , 11 id^ 11 -A 0 when R —> oo and we thus have 

| (/, wf) | < h(R) | (/, - A X1 f) | + (R + h(R)) (/, /), 

where h(R) can be made arbitrarily small by taking R large. In particular we can take 
h(R) < 1, and this leads to an operator inequality of the desired form (13.31) . 

□ 


The above assumptions ensure by well-known methods that Hn is self-adjoint on the 
domain of the Laplacian, and bounded below, see [160]. Our goal is to describe the ground 
state of (EH), i.e. a state achieving) 

E(N) = inf = inf (T, N . (3-6) 

11 ^ 11=1 J 

In the mean-field regime that concerns us here, we expect that any ground state satifies 

Ttv ~ u &N when N -a oo (3.7) 

l^The ground state might not exist, in which case we think of a sequence of states asymptotically 
achieving the inhmum. 
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in a sense to be made precise later, which naturally leads us to the Hartree functional 


1 


Sh M = N 1 (u® N , = (u,Tu)sj + ~{u<8> u,w u <g) u)^ 2 


1 2 

| 2 „„/ „.M„.O.U 2 , 


= f \S7u\ 2 + V\u\ 2 + \- If \u(x)\ 2 vj(x — y)\u(y)\ 2 dxdy. 
ii! 2 J JfixQ 


(3.8) 


We shall denote en and itn the minimum and a minimizer for £h respectively. By the 
variational principle we of course have the upper bound 

E(N) 


N 


< e H 


and we aim at proving a matching lower bound to obtain 

E(N) 


N 


en when IV —>• 00 . 


(3.9) 

(3.10) 


Remark 3.2 (Generalizations). 

All the main ideas can be introduced in the preceding framework, we refer to m for 
a discussion of generalizations. One can even think of the case where both V and w are 
smooth compactly supported functions if one wishes to understand the method in the 
simplest possible case. 

A very interesting generalization consists in substituting the Laplacian in (13.21) by a 
relativistic kinetic energy operator and/or including a magnetic field as described in Sec¬ 
tion Ol One must then adapt the assumptions (13.31) and ([3.4D but the message stays the 
same: the approach applies as long as the many-body Hamiltonian and the limit functional 
are both well-defined. 

Another possible generalization is the inclusion of interactions involving more than two 
particles at a time, to obtain functionals with higher-order non-linearities in the limit. It 
is of course necessary to consider a Hamiltonian in mean-field scaling by adding terms of 
the form e.g. 

A at ^2 w ( x i ~ Xj,Xi - x k ) 

with Aat oc N ~ 2 when N —» 00. It is also possible to take into account a more general 
form than just the multiplication by a potential, under assumptions of the same type 
as (|3.3p . □ 


3.2. Confined systems and the strong de Finetti theorem. 

By “confined system” we mean that we are dealing with a compact setting. We will 
make one the two following assumptions: either 

II C is a bounded set (3.11) 


or 

H = M d and V(x) —> 00 when |x| —> 00 (3-12) 

with V the potential appearing in (|3.2I) . We will also assumte that 

V G Af oc (H), max(l, d/2) < p < 00 . 

In both cases it is well-known m that 


T = — A + V has compact resolvent 


(3.13) 
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which allows to easily obtain strong convergence of the reduced density matrices of a 
ground state of m- One may then really think of the limit object as a quantum state 
with infinitely many particles. Taking inspiration from the classical setting discussed 
previously, the natural definition is the following: 

Definition 3.3 (Bosonic state with infinitely many particles). 

Let L) be a complex separable Hilbert space and, for all n € N, let $ j” be the corresponding 
bosonic n-particles space. We call a bosonic state with infinitely many particles a sequence 
(7^)neN of trace-class operators satisfying 

• y( n ) is a bosonic n-particles state: y^ € ©^fj”) is self-adjoint, positive and 

Tr^[ 7 W] = l. (3.14) 

• the sequence (y^) ra eN is consistent: 

lW[ 7 (n+1) ]= 7 (ri) (3-15) 

where Tr n _|_i is the partial trace with respect to the last variable in T) n+l . 

□ 


The key property is the consistency (13.151) . which ensures that the sequence under 
consideration does describe a physical state. Note that y® is just a real number and that 
consistency implies that Tr^nfyl™)] = 1 for all n as soon as y(°) = 1. 

A particular case of symmetric state is a product state: 


Definition 3.4 (Product state with infinitely many particles). 

A product state with infinitely many particles is a sequence of trace-class operators y( n ) € 
©'(££) with 

y (n) = y® n , (3.16) 

for all n > 0 where y is a one-particle state. A bosonic product state is necessarily of the 
form 

yW = (|u)(u|)® n = \u® n )(u® n \ (3.17) 

with u € ST), the unit sphere of the one-body Hilbert space. □ 


That bosonic product states are all of the form (13.171) comes from the observation 
that if y € © 1 (fj) is not pure (i.e. is not a projector), then y® 2 cannot have bosonic 
symmetry 1041 . 

The strong de Finetti theorem is the appropriate tool to describe these objects and 
specify the link between the two previous definitions. In the following form, it is due to 
Hudson and Moody [94] : 


Theorem 3.5 (Strong quantum de Finetti). 

Let Sj be a separable Hilbert space and (y^)neN a bosonic state with infinitely many 
particles on S). There exists a unique Borel probability measure /i € V(SSf) on the sphere 
SS) = {u € fj, ||it|| = 1} of Sj, invariant under the actioi^f \ of S 1 , such that 


for all n > 0. 


y(n) = f | u ®n)( u ®n| dM ( u ) 

Jss j 


(3.18) 


Multiplication by a constant phase e x6 ,9 £ R. 
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In other words, every bosonic state with infinitely many particles is a convex 
combination of bosonic product states. To deduce the validity of the mean-field 
approximation at the level of the ground state of 0 , it then suffices to show that the 
limit problem is set in the space of states with infinitely many particles, which is relatively 
easy in a compact setting. We are going to prove the following result: 


Theorem 3.6 (Derivation of Hartree’s theory for confined bosons). 

Under the preceding assumptions, in particular (13.111) or (13.121) 

r E(N) 

lim -= ew- 

TV—>• oo N 

Let be a ground state for Hn, achieving the infimum EH) and 

7? :=TV n+1 _,iv[|^)(^|] 


be its n-th reduced density matrix. There exists a unique probability measure p on Mu, 
the set of minimizers of £u (modulo a phase), such that, along a subsequence and for all 
n € N 

lim 7 ^ } = f dp{u) |u® n )(u® n | (3.19) 

N ->°° Jmh 

strongly in the S 1 (f} n ) norm. In particular, if £u has a unique minimizer (modulo a 
constant phase), then for the whole sequence 


= «•>« 


(3.20) 


The ideas of the proof are essentially contained in [701 115411159] , applied to a somewhat 
different context however. We follow some clarifications given in Section 3]. 


Proof. We have to prove the lower bound corresponding to (13.91) . We start by writing 


E(N) 

N 


= ^ (*n,H n * n )^ n = TV^[T 7 «] + ^Trtflvrytf] 


„(2)i 


N 

= ~ Try 
2 ■ 


( 2 ) 


(Ti + T 2 + w) Y n 


(3.21) 


and we now have to describe the limits of the reduced density matrices 7 ^ and 7 ^ . 

Since the sequences ( 7 ^)jvsn are by definition bounded in (5 1 , by a diagonal extraction 
argument we may assume that for all n € N 


(n) 

IN 


7 (n) e6‘(^ 


1 ( 


weakly-* in © 1 (i^ n ). That is, for every compact operator K n over Sj n we have 


Tin 




Tt, 




We are going to show that the limit is actually strong. For this it is sufficient (see [531 fTM] 
or Addendum H]) to show that 


Try 


7 


(n) 


= Ti'v 


(n) 

7iv 


= 1 , 


(3.22) 
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i.e. no mass is lost in the limit. We start by noting that, using Assumption (13.31) . 


eH > ^ Tr hi 


(Ti + T 2 + w) 

l - - * (2)1 


>(l-/3_)^Tr (T 1 + T 2 ) 7 ^ -C TV[ 7 ^ 2) ] 


„(!) 


-c. 


= (l-/3_)Tr[T 7 ] v 

Thus Try-,[T 7 ^ ] is uniformly bounded and, up to a possible further extraction, we have 

(T + C 0 ) 1 ' 2 7 {J } (T + Cq) 1 ' 2 (T + C 0 ) 1/2 7 (1) (T + O 0 ) 1/2 
for a certain constant Co- Consequently 

1 = Tr 7 [ 7 y } ] = Tr 7 [(T + Co )" 1 (T + Co) 1 ' 2 7 g } (T + O 0 ) 1/2 " 

Tr 7 [(T + Co )" 1 (T + C 0 ) 1/2 7 (1) (T + Co) 1 ' 2 ] = TY*[ 7 «] 

since (T + Co ) -1 is by Assumption (13.131) a compact operator. One obtains ([3.221) similarly 
by noting that 

1 

n 




(n) 


T. T >rS 

3 = 1 


is uniformly bounded in IV and that l Tj also has compact resolvent which allows for 
a similar argument. 

We thus have, for all n G N 

v( n ) 


(n) 

Tn 


7 


strongly in trace-class norm, and in particular, for all bounded operator B n on Sj n 

TVbS’Sn] -> TV^nlyW^]. 

Testing this convergence with B n+ \ = B n ( 8 ) 1 we deduce 

TY n+1 [ 7 ^ +1 )] = 7 W 

and thus the sequence ( 7 ^ nl )neN describes a bosonic state with infinitely many particles 
in the sense of Definition 13.31 We apply Theorem 13.51 which yields a measure /r G V(SSj). 
In view of Assumption (13.31) . the operator T\ + T 2 + w is bounded below on F) 2 , say by 
2Cr- Since Tr f -,2 7 1 2 * = 1 we may write 


lim inf — Tr f 

iV—>-oo 2 


( 2 )' 


(ii + r 2 + w) y n 


= lim inf — TV? 

N— >oo 2 

> - Tr 
“ 2 

= - Tr 
2 


( 2 ) 


(Ti + T 2 + w- 2 Or ) 7 )v 


+ Op 


h 2 


(Ti + T 2 + w; — 20p) 7 ® 
(Ti +T 2 + w) 7 (2) 


+ o T 


using Fatou’s lemma for positive operators. Using the linearity of the energy as a function 
of 7 ( 2 ) and ([3.181) 

1 Tr& [(T x + T 2 + w) \u® 2 )(u® 2 \] df3(u) = 


r E(N) 

lim inf —— > 

N—t 00 N 


'uGSF) 


' u£Sf) 


£n[u]dn(u) > e H , 





























50 


NICOLAS ROUGERIE 


which is the desired lower bound. The other statements of the theorem follow by noting 
that equality must hold in all the previous inequalities and thus that p may charge only 
minimizers of Hartree’s functional. □ 

It is clear from the preceding proof that the structure of bosonic states with infinitely 
many particles plays the key role. The Hamiltonian itself could be chosen in a very abstract 
way provided it includes a confining mechanism allowing to obtain strong limits. Several 
examples are mentioned in |113l Section 3]. 


3.3. Systems with no bound states and the weak de Finetti theorem. 

In the previous section we have used in a strong way the assumption that the system 
was confined in the sense of (K3. 1 111 - (13.12f) . These assumptions are sufficient to understand 
many physical cases, but it is highly desirable to be able to relax them. One is then 
lead to study cases where the convergence of reduced density matrices is not better than 
weak-*, and to describe as exhaustively as possible the possible scenarii, in the spirit of 
the concentration-compactness principle. A first step, before asking the question of how 
compactness may be lost, consists in describing the weak limits themselves. It turns out 
that we still have a very satisfying description. In fact, one could not have hoped for 
better than the following theorem, proven in |1131 : 

Theorem 3.7 (Weak quantum de Finetti). 

Let S) be a complex separable Hilbert space and (T n)n&n cl sequence of bosonic states with 
Tat € We assume that for all n € N 

T ^ 7 W (3.23) 

in Then there exists a unique probability measure p € 'P(Hi^) on the unit ball 

B. f) = {u € Sj, ||u|| < 1} of Sj, invariant under the action of S 1 , such that 


for all n > 0. 


» 



u® n ){u® n \dp{u) 


(3.24) 


Remark 3.8 (On the weak quantum de Finetti theorem). 

(1) Assumption (|3.23l) can always be satisfied in practice. Modulo a diagonal extrac¬ 
tion, one may always assume that convergence holds along a subsequence. This 
theorem thus exactly describes all the possible weak limits for a sequence of bosonic 
A-body states when N —>• oo. 

(2) The fact that the measure lives over the unit ball in (I3.24P is not surprising since 
there may be a loss of mass in cases covered by the theorem. In particular it is 
possible that 'yW) = o for all n and then p = 5 q, the Dirac mass at the origin. 

(3) The term weak refers to “weak convergence” and does not indicate that this result 
is less general than the strong de Finetti theorem. It is in fact more general. To 
see this, think of the case where no mass is lost, Tr 7 n[ 7 ( n )] = 1. The measure p 
must then be supported on the sphere and convergence must hold in trace-class 
norm. It is by the way sufficient to assume that Tr^u [yW] = 1 for a certain n € N, 
and the convergence is strong for all n since the measure p does not depend on n. 
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(4) Ammari and Nier have slightly more general results [HEEIES- In particular, it is 
not necessary to start from a state with a fixed particle number. One can consider 
a state on Fock space provided suitable bounds on its particle number (seen as a 
random variable in this framework) are available. 

(5) Uniqueness of the measure follows from a simple argument. Here we will mostly 
be interested in the existence, which is sufficient for static problems. For time- 
dependent problems on the contrary, uniqueness is crucial OT0I3SI. 

□ 


Coming back to the derivation of Hartree’s theory, let us recall that the energy upper 
bound (13.91) is always true. We are thus only looking for lower bounds. A case where 
knowning the weak-* limit of reduced density matrices is then that of a weakly-* lower 
semi-continuous functional. This remark may seem of marginal interest, but this case does 
cover a number of physically relevant systems, those with no bound states. 

We are going to prove the validity of Hartree’s approximation in this case, using Theo¬ 
rem [373 We shall work in and assume that the potential V in (|3.2I) is non-trapping in 
every direction: 

V € L p (M d ) + L°°(M d ), max 1, d/2 < p < 00 , V(x) —>• 0 when |x| —> 00 . (3.25) 

The assumption about the absence of bound states concerns the interaction potential w. 
It is materialized par the inequality^ 

w 

— A + — > 0 (3.26) 

as operators. This means that w is not attractive enough for particles to form bound 
states such as molecules etc ... Indeed, because of Assumption (13.41) . — A + ^ can have 
only negative eigenvalues, its essential spectrum starting at 0. In view of (13.261) . it can in 
fact not have any eigenvalue at all, and thus no eigenfunctions which are by definition the 
bound states of the potential. A particular example is that of a purely repulsive potential 
w > 0. 


Under Assumption (13.261) . particles that might escape to infinity no longer see the one- 
body potential V and then necessarily carry a positive energy that can be neglected for a 
lower bound to the total energy. Particles staying confined by the potential V are described 
by the weak-* limits of density matrices. We can apply the weak de Finetti theorem to 
the latter, and this leads to the next theorem, for whose statement we need the notation 


eH(A):= 



0 < A < 1 


(3.27) 


for the Hartree energy in the case of a loss of mass. Under assumption (13.261) it is not 
difficult to show that for all 0 < A < 1 


ch(A) > e H (l) = e H 

by constructing trial states made of two well-separated pieces of mass. 


(3.28) 


an appropriate variant when one considers a different kinetic energy, cf Remark 13.21 
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Theorem 3.9 (Derivation of Hartree’s theory in the absence of bound states). 

Under the previous assumptions, in particular (13.251) and (13.261) we have 


lim 

N—> oo 


EjN) 

N 


= e H - 


Let ^n be a sequence of approximate ground states, satisfying 

(^ n ,H n ^ n ) = E(N) + o(N) 


in the limit N —>■ oo and 

7 & } :=Tr n+1 ^ N [\y N ){* N \} 

be the corresponding n-th reduced density matrices. There exists a probability measure p 
supported on 

Mu = ju G BS), 8u[u\ = eu (j|tt|| 2 ) = en(l)| 
such that, along a subsequence and for all n € N 


(n) 

Tn 


Wh 


dp[u) |u® n )(u® n | 


in ©^fj"). In particular, if eu has a unique minimizer uu, with ||uh| 
whole sequence, 


N lim 7"> = |»g”> («r 


strongly in trace-class norm. 


(3.29) 
1, then for the 

(3.30) 


Remark 3.10 (Case where mass is lost). 

Note that it is possible for the convergence in the above statement to be only weak-*, 
which covers a certain physical reality. If the one-body potential is not attractive enough 
to retain all particles, we will typically have a scenario where 

f ch(A) = e H (l) for A c < A < 1 
\e H (A) < e H (l) for 0 < A < A c 

where A c is a critical mass that can be bound by the potential V. In this case, en(A) will 
not be achieved if A c < A < 1 and one will have a minimizer for Hartree’s energy only 
for a mass 0 < A < A c . If for example the minimizer uu at mass A c is unique modulo a 
constant phase, Theorem 13.91 shows that 

iP 1-rx-n 

and one should note that the limit has a mass A” < 1. This scenario actually happens in 
the case of a “bosonic atom”, see Section 4.2 in |TT3] . □ 

Theorem 13.91 was proved in [ 113] . In order to be able to apply Theorem 13.71 we start 
with the following observation: 


Lemma 3.11 (Lower semi-continuity of an energy with no bound state). 

Under the previous assumptions, let 7^,7^ > 0 be two sequences satisfying 

Ti >>2 = 1, tJ? = Tr 2 7Jv } 
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as well as 7 ^ —s 7 ^) weakly-* in 6 1 (Sj k ) for k = 1,2. Then 

liminf (t^[T^) + > T^[T 7 «] + ±Trtf[vry®]. (3.31) 

TV—>-oo \ Z J Z 

Proof. We shall need two trunctation functions 0 < XRi Vr < 1 satisfying 

Xr + Vr = l.supp (xr) C i?(0, 2R), supp(rm) c B(0, R) c , |Vxh| + |V?? R | < CR~ l . 

It is easy to show the IMS formula: 

- A = xr(~A)xr + vr(-A)t] R ~ |Vxr| 2 - |V77^| 2 

> Xr(-A)xr + rj R (-A)riR - CR~ 2 (3.32) 

as an operator. For the one-body part of the energy we then easily have 

> Tr ^XR^nXr] +Tr j --,[-A?? i? 7^ ) r/ jR ] +n(N,R) 


with 

n(N, R) > -CR - 2 + TrninkV'yff] 


and thus 


lim inf lim inf 77 ( N , R) = 0 

R—y oo N^-oo 

because in view of (13.2511 we have 

Tr ^[VrYIn^] = / Vr( x )V( x )p^n ( x )^ x 0 when R —>• 00 (3.33) 

JsL d 


uniformly in N. Here we have denoted pff the one-body density of 7 ^* , formally defined 
by 


(1)/ \ (1)/ X 

P N ( x ) = Tn (®>®)> 


where we identify 7 ^ and its kernel. The details of the proof of (|3.33l) are left to the 
reader. One should use the a priori bound TrfT^j^] < C to obtain a bound on pff in 
some appropriate L q space via Sobolev embeddings, then use Holder’s inequality and the 
fact that, for the conjugate exponent p, 


hR V \\LP(R d ) 0 


by assumption. 

To deal with the interaction term we introduce 


Pn\ x ,V ) := Tn\ x i Vi x i v) 
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( c y\ (n\ 

the two-body density of 7 ]^ (identified with its kernel r f N {x ,y \ x,y)) and write 
T rrA w 'y < N ) ]= II w{x - y)p$(x,y)dxdy 

' <K d 

w{x - y)x 2 R(x)PN\x,y)x 2 R(y)dxdy 

<K d 

+ // w(x - y)r] R (x)p^ ] (x,y)rj 2 R (y)dxdy 

:K d 

+ 2 // w(x - y)xn(x)/„ ] (x,y)y\(ji)ixiy 

+ r 2 (N,R ) 

(2) (2) 

= Tiy.,2 [w XR ® XRTn XR ® Xr\ + Tr [w r) R (S> VrT n VR ® %] 


+ r 2 (iV,i?) 


where 


lim inf lim inf r 2 (N, R ) = 0 

R—»oo AT—xoo 


by a standard concentration-compactness argument that we now sketch. One has to show 
that 

„( 2 ), 


lim lim 

R—^oo N —>oo 


- y)\x R (x)p N ( x ,y)riR(y)dxdy. 


Since x R (x)\w(x - y)\yj R (y) < ^s-j^rM 21 ~ 2/)|, the term 

\w{x - 2 /)|x|(x)p^ ) (x,y)r 7 | R (y)dxdy 

can be dealt with similarly as (13.331) . We leave again the details to the reader. There thus 
remains to control 

[[ \w(x - y)\x 2 R (x)p%\x,y) (: r? R {y ) - vIr(v)) dxdy. 

J JR d xl d 

for which we claim that 

K® - y)\XR(x)[v R (y) ~ yl R {y)\p { x > k (x,y)dxdy 

\w(x - y)\t(R < \y\ < 8 R)p ( $ k (x,y) dxdy) = 0 
for an appropriate subsequence (iV*.). Let us introduce the concentration functions 
Qn{R) ■= [[\w(x-y)\t(\y\>R)p < $(x,y)dxdy. 


lim lim 

R—> oo k—yoo 


< lim lim 

R —>-00 fc—>-oo 
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For every N, the function R i-a Qn(R) is decreasing on [0, oo). Moreover, 

o < Qn(R) < Tr fl2 [|«;| 7 ^ ) ] < C 0 

using our assumptions on w and the a priori bound Tr T 7 j^ < C again. Therefore, 
by Helly’s selection principle, there exists a subsequence and a decreasing function 
Q : [ 0 , oo) —> [0, Co] such that QN k {R ) —> Q{R) for all R € [ 0 , oo). Since limR^. 00 Q(R) 
exists, we conclude that 

lim lim (Q Nk (R) — Q]y k (8R )) = lim (Q(R) - Q{8R)) = 0 . 

R —»oo k— too R —»oo 

This is the desired convergence. 

At this stage we thus have 


HminfTr j j[T 7 jJ ) ] + \ Tr^^JJ 0 ] > 

iv—>-oo Z 

iim inf lim inf Tr$ [Txr^nXr] + \ Tr f:i 2 [w XR ® XRJ^nXR ® Xr] 

R—^oo iV—>-oo Z 

+ lim inf lim inf Tr fl [- Ar/Rj^rj R ] + - Tr fl2 [w r) R ® r] R 'y ( Z ) VR ® Vr] (3-34) 

R—s- oo N^oo Z 

The terms on the second line give the right-hand side of (13.311) . Indeed, recalling that 
T = —A + V, using Fatou’s lemma for —A > 0 and the fact that 

XrJnXR -> XR7 (1) XR 

in norm since xjt has compact support, we have 

lim inf TYq [T XrI^n Xr\ > ^[TxrX^Xr]- 

N—>oo 

It then suffices to recall that xr 1 pointwise to conclude. The interaction term is dealt 
with in a similar way, using the strong convergence 

xr ® xrIn^xr ®xr^ xr® xri^xr ® xr 

and then the pointwise convergence xr 1 - 

The terms on the third line of (I3.34p form a positive contribution that we may drop 
from the lower bound. To see this, note that since 0 < r] R < 1 we have r/ R <g)l > t) R (% rq R 
and (with Ti '2 the partial trace with respect to the second variable) 

vrIn^vr ^ Tr 2 [m ® vrtnVr ® vr], 

( 2 ) 

which gives, by symmetry of 7 b' , 


TXsj [-^VRX^nVr] + 7} TV [w VR ® VRI^nVR ® Vr] 


> - Tr f 


((-A) < 8 > 1 + 1 < 8 ) (—A) +w) r] R <S) VRlffvR ® VR 


and it is not difficult to see that Assumption (13.261) implies 

(—A) < 8 ) 1 + 1 ( 8 ) (—A) + w > 0 


l^Just decouple the center of mass from the relative motion of the two particles, i.e. make the change 
of variables ( X\,X 2 ) eA (x\ + X 2 ,x\ —* 2 ). 
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which ensures the positivity of the third line of (13.341) and concludes the proof. □ 

We now conclude the 

Proof of Theorem \S. 91 Starting from a sequence V at = |4 / at}('I , aH °f ^V-body states we 
extract subsequences as in the proof of Theorem 13.61 to obtain 

L N *7 

Using Lemma 13.Ill we obtain 

lim inf IV -1 Tr^jv [H^T n ] = linr inf ( Try-, [TT^] + ^ Tr f -,2 [wT ( f ’]) 

W—>-oo IV—kx> y 2 J 

>TV fl [T 7 «] + lTV fl2 [u; 7 ( 2 )] 

and there remains to apply Theorem 13.71 to the sequence ( 7 ^ n ^) n eN to obtain 

liminf 77 _1 Tr f jv^atTat] > / £u[u\dfi(u) > en 
Af^oo J BS) 

using (13.281) and the fact that f BS , dfr(u) = 1. Once again, the other conclusions of the 
theorem easily follow by inspecting the cases of equality in the previous estimates. □ 

For later use, we note that during the proof of Lemma 13.111 we have obtained the 
intermediary result (I3.34|) without using the assumption that w has no bound state. We 
write this as a lemma that we shall use again in Chapter [6l 

Lemma 3.12 (Energy localization). 

Under assumptions m and (13.251) , let Tat be a sequence of almost minimizers for E(N): 

(^ n ,H n ^ n ) = E(N)+o(N) 

(k) 

and 7 yy the associated reduced density matrices. We have 

bminf = lim inf Tr^Tyj^] + \ Tr fl2 > 

iV—>• OO iV N —>-oo Z 

lim inf lim inf TV 7 [Txrt'nXr] + \ Tr fl2 [w Xr 2 '7n ) Xr 2 ] 

H —>-oo iv— >-oo Z 

+ lim inf lim inf Tr fl [-A tirTn^r] + \ Tr 7 2 [w V^T^Vr 2 ] (3-35) 

R—>oc iv—>oo Z 

where 0 < \R < 1 is C 1 , with support in B( 0, R), and r/n = \J 1 — x\- 

3.4. Links between various structure theorems for bosonic states. 

We have just introduced two structure theorems for many-particles bosonic systems. 
These indicate that, morally, if Tat is a IV-body bosonic state on a separable Hilbert 
space f), there exists a probability measure /u € V{$)) on the one-body Hilbert space such 
that 

f \u® n ){u® n \dpi(u) (3.36) 

J udS) 

when N is large and n is fixed. Chapters [5] and [5] of these notes are for a large devoted to 
the proofs of these theorems “a la de Finetti”. As mentioned in Remark 13.81 Theorem 13.71 
is actually more general than Theorem 13.51 and we will thus prove the former. 
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In order to better appreciate the importance of the weak theorem in infinite dimensional 
spaces, we emphasize that the key property allowing to prove the strong theorem is the 
consistency (13.151) . When starting from the reduced density matrices of a JV-body 
state T n and extracting weakly-* convergent subsequences to define a hierarchy (7 n ) ngN > 


~p( n ) 

1 N 



one only has 


Tr 


71+1 


(nH-3.) 


< lim Tr n+ i 
iV —>-00 


tST 1 


= lim iP = 7 <“>. 


TV—>-oo 


because the trace is not continuous^ for the weak-* topology , only lower semi-continuous. 
Obviously, the relation 


Tr 


71+1 


7 


O+i) 


< 7 


(n) 


(3.37) 


is not sufficient to prove a de Finetti theorem. A simple counter-example is given by the 
sequence of reduced density matrices of a one-body state v € SS)\ 

7 ^ 0) = 1 , 7 *- 1 ’ = |u)(u|, 7 ^ = 0 for n > 2 . 


In these notes we have chosen to prove the weak de Finetti theorem as constructively 
as possible. Before we announce the plan of the proof, we shall say a few words of the 
much more abstract approach of the historical references [192 . I9i] , These papers actually 
contain a version of the theorem that is even stronger than Theorem 13.71 This result 
applies to “abstract” states that are not necessarily normal or bosonic. These may be 
defined as follows: 

Definition 3.13 (Abstract state with infinitely many particles). 

Let h be a complex separable Hilbert space and for all n £ N, fj n = (£) n Sj be the 
corresponding n-body space. We call an abstract state with infinitely many particles a 
sequence (+ T + ne N where 

• u/ n ) is an abstract n-body state: u +) € (23(f) n ))*, the dual of the space of bounded 
operators on $j n , > 0 and 

+ n )(l7») = 1. (3.38) 

• £>) is symmetric in the sense that 

J n) (Hi (gi... ® B n ) = (B a{x) ® ... ® B a(n] ) (3.39) 

for all Hi,..., B n £ 23(f)) and all permutation o £ E n . 

• the sequence (u+)) ne N is consistent: 

J n+ V (Hi ® ... ® H„ ® 1+ = + n) (Hi (gi... (8) B n ) (3.40) 

for all Hi,..., B n £ 23(f)). 

□ 

An abstract state is in general not normal, i.e. it does not fit in the following definition: 
l^This is a fancy way of characterizing infinite dimensional spaces. 
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Definition 3.14 (Normal state, locally normal state). 

Let fj be a complex separable Hilbert space and (a/ n ^) n eN be an abstract state with 
infinitely many particles. We say that (u/ n )) raGN is locally normal if is normal for all 
n € N, i.e. there exists €E © 1 (i5 n ) a trace-class operator such that 

J n \B n ) = Tr^[^B n ] (3.41) 

for all B n € ©(£ n ). □ 

Identifying trace-class operators with the associated normal states we readily see that 
Definition 13.31 is a particular case of abstract state with infinitely many particles. Note 
that (by the spectral theorem), the set of convex combinations of pure states (orthogonal 
projectors) coincides with the trace class. A non-normal abstract state is thus not a mixed 
state, i.e. not a statistical superposition of pure states. The physical interpretation of a 
non-normal state is thus not obvious, at least in the type of settings that these notes are 
concerned with. 

The consistency notion (13.401) is the natural generalization of (13.151) but it is important 
to note that the symmetry (13.391) is weaker than bosonic symmetry. It in fact corresponds 
to classicaly indistinguishable particles (the modulus of the wave-function is symmetric 
but not the wave-function itself). One may for example note that if u >^ is normal in the 
sense of (13.411) . then 7 ^ € 6 1 (L) n ) satisfies 

u^u; = 7 (n) 

where U a is the unitary operator permuting the n particles according to a € H n . Bosonic 
symmetry corresponds to the stronger constraint 

u a 7 (n) = 7 H u; = 7 (ti) , 

cf Section 11.21 

We also have a notion of product state that generalizes Definition 13.41 

Definition 3.15 (Abstract product state with infinitely many particles). 

We call abstract product state an abstract state with infinitely many particles such that 

uW = u® n (3.42) 

for all n € N, where oj € (03(f)))* is an abstract one-body state (in particular u > 0 and 

= !)• Q 

The most general form of the quantum de Finetti theorenJ*^ says that every abstract 
state with infinitely many particles is a convex combination of product states: 

Theorem 3.16 (Abstract quantum de Finetti). 

Let S) be a complex separable Hilbert space and (u/ n )) nG N on abstract bosonic state 
with infinitely many particles built on f). There exists a unique probability measure 
[i G V ((5S(ij))*) on the dual of the space of bounded perators on Sj such that 

pi ({w € (»(«))*, w > 0, w(l fl ) = 1}) = 1 (3.43) 

and 

J uj® n dfi(io) (3.44) 


W ("> = 


least the most general form the author is aware of. 
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for all n > 0. 

Remark 3.17 (On the abstract quantum de Finetti theorem). 

(1) This result was first proven by Stprmer [ 192] , Hudson and Moody [91] then gave a 
simpler proof by adapting the Hewitt-Savage proof of the classical de Finetti theo- 
rem l2.ll they prove that product states are the extremal points of the convex set of 
abstract states with infinitely many particles. The existence of the measure is then 
a consequence of the Choquet-Krein-Milman theorem. This approach does require 
the notion of abstract states and does not give a direct proof of Theorem 13.51 

(2) Hudson and Moody [M] deduce the strong de Finetti theorem from the abstract 
theorem. An adaptation of their method (see Appendix A]) shows that the 
weak de Finetti theorem is also a consequence of the abstract theorem. 

(3) This theorem has been used to derive Hartree-type theories for abstract states 
without bosonic symmetry in [701115911154] . To recover the usual Hartree theory 
one must be able to show that the limit state is (locally) normal. In finite dimen¬ 
sional spaces, 03 (fj) of course coincides with the space of compact operators, which 
implies that any abstract state is normal. This difficulty thus does not occur in 
this setting. 

(4) There also exist quantum generalizations of the probabilistic versions of the clas¬ 
sical de Finetti theorem (that is, those dealing with sequences of random vari¬ 
ables 13 isnj) , see e.g. HUT]. These are formulated in the context of free probability 
and require additional symmetry assumptions besides permutation symmetry. 

□ 


At this stage we thus have the scheme (“deF” stands for de Finetti) 


abstract deF =>• weak deF strong deF 


(3.45) 


but the proof of the weak de Finetti theorem we are going to present follows a different 
route, used in [11311116] , It starts from the finite dimensional theorem^: 


finite-dimensional deF weak deF strong deF 


(3.46) 


This approach leads to a somewhat longer proof than the scheme (I3.45p starting from the 
Hudson-Moody proof of Theorem 13.161 This detour is motivated by five main practical 
and aesthetic reasons: 


(1) The proof following (13.461) is simpler from a conceptual point of view: it requires 
neither the notion of abstract states nor the use of the Choquet-Krein-Milman 
theorem. 

(2) Thanks to recent progress, due mainly to the quantum information community [391 
iMimnm96ii7Hfii6]. we have a completely constructive proof the finite dimensional 
quantum de Finetti theorem at our disposal. One first proves explicit estimates 
thanks for a construction for finite N, in the spirit of the Diaconis-Freedman 
approach to the classical case. Then one passes to the limit as in the proof of the 
Hewitt-Savage theorem we have presented in Section 12.21 


16 


In finite dimension there is no need to distinguish between the weak and the strong version. 
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(3) The first implication in Scheme (13.461) is also essentially constructive, thanks to 
Fock-space localization techniques used e.g. in [3l 15511112] . These tools are inher¬ 
ited from the so-called “geometric” methods [66, 63 118011182] that adapt to the 
iV-body problem localization ideas natural in the one-body setting. These allow 
(amongst other things) a fine description of the lack of compactness due to loss of 
mass at infinity, in the spirit of the concentration-compactness principle [ 136 ] . 137 ]. 

(4) In particular, the proof of the first implication in (I3.4t>|) yields a few corollar¬ 
ies which will allow us to prove the validity of Hartree’s theory in the general 
case. When the assumptions made in Section T3. 31 do not hold, the weak de Finetti 
theorem and its proof according to (13.45]) are not sufficient to conclude: Parti¬ 
cles escaping to infinity may form negative-energy bound states. The localization 
methods we are going to discuss will allow us to analyze this phenomenon. 

(5) In Chapter [7] we will deal with a case where the interaction potential depends on 
N to derive non-linear Schrodinger theories in the limit. This amounts to taking 
a limit where w converges to a Dirac mass simultaneously to the IV —» oo limit. In 
this case, compactness arguments will not be sufficient and the explicit estimates 
we shall obtain along the proof of the finite dimensional de Finetti theorem will 
come in handy. 

An alternative point of view on the proof strategy ([3.461) is given by the Ammari-Nier 
approach [H El El E]> based on semi-classical analysis methods. The relation between the 
two approaches will be discussed below. 

As noted above, the second implication in (13.461) is relatively easy. The following chap¬ 
ters deal with the first two steps of the strategy. They contain the proof of the weak de 
Finetti theorem and several corollaries and intermediary results. The finite dimensional 
setting (where the distinction between the strong and the weak theorems is irrelevant) is 
discussed in Chapter [U The localization methods allowing to prove the first implication 
in (13.461) are the subject of Chapter [5] 
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4. The quantum de Finetti theorem in finite dimensonal spaces 

This chapter deals with the starting point of the proof strategy (13.4611 , that is a proof of 
the strong de Finetti theorem in the case of a finite dimensional complex Hilbert space S), 

dimf) = d. 

In this case, strong and weak-* convergences in S 1 ^ 71 ) are the same and thus there is 
no need to distinguish between the strong and the weak de Finetti theorem. The main 
advantage of working in finite dimensions is the possibility to obtain explicit estimates, in 
the spirit of the Diaconis-Freedman theorem (with a completely different method, though). 
We are going to prove the following result, which gives bounds in trace-class norm: 


Theorem 4.1 (Quantitative quantum de Finetti). 

Let Tv € S 1 ^^) be a bosonic state over and yjy its reduced density matrices. There 
exists a probability measure € V(SS)) such that, denoting 

f N = [ \u® N )(u® N \dp N (u) (4.1) 

J u£S$) 

(n) 

the associated state and ' its reduced density matrices, we have 


for all n = 1... N. 


Tr^n 




2 n(d + 2 n) 
N 


(4.2) 


Remark 4.2 (On the finite dimensional quantum de Finetti theorem). 

(1) This result is due to Christandl, Konig, Mitchison and Renner [39], important 
earlier work being found in [106] and [71j . One may find developments along 
theses lines in [381I901HT6]. The quantum information community also considered 
several variants, see for example [33], [39] [40] [44] 11621 [27] , 

(2) One can add a step in the strategy (13.4611 : 


quantitative deF => finite dimensional deF => weak deF => strong deF 


(4.3) 


Indeed, in finite dimension one may identify the sphere SS) with a usual, compact, 
Euclidean sphere (with dimension 2d — 1, i.e. the unit sphere in M? d ). The space of 
probability measures on SS) is then compact for the usual weak topology and one 
may extract from /cv a converging subsequence to prove Theorem 13.51 in the case 
where dim S) < oo, exactly as we did to deduce Theorem 12.11 from Theorem 12.21 in 
Section 12.21 

(3) The bound (14.211 is not optimal. One may in fact obtain the estimate 


Trsyi 




< 


2nd 

1 \T’ 


(4.4) 


with the same construction, see [Ml EH] EES]. The proof we shall present only 
gives (14.211 but seems more instructive to me. For the applications we have in 
mind, n will always be fixed anyway (equal to 2 most of the time), and in this 
case (14. 21) and (14.41) give the same order of magnitude in terms of N and d. 
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(4) The bound (14.41) in the quantum case is the equivalent of the estimate in dn/N 
mentioned in Remark IQ for the classical case. One may ask if this order of 
magnitude is optimal. It clearly is with the construction we are going to use, but 
it would be very interesting to know if one can do better with another construction. 
In particular, can one find a bound independent from d, reminiscent of ([2.141) in 
the Diaconis-Freedman theorem ? 

□ 

The construction of Tj\f is taken from [39]. It is particularly simple but it does use 
in a strong manner the fact that the underlying Hilbert space has finite dimension. The 
approach we shall follow for the proof of Theorem 14.11 is originally due to Chiribella [38]. 
We are going to prove an explicit formula giving the density of Tat as a function of those of 
T]\r, in the spirit of Remark 12.31 This formula implies ([4.2D in the same manner as ([2.19D 
implies (12.141) . 

In Section |TT] we present the construction, state Chiribella’s explicit formula and obtain 
Theorem im as a corollary. Before giving a proof of Chiribella’s result, it is useful to 
discuss some informal motivation and some heuristics on the Christandl-Konig-Mitchison- 
Renner (CKMR) construction, which happens to be connected to well-known ideas of 
semi-classical analysis. This is the purpose of Section 14.21 Finally, we prove Chiribella’s 
formula in Section [4.31 following the approach of |116j . It has been independently found 
by Lieb and Solovej [ 129] (with a different motivation), and was inspired by the works 
of Ammari and Nier El El E]. Related considerations appeared also in [104 . Chapter 3]. 
Other proofs are available in the literature, cf [3!] and [90]. 

4.1. The CKMR construction and Chiribella’s formula. 

We first note that the Diaconis-Freedman construction introduced before is purely classi¬ 
cal since it is based on the notion of empirical measure, which has no quantum counterpart. 
A different approach is thus clearly necessary for the proof of Theorem 14.11 

In a finite dimensional space, one may identify the unit sphere SSj = {u € Sj, ||rt|| = 1} 
with a Euclidean sphere. One may thus equip it with a uniform measure (Haar measure of 
the rotation group, simply the Lebesgue measure on the Euclidean sphere), that we shall 
denote du , taking the convention that 

I du = 1. 

Jssj 

We then have a nice resolution of the identity as a simple consequence of the invariance 
of du under rotations (see Section IT31 below for a proof). We state this as a lemma: 

Lemma 4.3 (Schur’s formula). 

Let Sj be a complex finite dimensional Hilbert space and S) N the corresponding bosonic 
N-body space. Then 

dimf}^ f \u® N )(u® N \du = tfiN. (4.5) 

Jssj 

The idea of Christandl-Konig-Mitchison-Renner is to simply define 
dpisr(u) := dimfj^Tr^iv [Fjv|u® Ar )(u l8,Ar |] du 
= dimf)]^ Tr^jv (u ® N , Tnu (Z> n ) du, 


(4.6) 
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i.e. to take 


r N = dim I \u® N )(u® N \ (u® N , T n u ® n > du. 
Jssj 

Chriibella’s observatior0 is the following: 

Theorem 4.4 (Chiribella’s formula). 

With the previous definitions, it holds 


(4.7) 


~(n) _ (N + n + d - 1 \ ST ( N \ M ~ ^ 

T'jv - ( ) ( / ] 7jV 

V ' f=0 v 7 


(4.8) 


with the convention 

Tn ®s 1 S) n ~ e — (7jv)<t( 1),...,<7(£) ® (l^n—«)o-(^+l),...,(T(n) 

(TtzSn 

where (7Ar)o-(i),...,o-(Q «c£s 072 <r(l) • • • >cr(^) variables. 

From this result we deduce a simple proof of the quantitative de Finetti theorem: 

Proof of Theorem \ f.l\ We proceed as in (12.2011 . Only the first term in the sum (|4.8f) is 
really relevant: 

7 N ] - lN ] = (C(d, n, N ) - 1)7^ + B = —A + B 

where 

and A, B are positive operators. We have 


(4.9) 


Tr^n [—A + B] = Tr 
and thus, by the triangle inequality, 


~ (n) (n) 

Tn ~Tn 


= 0, 


Tr 


~ (n) (n) 

Tn ~Tn 


<TrA + TrB = 2TrA = 2(l- C(d,n,N )). 


Next, the elementary inequality 


n— 1 


C(d,n,N) = \ 


3=0 


N-j 
N+j + d 


> 1 - 


2n + d — 2 
N + d + n — 1 


> 1 — n 


2n + d — 2 
N + d + n — 1 


gives 


which is the desired result. 


Tr 


(n) 

Tn 


~ (n) 
IN 


< 


2n(d + 2 n) 


N 


(4.10) 

□ 


Everything now relies on the proof of Theorem 14.41 which is the subject of Section ll.ol 
Before we give it, some heuristics regarding the relevance of the construction (14.71) shall 
be discussed. 

'On the quantum information vocabulary this is formulated as a relation between “optimal cloning” 
and “optimal measure and prepare channels”. 
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4.2. Heuristics and motivation. 

Schur’s formula (14.511 expresses the fact that the family (u® N ^ forms an over¬ 
complete basis of . Such a basis labeled by a continuous parameter is reminiscent of a 
coherent state decomposition (1031.1197] , This basis in fact turns our to be “less and less 
over-complete” when N gets large. Indeed, we clearly have 

(u® N ,v® N )fiN = 0 when IV —> oo (4-11) 

as soon as u and v are not exactly colinear. The basis (u® N ) ueS ^ thus becomes “almost 
orthonormal” when N tends to infinity. 

In the vocabulary of semi-classical analysis [1201 [1841 [T6] , trying to write 



amounts to looking for an upper symbol representing Tjv- In fact, it so happens [184 ] 
that one may always find such a symbol, only ///v is in general not a positive measure. The 
problem we face is to find a way to approximate the upper symbol (for which no explicit 
expression as a function of the state itself exists, by the way) with a positive measure. 

Note that if the coherent state basis were orthogonal, then the upper symbol would 
be positive, simply by positivity of the state. Since we noted that the family ( u ® N ) ueSs - :i 
becomes “almost orthonormal” when N gets large, it is very natural to expect that the 
upper symbol may be approximate by a positive measure in this limit. 

On the other hand, the measure introduced in (14.611 is exactly what one calls the lower 
symbol of the state Tv- One of the reasons why lower and upper symbols were introduced 
is that these two a priori different objects have a tendency to coincide in senri-classical 
limits. But the N —>• oo limit we are concerned with may indeed be seen as a senri-classical 
limit and it is thus very natural to take the lower symbol as an approximation of the upper 
symbol in this limit. 

One can motivate this choice in a slightly more precise way. Assume we have a sequence 
of IV-body states defined starting from an upper symbol independent of IV, 

T JV = dim(iof) [ fi sup {u)\u® N ){u® N \du, 

J uOSS) 

and let us compute the corresponding lower symbols: 

Vn(v) = (v m ,r N v® N ) = dim(iof) [ // u » \(u® N ,v® N )\ 2 du 

J uOSSj 

= dim (ft?) [ p, sup {u)\(u,v)\ 2N du. 

J uOSSj 

In view of the observation (14.IIP and the necessary invariance of /r sup under the action of 
S 1 it is clear that we have 

/^ f (u) -A // up (u) when N ^ oo. 

In other words, the lower symbol is, for large IV, an approximation of the upper sym¬ 
bol, that has the advantage of being positive. Without consituting a rigorous proof of 
Theorem I4TI this point of view shows that the CKMR constuction is extremely natural. 
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4.3. Chiribella’s formula and anti-Wick quantization. 

The proof of Theorem 14.41 we are going to present uses the second quantization formal¬ 
ism. We start with a very useful lemma. In the vocabulary alluded to in the previous 
section it says that a state is entirely characterized by its lower symbol, a well-known 
fact pm. 

Lemma 4.5 (The lower symbol determines the state). 

If an operator 7 ^ on Sj k satisfies 

{u® k ,^u® k ) = 0 for all u€Sj, (4.12) 

then 7 ^) = 0 . 


Proof. We use the symmetric tensor product 

1 


T fc (g) s = = V ^(x CT (i) ,-,x a m)'l! k - i (x a r e+ i),...,x CT t k )) 

y/l\{k~l)\k\^ k 

for two vectors NE 7 € and 'I'k-e € and assume that (14.1211 holds for all u £ $). 

Picking two unit vectors u, v, replacing u by u + tv in (14.1211 and taking the derivative 
with respect to t (which must be zero), we obtain 

x(u 9k - 1 <g> a v,'yWu 9k )=0. 

Next, replacing u by u + itv and taking the derivative with respect to t we get 

= 0 . 

Hence for all u and v 

u® k ~ 1 ® s v, 7 W u ® k ^ = 0. (4.13) 

Doing the same manipulations but taking now second derivatives with respect to t we also 

v, 7 ( fc ) u® k ~ l ® s + 2 (u® k ~ 2 ® s v® 2 ,^u® k ^ = 0. 

But, on the other hand, replacing u by u + tv and u + itv in (14.131 and taking first 
derivatives we get 

( 1 i® k ~ 1 ® s v,^ k) u® k ~ 1 

Combining the last two equations we infer 

u,7 (fc) u® fc_1 <g) s = 0 


have 


u ® k ~ 1 ® s 


i s v) + (u® k ~ 2 ® s v ® 2 , 7 ^) u® k ) = 0 . 


U® k ~ l ®s 


for all unit vectors u and v. Taking v of the form v = V\ ± v\ then v = v\ ± v6\ and 
iterating the argument, we conclude 

(vi ® s v 2 ® s ... ® s v k , 7 (fc) ui ® s v 2 ® s ■ ■ • v k ) = 0 

for all Vj,Vj £ S). Vectors of the form v\ ® s v 2 ® s ■ ■ ■ <S>s v k form a basis of Sj k , thus the 
proof is complete. □ 


For self-containedness we use the previous lemma to give a short proof of Schur’s for¬ 
mula (14.51 : 
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Proof of Lemma f.3. In view of Lemma 14.51 it suffices to show that 

[ I/.. _.\ 12JV 


dimly 


(Mi V) 


du = 1 


ISf) 


for all v € SSj. Pick v and v in SSj and U a unitary mapping such that Uv = v. Then, by 
invariance of du under rotations, 


dim I} 


N 


[ IMI 

JSf I 


2N du = dimfj^ 


/ \(Uu,Uv)\ du = diming / |(u, u)|“' du 

'ssj Jss) 


and thus the quantity 


N 


dim 


[ \(u,v)\ 2N du 
Jssj 


does not depend on v. By Lemma 14.51 this implies that 

dimf)^ f \u® N )(u® N \du = cL^jv 
Jss) 

for some constant c. Taking the trace of both sides of this equation shows that c = 1 and 
the proof is complete. □ 


In the sequel we use standard bosonic creation and annihilation operators. For all 
/ € Sj, we define the creation operator a*(f) : &s by 

a* if) ( X Mp ® - ® fa(k- 1 ) ) = ( k )~ 1/2 X /«r(i) ® ® /o-(fc) 

\aeSk -1 / o-eS/b 

where in the right-hand side we set fk = /. The annihilation operator a(/) : L ) fc+1 —> S) k 
is the formal adjoing of a*(f) (whence the notation), defined by 


a if ) ( X 0 - 0 /o-(fc+i) I = (* + i) 172 X] (/’ /fff 1 )) /«-(2) ® ® f*(k) 

yrr£Sk+i ) cr£Sk +1 

for all /, fi, ..., fk in fj. These operators satisfy the canonical commutation relations (CCR) 
[«(/), 0 ( 5 )] = 0, [a*(f),a*(g)] = 0, [a(/), a* ( 5 )] = (/, 5 ) 3 . (4.14) 

(n) 

One of the uses of these objects is that the reduced density matrices 7 ]^' of a bosonic 
state Tat are characterized by the relational 

(»®”, A"L®”) = ~" )! Tr^ [ 0 *W” 0 M”r„]. (4.15) 

Lemma 14.51 guarantees that this determines 7 ^ completely. The definition above is called 
a Wick quantization: Creation and annihilation operators appear in the normal order, all 
creators on the left and all annihilators on the right. 

The key observation in the proof of Theorem 14.41 is that the density matrices of the 
state (HU) can be alternatively defined from Tjy via an anti-Wick quantization where 
creation and annihilation operators appear in anti-normal order: All annihilators on the 
left and all creators on the right. 


18We recall our convention that Trpyj^] = 1. 
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Lemma 4.6 (The CKMR construction and anti-Wick quantization). 

Let f jv be defined by m and 7 ^ be its reduced density matrices. We have 

J + Z + [«(»)”«*(»)"I>] (4.16) 

for all v £ $). 


Proof. It suffices to consider the case of a pure state Tat = | th jv) iv | and write 
(v® k ,^v® k ) = dimi^ [ du\(u® N ,^ N )\ 2 \(u® k ,v® k )\ 2 


iss j 


= dim 15 
N\ 


N 


[ du\(u®^ N+k \v® k ®^ N )\^ 
Jssi 


(N + k)\ 


dimf) 


N 


'S 7 




IV! dirrifj;^ 

(N + k)\ dimi5f +fc 

(JV + <£ — !)! 


(a(u) fc 4'Ar,a(u) A: 4'Ar) 


- { N + k + d - 

using Schur’s lenrnra (14.51) in in the third line, the fact that a*(v) is the adjoint of 

a(v) in the fourth line and recalling that 

W + d - r 


dim ijl' = 


d- 1 


(4.17) 


the number of ways of choosing iV elements from d, allowing repetitions, without taking 
the order into account. This is the number of orthogonal vectors of the form 

u il ■ ■ ■ ®s Ui N , (ir,..., ijy) G {1 ,... ,d} 

one may form starting from an orthogonal basis (u\,... ,Ud) of Sj, and these form an 
orthogonal basis of S) N . □ 


The way forward is now clear: we have to compare polynomials in a* ( v ) and a(v) written 
in normal and anti-normal order. This standard operation leads to the final lemma of the 
proof: 


Lemma 4.7 (Normal and anti-normal order). 

Let v G Sf j. We have 

a{v) n a*(v) n = ^ for all n € N. (4.18) 

k =0 ' ' 


Proof. The computation is made easier by recalling the expression for the n-th Laguerre 
polynomial 


L n{x) = 

k =0 


t x e 


These polynomials satisfy the recurrence relation 

(n + l)L n+ i(x) = ( 2 n + 1 )L n (x) - xL n (x ) - nL n _ fix) 
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and one may see that (14.181) may be rewritten 


n 

a(v) n a*(v) n = Y J ^ k a*{v) k a{v) k (4.19) 

k =0 

where the c n ^ are the coefficients of the polynomial 

L n (x) := n\ L n (—x). 

It thus suffices to show that, for any n > 1, 

a{v) n+l a*(y) n+1 = a*(v)a(v) n a*(v) n a(v ) + (2 n + 1 )a(v) n a*(v) n - n 2 a(v) n - 1 a*(u) n - 1 . 

Note the order of creation and annihilation operators in the first term of the right-hand 
side: knowing a normal-ordered representation of a{v) n a*(y) n and a(v) n ~ 1 a*(v) n ~ 1 we 
deduce a normal-ordered representation of the left-hand side. 

A repeated application of the CCR (14.141) gives the relations 

a(v)a*(v) n = a*(v) n a(v) + na*(y) n ~ l 

a(v) n a*(v ) = a*{y)a(y) n + na{v) n ~ l (4.20) 

Then 

a*(v)a(v) n a*(v) n a(v ) = a(v) n a* (v) n+1 a(v) - na^^a* (v) n a(v) 

= a(v) n+1 a*(v) n+1 - (n + 1 )a(v) n a*(v) n 
- na(v) n a*{v) n + n 2 a(y) n ~ l a*{v) 71 - 1 , 

and the proof is complete. □ 


The final formula (14.81) is deduced by combining Lemmas 14.5114.Hi and 14771 with (14.151) . 
simply noting that for j < n 


Tr [a*(u) J a(u) J r 7 v] = 


N\ 


(- N-j)\ 
N\ 


-(vV^rflvW) 


—{v® n ,r$®i® n - j v® n ) 


( N-j)\ 

' ,.<Sin pO) 


(N-j)\ 


^{v® n ,T%> ® s 


The first two equalities are just definitions and the third one comes from the bosonic 
symmetry of v® n 
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5. Fock-space localization and applications 

We now turn to the first implication in the proof strategy (13.461) . We shall need to 
convert weak-* convergence of reduced density matrices into strong convergence, in order 
to apply Theorem 13.51 The idea is to localize the state T n one starts from using either 
compactly supported functions or finite rank orthogonal projectors. One may then work 
in a compact setting with 6 1 -strong convergence, apply Theorem 13.51 and then pass to 
the limit in the localization as a last step. More precisely, we will use localization in finite 
dimensional spaces, in order to show that the general theorem can be deduced from the 
finite-dimensional constructive proof discussed in the previous chapter. 

The subtlelty here is that the appropriate localization notion for a iV-body state (e.g. a 
wave-function € L 2 (R dN )) is more complicated than that one is used to for one-body 
wave-functions (R d ). One in fact has to work directly on the reduced density matri¬ 

ces and localize them in such a way that the localized matrices correspond to a quantum 
state. The localization procedure may lead to particle losses and thus the localizaed state 
will in general not be a Wbody state but a superposition of k -body states, 0 < k < N, 
i.e. a state on Fock space. 

The localization procedure we shall use is described in Section 15.21 We shall first give 
some heurtistic considerations in Section EH in order to make precise what has been 
said above, namely that the correct localization procedure in L 2 (R dN ) must differ from 
the usual localization in L 2 (R rf ). Section [5.31 contains the proof of the weak quantum de 
Finetti theorem and a useful auxiliary result which is a consequence of the proof using 
localization. 


5.1. Weak convergence and localization for a two-body state. 

The following considerations are taken from m■ Let us take a particularly simple 
sequence of bosonic two-body states 


tkn .— 1p n &) s (fin 




' 4>n + (f’n ® Ipn) €E L 2 (R 2d ) 


(5.1) 


with i\j n and (j) n being normalized in L 2 (R d ). This corresponds to having one patricle in 
the state i[) n and one particle in the state (f> n . We will assume 


— 0 , 

which ensures that ||\L n || = 1. Extracting a subsequence if need be we have 

weakly in L 2 (R 2d ) 

and the convergence is strong if and only if ||T|| = |||| = 1. In the case where some 
mass is lost in the limit, i.e. ||\I/|| < 1, the convergence is only weak. 

We will always work in a locally compact setting and thus the only possible source for 
the loss of mass is that it disappears at infinity [ 13611137111381 139] , A possibility is that 
both particles cj) n and il> n are lost at infinity 

Vvi —*• 0, cj) n — ^ 0 in L 2 (R 2d ) in n — >• oo, 

in which case 1 0 in L 2 (R 2rf ). In T 2 (R 2rf ) this is the scenario that is closest to the 
usual loss of mass in L 2 ( R rf ), but there are other possibilities. 
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A typical case is that where only one of the two particles is lost at infinity, which we 
can materialize by 

ip n —*• 0 weakly in L 2 (R d ), <p n —» tp strongly in L 2 (R 2d ). 

For the loss of mass of ip n one may typically think of the example 

Ipn = *P (• + Xn) (5-2) 

with \x n \ —>• oo when n —>• oo and ip say smooth with compact support. We have in this 
case 

->• 0 in L 2 (R 2d ) 

but for obvious physical reasons we would prefer to have a weak convergence notion en¬ 
suring 

© -j=<p. (5.3) 

In particular, since only the parrticle in the state ip n is lost at infinity it is natural that the 
limit state be one with only the particle described by cp left. We denote this convergence —*• g 
because this is precisely the geometric convergence discussed by Mathieu Lewin in HE]- 
The difficulty is of course that the two sides of (15.31) live in different spaces. 

To introduce the correct convergence notion, one has to look at the density matrices 
of 

7 = I (pn Ipn) ((pn ®s 1pn\ 0 in 6 1 (l 2 (R M )) 

7^ = h I'M (fin | + \ | Ipn) (lpn\ * \ \P) (P\ hi 6 1 (L 2 (R d )) . (5.4) 

One then sees that the pair ^7^* , 7^J converges to the pair (0, \ \ <p) (^>|) that corresponds 

to the density matrices of the one-body state y/2 1 <p € L 2 (R d ). More precisely, the 
geometric convergence notion is formulated in the Fock space (here bosonic with two 
particles) 

Tf 2 (L 2 (R d )) := C © L 2 (M d ) © L 2 (M m ) (5.5) 

and we have in the sense of geometric convergence on (5 1 (J~f 2 ) 

0 0 0© |T n ) (T n | -^ g i © i |^) ( (P\ 0 0, 

which means that all the reduced density matrices of the left-hand side converge to those 
of the right-hand side. We note that the limit does have trace 1 in S 1 (j 7 ^ 2 ), there is 
thus no loss of mass in J-pp 2 . More precisely, in J-f N the loss of mass for a pure 
IV-particles state is materialized by the convergence to a mixed state with less 
particles. 

Just as the appropriate notion of weak convergence for A r -body problems is different 
from the usual weak convergence in L 2 (R dN ) (a fortiori when taking the limit N —» 
oo), the appropriate procedure to localize a state and turn weak convergence into strong 
convergence must be thought anew. Given a self-adjoint positive localization operator A, 
say A = P a finite rank projector or A = x the multiplication by a compactly supported 
function y, one usually localizes a wave-function ip € L 2 (M rf ) by defining 

ip a = Aip 
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which amounts to associate 

\ip) (^| CA 1^4-0) (M\ ■ 

Emulating this procedure for the two-body state (EED one might imagine to consider a 
localized state defined by its two-body density matrix 

= | A®A^ n ){A®Am n \. 

It is then clear that 

In A - 1 ‘ * 0 in S 1 ^), 

which was to be expected, but it is more disturbing that we also have for the corresponding 
one-body density matrix 

InA 0 in S 1 ^) 
whereas, in view of (15.41) one would rather like to have 

InA {M\ strongly in S 1 ^). 

The solution to this dilemma is to define a localized state by asking that its reduced 
density matrices be A®A^^A®A and Aj^^ A. The corresponding state is then uniquely 
determined, and it so happens that it is a state on Fock space, as we explain in the next 
section. 

5.2. Fock-space localization. 

After the preceding heuristic considerations, we now introduce the notion of localization 
in the bosonic Fock soacJ^I 

•F s (£) = C © Sj © ... © fj™ © ... 

j;(L 2 (M d )) = C 0 L 2 (R d ) © ... © L 2 s (R dn ) 0 .... (5.6) 

In this course we will always start from IV-body states, in which case it is sufficient to 
work in the truncated Fock space 

J t ? N (S)) = C 0i50...0i5"0...0iof 

jf N (L 2 (R d )) = C ffi L 2 (R d ) 0 ... 0 L 2 (M dn ) 0 ... 0 L 2 (R dN ). (5.7) 

Definition 5.1 (Bosonic states on the Fock space). 

A bosonic state on the Fock space is a positive self-adjoint operator with trace 1 on Jj. 
We denote 5(77,.(if)) the set of bosonic states 

sfcm = {re = r*,r > o,T^ } [r] = 1 }. ( 5 . 8 ) 

We say that a state is diagonal (stricto sensu, block-diagonal) if it can be written in the 
form 

T = Go © Gi 0 ... 0 G n 0 ... (5.9) 

with G n € (S 1 ^™)- A state T^ on the truncated Fock space J-j^fj), respectively a 
diagonal state on the truncated Fock space, are defined in the same manner. For a diagonal 
state on j), of the form 

r = Go,TV 0 Gi t TV 0 ... © Gtv,TV, 

l^The procedure is the same for fermionic particles. 
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(n) 

its n-th reduced density matrix T K N ' is the operator on Sjf defined by 

r » = (T) 1 E (*) Tr » + i-* (5.10) 

^ ' k=n ^ ' 

□ 

The acquainted reader will notice two things: 

• We introduce only those concepts that will be crucial to the sequel of the course. 

One may of course define the density matrices of general states, but we will not 
need this hereafter. For a diagonal state, the reduced density matrices 115.101) 
characterize the state completely. For a non-diagonal state, one must also specify 
its “off-diagonal” density matrices i-a fjf for p / q. 

• The normalization we chose in (15.101) is not standard. It is chosen such that, in 
the spirit of the rest of the course, the n-th reduced matrix of a N- particles state 
(i.e. one with Go ,at = • • • = Gv-i,ac = 0 in (15.101) 1 be of trace 1. The standard 
convention would rather be to fix the trace at ( ), which is less convenient to 
apply Theorem 13.51 In other words the normalization takes into account the fact 
that throughout the course we work we a prefered particle number N. 

We may now introduce the concept of localization of a state. We shall limit ourselves 
to IV-body states and self-adjoint localization operators, which is sufficient for our needs 
in the sequel. The following lemma/definition is taken from [112] , Other versions may be 
found e.g. in [3l [55; [89]. 

Lemma 5.2 (Localization of a IV-body state). 

Let Tat € S(Sj^) be a bosonic N-body state and A a self-adjoint operator on 15 with 
0 < A 2 < 1. There exists a unique diagonal state Tjy G S(iFf N ($))) such that 

(r^) (n) = A® n T^ ] A® n (5.11) 

for all 0 < n < N. Moreover, writing Tjj in the form 

= Gq n © Gf N © ... © 
we have the fundamental relation 

Tr 5? [G)v,n] = Tr f D ^- TI 
Remark 5.3 (Fock-space localization). 

(1) The uniqueness part of the lemma shows that one has to work on Fock space. The 
localized state is in fact unique in S(J- s (S))), but to see this we would need slightly 
more general definitions, cf m- 

(2) The relation (15.121) is one of the cornerstones of the method. Loosely speaking it 

expresses the fact that, in the state T n, the probability of having n particles 
A-localized is equal to the probability of having N — n particles \/l — A 2 
localized. Think of the case of a very simple localization operator, A = 1^(0, R)> 
the (multiplication by the) indicative function of the ball of radius R. We are then 
simply saying that the probability of having exactly n particles in the ball equals 
the probability of having exactly N — n particles outside of the ball. Indeed, in 
probabilistic terms, these two possibilities correspond to the same event. □ 


G 


vT=54^ 
N,N—n 


(5.12) 
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Proof of Lemma HOI Uniqueness, at least amongst diagonal states on the truncated Fock 
space, is a simple consequence of the fact that the reduced density matrices uniquely 
characterize the state. Details can be found in [ 112] . 

For the existence, one can use the usual identification (in the sense of unitary equiva¬ 
lence) 

T s {AS) © \/l - A 2 S)) ~ F s (A$f) © J'si's /1 - A 2 Sj) 

and dehne the localized state by taking a partial trace with respect to the second Hilbert 
space in the tensor product of the right-hand side. We shall follow a more explicit but 
equivalent route. To simplify notation we shall restrict to the case where A = P is an 
orthogonal projector and thus y/l — A 2 = P±. 

First note that for any n-body observable O n 


Tr 


q p®np( n )p®n 


Tr 


p®riQ p®n-p( n ') 


Tr [(. p® n O n P ® n ) © l 9N ~ n T N ] 

Tr [O n © t® N ~ n (P® n © t® N ~ n T N P® n 




so that 

p®n r (n) p ® n = 0 0 _ 

Thus, by cyclicity of the trace, we obtain 


p>®7ip( n ) p®n 


— Tr^+i-^jy 
N-n 


p®n ^ -^®{N-n)^p®n ^ p®(]V-n)p^p®7i ^ p®(A-n) 

p®n+k ^ p(8>(AT—n—fc)p^p(g>n ^ p<g>(JV—7i) 


k =0 ' 
N-n 


E l N — ri\ 

( ) TV n+1 ^jv 

k =0 V 

J2 ( N U _ r ) ^n+l^IV \p® k ® P®( N ~ k) T N P® k © pf (JV ~ fc) 


p®n+k ^ p®(-/V—n—fc)p^p(g)n+fc ^ p®(N—n—k) 


k=n 


k — n 


Here we inserted 1 = P+P± in the first occurence of in the first line and expanded. 

Since the state P® n © x^f^-nlp iyP® n © i®(W-n) we act on is symmetric under exchange 
of the last N — n variables, we can reorganize the terms in the expansion to obtain the 
sum of the second line, involving binomial coefficients. It then suffices to note that 

N - n\ _ (N\ / N\ _1 fk 
k — n J \k J \n J \n 


to deduce 


N 


k=n 

with (cf Definition (15.101) 1 

'If 


-l 


n 


Tr n+ i^.7v [GN,k] — {Gn) 


p \ ( n ) 


n P 

^N,k ~ 


k 


TV, 


k+l-tN 


P ®k 0 p®( N ~ k )r N p* 


®fc ^ p®(N-k) 


(5.13) 
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and 

Gn = G%f } o © ■ ■ ■ © Gn > n - 

This is indeed an operator on P^ N {PS)): in (|5.13l) . one can interchange the first k particles 
and they live on PS), whereas the last N — k particles are traced out. 

There remains to show that G^ is indeed a state, i.e. that its trace is 1. To see this, 
we write 


i = t v[r*] = TV [(p + p ±f N {p + P±f N 

N 


E 

fc =0 


Tr 


-Sr 


P® k ig) p^ N ~ k ^Y N P® k <g> pf( N - k ) 


N 


= £>«* [G? jk ] =Tr m) [G^}. 

k =0 


The relation (15.121) is an immediate consequence of (|5.13D and the symmetry of Tjv- □ 


5.3. Proof of the weak quantum de Finetti theorem and corollaries. 

We are now going to use the localization procedure just described to prove the first 
implication of the strategy (13.461) . The idea is to use a finite-rank projector P, and then 
combine (15.101) with (15.111) to write (with n € N fixed) 

= £ f N ) 

k=n V 7 V 7 
N / , \ n 

V+i ^ kG N,k- (5-14) 

k=n ^ 7 


Here we inserted the simple estimate (see the computation in |113j , Equation (2.13)) 

^ / A C \ / K‘ \ 

(5.15) 


= (^) +o(iv- 1 ' 


We then argue as follows: The terms where k is small contribute very little to the 
sum ()5.14p because of the factor (-^) . For the terms where k is large, we note that, 
up to normalization, G^ k is a /e-particles bosonic state over PS). One may thus apply 
the de Finetti theorem discussed in Chapter []] to it, without worrying about compactness 
issues since PS) has finite dimension. Since k is large in these terms, and n is fixed, we 
obtain (formally) 


TWi-fcGft,* « T [ dMu)\u® n )( 

J uGSPS) 


for a certain measure v^, and thus 

N 


« E TVPCl (4) / 

} ~Ztr V iV / Ju&SPS) 


In the limit N —>■ oo, the discrete sum should become an integral in A = k/N. Using the 
fact that Gft is a state to deal with normalization (recall that Tr[G^ fc ] = 1), it is 
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natural to hope that we can obtain 

p®n^n)p®n^ [\ n [ du(X, u) (u) |u® n ) (u® n |, 

JO J u£SPS) 

which may be rewritten in the form (13.2411 by defining (in “spherical” coordinates on the 
unit ball BPS)) 

d/i(u ) = dfi (jM| , n—jj-^ := d ||u || 2 X dv^2 

There remains, as a last step, to apply this procedure for a sequence of projectors Pi —» 1 
and to check some compatibility relations to conclude. The final measure might not be a 
probability measure, which one can compensate by adding a delta function at the origin 
without changing any of the previous formulae. 

Note that this proof, which combines the methods of Chapter [4] and Section [5.21 gives 
a recipe to construct the de Finetti measure “by hand” (up to the fact that we have to 
pass to the limit at some point). This is very useful in practice, see Chapters [ 6 ] and 0 
The spirit of the proof using localization is reminiscent of some aspects of the method of 
Ammari and Nier mm who define cylindrical projections of the measure first. 

Let us now present the details of the proof, following [1131 Section 2]. 



Proof of Theorem 3.1 Existence. We carry on with the above notation. The following 
formula defines as a measure on [ 0 , 1 ], with values in the positive hermitian matrices 

of size dim (<g>™(PS))): 

N 

dM^(A) := ]T S k/N ( A) TWi-fcGft.fc. 


k=n 


We then have, using ([5.151) 


Tr 


/* 1 p N 

p®n 7 (n)p®n _ / yi d M^(A) < — E Tr G^ k -A 0 when N ^ 00 . (5.16) 

" / ° V k=n 


(n) 

Since P is a finite rank projector and 7 ]^' converges weakly-* by assumption, 

p®n^{n)p®n p®n^(n)p®n 

strongly in trace-class norm. On the other hand 
-l 


Ti f 


f dMp],(A) 
Jo 


N N 

E Tl V °N,k < E ^ G Pk = 1, 

k=n+l k=0 


(5.17) 


(n) 

so Mp N is a sequence of measures with bounded total variation over a compact finite 
dimensional space (positive hermitian matrices of size dim PS) having a trace less than 

1). One may thus extract from it a subsequence converging weakly as measures to Mp l \ 
Combining with (15.161) and (15.171) we have 

/*1 

p®n ry (n)p®n = / X n dMf\\). 

Jo 


(5.18) 
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An) 


neN 


that we just obtained 


We now have to show that the sequence of measures (M p 

is consistent in the sense that, for all n > 0 and for all continuous functions / over [0,1] 
vanishing at 0, 

^ /(A) Tr n+1 dM^ +1 \ A) = f 1 f(X)dM p '\x). (5.19) 

IQ Jo 

Here Tr n +i denotes partial trace with respect to the last variable. We have 


f 


N 

Tr n+1 dM < pJ f 1) (X) = Y, S k/N {X) Tr n ^ k G N , k 

k=n -\-1 


= dMgr(A) - 6 n/N (X)G^ n 


and thus 


f 1 /(A) Tr^n Tr n+1 dM^+ 1 \x)-dMP N (X) < f f(X)6 n/N (X) Tr fl 
Jo ’’ Jo 


u N,n 


0 

<f( — 
- J \n 


since Tr^n G^j n < 1. Passing to the limit we obtain (15.1911 for all function / such that 

/( 0 ) = 0 . 

We now apply Theorem 13.51 in finite dimension. Stricto sensu, this result applies only 

(n) 

at fixed A, but approaching dM p by step functions and then passing to the limit, we 
obtain a measure up on [0,1] x SS) n (PS)) such that 


f 1 f(X)dM < j?\x)= I [ l f(X)du P (X,u)\u® n )(u® n 
Jo Jssj Jo 

for all continuous functions / vanishing at 0. We thus obtain 


p®n rr Hp®n = / / dup(X,u)X n \u® n )(u® n \ 

Jo Jssj 

= 11 dup(X,u)\(VXu)® n )((VXu)® n \ 
Jo Jss) 


[ dnp(u)\u® n )(u* 
Jbs j 


defining the measure /rp in spherical coordinates. We are free to add a Dirac mass at the 
origin to turn /rp into a probability measure. 

The argument can be applied to a sequence of finite rank projectors converging to the 
identity. We then have a sequence of probability measures fj k on BS) such that 


p®n ry (n)p®n = f d ^ u ) \ u ®n^ u ®ny 

Jbs j 

Taking an increasing sequence of projectors (i.e. P k S) C P k+ \S)) it is clear that ji k coincides 
with on P^S) for £ < k (cf the argument for measure uniqueness below). Since all these 
measures have their supports in a bounded set, there exists (see for example [186 . Lemma 
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1]) a unique probability measurco h on BS) which coincides with / 1 k on P k S) in the sense 
that: 


/ dn k (u)\u® n )(u® n \= [ dfi(u)\(P k uf n )((P k uf 

IBS 3 JBS) 


We thus conclude 


P® n ^ n )p® n = p® n ( / d^{u)\u® n ){u® n \ ) PI 


'BSj 


and there only remains to take the limit k —» oo to deduce the existence of a measure 
satisfying (13.241) . 

Uniqueness. Let us now prove that the measure just constructed must be unique. Let 
fi and // satisfy 

[ \u® k )(u® k \dp(u)- [ \u® k ){u® k \dp'{u) = 0 (5.20) 

JBS) JBS] 

for all k > 1. Let V = vect(ei,...,e^) be a finite-dimensional subspace of fj, Pi = 
|ej)(ej| the associated projectors and fiy, p' v the cylindrical projections of p and // over V. 


Applying P h 


Pik on the left and Ph 


Pj k on the right to ()5.20|) . we obtain 


/ u n ■ ■ ■ u ik u h ■ • • u jk d{n v - A i' v ){u) = 0 
Jbv 


for all mutli-indices i\, ...,i k and ji, ■■■■,j k (we denote u = Y^ d j=\ u j e j )• O n the other hand, 


by 5 1 invariance of both measures it is clear that if k ^ 

Hi ■ ■ ■ u ik Uj[ ■ ■ ■ uf t d{fly - A iy)(u) = 0. 


f U, 

JBV 


Since polynomials in m and uj are dense in C°(BV., C) (continuous functions on the unit 
ball of V), we deduce that the cylindrical projections of n and [P over V coincide. That 
being true for all finite dimensional subspace V, we conclude that the two measures must 
coincide everywhere. □ 

We now give a particular case, adapted to our needs in the next chapter, of a very useful 
corollary of the above method of proof (see m\ Theorem 2.6] for the general statement): 

Corollary 5.4 (Localization and de Finetti measure). 

Let (r^AeM be a sequence of N-body states over Sj = L 2 (R d ) satisfying the assumptions 
of Theorem \ 3. 7| and p be the associated de Finetti measure. Assume that 


TV 




= Tr 


|V| 7 ff|V 


< c 


(5.21) 


for some constant independent of N. Let x be a localization function with compact support 
with 0 < x < 1 and G ^ be the localized state defined in Lemma 15. A Then 

N 


lim 

TV—>• oo 


E/( 


k =0 


k 
\N 


Ti' fD fc G^ k — 


IBS) 


dn(u) f(\\xu\\ 2 ) 


(5.22) 


for all continuous functions f on [0,1]. 


2 ^To construct it, note that the cr-closure of the union for k > 0 of the borelians of P^S) coincides with 
the borelians of $). 
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Proof. Since polynomials are dense in the continuous functions on [0,1] it is sufficient to 
consider the case /(A) = A” with n = 0,1,.... We then use (15.101) . (15.111) and (15.151) again 
to write 


N 


E(4)" Tr »*°^- Tr »- 


k =0 


X® n 7S‘ ) X® n 


0 when N —>• oo. 


Assumption (I5.21|) ensures 


TV 


(n) 


E ~ A i ] IN 
vl =1 


< Cn 


and we may thus assume that 


E ih I IE iv 


EiH- V”' EiH whenJV 


oo. 


,l=i 


i.l=i 




, 1=1 


The multiplication operator by x is relatively compact with respect to the Laplacian, thus 

Di := X s " I E IH ' 
is a compact operator. We then have 


TV, 


X® n 7S‘ ) X®” 


= Tr f 


TV, 


= Ti'f 




,i =i 


,i =i 


^ Ein h'"' Ein ^ 


,1=! 

(n) ®n 


,1 =1 


We conclude by using (I3.24P that 


N 

E 

fc =0 


Tr^fc 


JbSj 

= f dn(u) \\xu\\ 2n = [ d^u)f (\\ X 1 
J BSx Jb?) 


□ 


Remark 5.5 (Weak de Finetti measure and loss of mass.). 

(1) Assumption (|5.21l) is used to ensure strong compactness of density matrices, see the 
proof. It is of course very natural for the applications we have in mind (uniformly 
bounded kinetic energy). The convergence in (15.221) means that the mass of the de 
Finetti measure [i on the sphere {||u|| 2 = A} corresponds to the probability that a 
fraction A of the particles does not escape to infinity. 
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(2) We will use this corollary to obtain information on particles escaping to infinity in 
the following manner. We define the function 


V 


= n/i^ 


X^ 


which localizes “close to infinity”. One of course cannot apply the result directly 
to the localized state G V N . Instead one may use Relation (|5.12l) to obtain 


N 


lim S'' f ( — 

N—? OC ' X 


k =0 


k 


N 


Tr fi fc G v n k = lim 


N 

— V 

AT—>oo ' 
k =0 




Tr fi , Gl 


N,k 


[ dn(u ) / (1 - llx^ll 2 ) 

./ RS) 


(5.23) 


which gives some control on the loss of mass at infinity, encoded in the de Finetti 
measure. This is a bit surprising since the latter by definition describes particles 
which stay trapped. Typically we will use (15.231) with /(A) ~ en(A), the Hartree 
energy at mass A, see the next chapter. 

□ 
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6. Derivation of Hartree’s theory: general case 

We now turn to the general case of the derivation of Hartree’s functional as a limit of the 
iV-body problem in the mean-field regime. In Chapter [3 we saw that, under simplifying 
physical assumptions, the result is a rather direct consequence of the weak and strong de 
Finetti theorems. The general case requires a more thorough analysis and we shall use 
fully the localization tools introduced in Chapter [5J 

The setting is now that of particles in a non-trapping external potential V, interacting 
via a potential that may have bound states. Let us recap the notation: The Hamiltonian 
of the full system is 

N 1 

h n = ^2 T j + w ( Xt ~ ( 6 ‘ 1 ) 

j =1 1 1 <i<j<N 

acting on the Hilbert space f), with fj = L 2 (R rf ). The one-body Hamiltonian is 

T = - A + V, (6.2) 

that we assume to be self-adjoint and bounded below. We have emphasized the dependence 
on the potential V in the notation (16.11) because we will be lead to consider the system 
where particles are lost at infinity described by the Hamiltonian Hjy where one sets V = 0. 
One can generalize in the same directions as mentioned in Remark 13.21 but we shall for 
simplicity stick to the above model case. 

The interaction potential w : R i— > M will be asumed bounded relatively to T: for some 
0 < (3- , (3 + < 1 and C > 0 

- /3_(Ti + T 2 ) - C < w(x i - x 2 ) < /3+(Ti + T 2 ) + C. (6.3) 


We also assume symmetry 


w(—x) = w(x), 


and some decay at infinity 

w € L p (M d ) + L°°(R d ), max(l, d/2) < p < oo —>• 0, w(x) —> 0 when \x\ —> oo. 


(6.4) 


Again, it is rather vain to consider partially trapping one-body potentials, and we thus 
assume that V is non-trapping in all directions: 

V € L p (R d ) + L°°(R d ), max 1, d/2 < p < oo, V(x) —> 0 when |x| —> oo. (6.5) 

This ensures that Hn is self-adjoint and bounded below. The ground state energy 
of (13.11) is always given by 

E V {N) = inf<7 flJ vi$= inf (6.6) 

Finally let us recall what the limit objects are. The Hartree functional with potential V 
is given by 

4H : = / |Vu| 2 + V\u\ 2 + i [[ \u(x)\ 2 w(x - y)\u{y)\ 2 dxdy (6.7) 

JR d z J JR d xS. d 
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and we shall use the notation £^ for the translation-invariant functional where V = 0. 
The Hartree energy at mass A is given by 

eg (A) := inf £^ [it], 0 < A < 1. (6-8) 

||w|| 2 =A 

Under the previous assumptions we will always have the binding inequality 

en(l) < en(A) + en(l — A) (6.9) 

which is easily proved by evaluating the energy of a sequence of functions with a mass A 
in the well of the potential V and a mass 1 — A escaping to infinity. We will prove the 
following theorem, extracted from [113] (particular case of Theorem 1.1 therein). 


Theorem 6.1 (Derivation of Hartree’s theory, general case). 

Under the preceding assumptions, we have: 


(i) Convergence of energy. 


lim 

N —^oo 


E V (N) 

JT~ 



( 6 . 10 ) 


(ii) Convergence of states. LeC$>N be a sequence of L 2 (R dN )-normalized quasi-minimizers 
for HXj: 

($> N , H%V N ) = E V (N) + o(N) when N —»• oo, (6.11) 

and 'yff be the corresponding reduced density matrices. There exists /i € V(B£>) a proba¬ 
bility measure on the unit ball of S) with p(M v ) = 1, where 

M v =| u€BS) : £^[u] = e^(|M| 2 ) = e#(l) - e° H (l - |M| 2 )}, (6.12) 

such that, along a subsequence, 


Tn] [ \u® k )(u® k \dp,{u) 

3 J M v _ 

weakly-* in & 1 (Sj k ), for all k > 1. 


(6.13) 


( Hi ) If in addition the strict binding inequality 


e ti(l) < e H(^) + e H(l — 'M 


(6.14) 


holds for all 0 < A < 1, the measure pL has its support in the sphere SSj and the limit (16.131) 
holds in trace-class norm. In particular, if e^( 1) has a unique (modulo a constant phase) 
minimizer uh, then for the whole sequence 


7 jy \u^ k )(u^ k \ strongly in 6 1 (Sj k ) 

for all fixed k > 1. 


(6.15) 


Remark 6.2 (Characterization of the limit set). 

It is a classical fact that 

JCO' = e BSj | 3(u n ) minimizing sequence for ejj(l) such that u n —*■ u weakly in L 2 (K d 

This follows from the usual concentration compactness method for non-linear one-body 
variational problems. One can thus interpret (16.131) as saying that the weak limits for 
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the IV-body problem are fully characterized in terms of the weak limits for the one-body 
problem. □ 

The proof proceeds in two steps. In Section 16.11 we first consider the completely 
translation-invariant case where V = 0, which will describe particles escaping from the po¬ 
tential well V in the general case. We show that the energy e^(l) is the limit of N 1 E°(N). 
In this case one cannot hope for much more since there always exist minimizing sequences 
whose density matrices converge to 0 in trace-class norm. In addition to the tools already 
introduced we shall rely on ideas of Levy-Leblond, Lieb, Thirring and Yau [1111 113211133] 
to recover a bit of compactness. 

We then use fully the localization methods of Chapter [5] to treat the general case 
in Section 16.21 In the spirit of the concentration-compactness principle we will localize 
minimizing sequences inside and outside of a ball. The inside-localized part is described 
by the weak-* limit of reduced density matrices and we can thus use the weak de Finetti 
theorem. The outside-localized part no longer sees the potential V and we may thus apply 
to it the results of Section [6T1 


6.1. The translation-invariant problem. 

Here we deal with the case where V = 0. It is then possible to construct quasi¬ 
minimizing sequences for E°(N) with yjy (. — yjy) 0 for any sequence of translations 
xn. One may thus have vanishing in Lions’ terminology [ 136 , 137] and without any specific 
trick one cannot hope for more than the convergence of the energy. Indeed, it is possible 
to construct a state where the relative motion of the particles is bound by the interaction 
potential w but where the center of mass vanishes. This implies vanishing for the whole 
sequence. We shall thus be content with proving the convergence of the energy: 


Theorem 6.3 (Energy of translation-invariant systems). 

Under the preceding assumptions, we have 


to *£Q = 4 ( 1 ). 

N^-oo N 


(6.16) 


The first step is to use part of the interaction potential to create an attractive one-body 
potential (this roughly amounts to taking out the center of mass degree of freedom). This 
way we will define an auxiliary problem whose energy is close to the original one, but 
for which particles always stay trapped. This is done in the following lemma, inspired 
from [T32lfT33|: 

Lemma 6.4 (Auxiliary problem with binding). 

We split w = w + — w~ into positive and negative parts and define, for some e > 0 

w £ (x) = w(x) + ew~(x). 

Consider the auxiliary Hamiltonian 

N 1 

h n = Y { Ki ~ £W -( X i)) + Y W e{xi-Xj) ( 6 . 17 ) 

i =1 1 l<i<j<N 

and E £ (N) the associated ground state energy. Then 

E e (N ) , E°(N) 

a £ := lim ——— < lim ——— 

N—too N N— >oo N 


= : a. 


(6.18) 
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Proof. Using symmetry, we can for all T' € Lj(\ R“ JV ) write 


N 


N 


/N—l 




^ 1=1 


Similarly 
2 

N(N-l) , 

V \ l<i<j<N 

with 7 ^ = |T)('I'|. Then 


T, X w ( x i ~ x i / = Tr f 


( 2 ) 


, 2—1 


= TV, 


( 2 ) 


— eTr, 


- ( 2 )' 

W 7^, 


TV, 


and 


( 2 )' 

W £ %' 


eTr, 


(N-l)(N-2) 


TV 


x w £ (xi~xj) 7 $ 


- ( 2 ) 


N — l 


■TV 


'TV-1 


X W (Xi- X N )~jq, 


2—1 


All this implies 


N - 1 
N 


(*,<*) 


' N—l 


2 — 1 


T', X ( - Aj - eu; (37 - x N )j + 


)) + ^ 7 ^ X] - x j) ) ^ V (6-19) 

“ l<i< 7 <Af-l / / 


In the preceding equation the Hamiltonian between parenthesis depends on xjy via the 
one-body potential ew~{xi — xn) but since the other terms are translation-invariant the 
bottom of the spectrum is in fact independent of xn- We thus have 

^L^±(*,H 0 n V)>E £ (N- 1)<T,T) (6.20) 

for all T e L^(M. dN ), which implies 

E°(N ) ^ E e (N - 1) 

N ~ N — l 


The sequences E°(N)/N and E £ (N)/N are increasing since they are the infimum of vari¬ 
ational problems set on decreasing sets (cf (II. 551) 1. Using trial states spread over larger 
and larger domains method, one may on the other hand show that 


E°(N) 

N 


< 0 , 


E £ (N ) 
N 


< 0 


and thus that the limits a and a £ exist. Then (16.201) clearly implies (16.181) . □ 


We next derive a lower bound to E £ (N ) for e small enough. It is much easier to work on 
this energy because the corresponding sequences of reduced density matrices are strongly 
compact in S 1 . Indeed, compactness or its absence (physically, binding or its absence) 
resuts from a comparison between the attractive and repulsive parts of the one- and two- 
body potentials. Here the one-body potential ew~ and the two-body potential w £ are 
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well equilibrated because they have been built precisely for this purpose, starting from the 
original two-body potential w. This trick will allow us to conclude the 

Proof of Theorem 1 6'. -71 As usual, only the lower bound is non-trivial. Since ejj(l) < 0 one 
may assume that a < 0 since otherwise a = e^( 1) = 0 and there is nothing to prove. We 
are going to prove the lower bound 


\u(x)\ 2 w £ (x - y)\u{y)\ 2 dxdy \ , ( 6 . 21 ) 


a £ > efj(l) := inf { (u, (-A - ew-)u) + ^ 

IMI 2 =i l 1 

and it is then easy to show that efj(l) —>• e^(l) when e —>• to obtain (16.161) by combining 
with (16.181) . 

Let be a sequence of wave-functions such that 


{*N, H £ n V n ) = E £ (N) + o{N). 


(k) 

and 7^7 the corresponding reduced density matrices. Then 


a F = lim = Hm j ^ 




TV— >oo N TV—voo 

After the usual diagonal extraction we can assume that 

v O) 


+ 2 ^ 


( 2 ) 

WsT N 


(fc) 

Tn 


T 


and we are going to show that the convergence is actually strong. We shall afterwards use 
the strong de Finetti theorem to obtain Hartree’s energy as a lower bound. 

We pick a smooth partition of unity Xr + 77 ^ = 1 and use Lemma 13.121 to get 


a £ > lim inf lim inf l Tr f , (-A - ew^xr'Jn’xr 

ix —>00 N—> oo 1 


+ i TVs. 


®2 (2) <g)2 

w eXR 7 n Xr 


+ Tr^ 


-Arms 1 nVR 


+ 2 TV„ a 


Igi2 (2) ig)2 

w sVr in Vr 


. ( 6 . 22 ) 


We define the xr~ and r?_R-localized states and G V N by applying Lemma mi to ’Fat. 
We will use those to estimate separately the two terms in the right side of (16.221) . 

The XR-localized term. Using (15.101) and (15.111) we have 


Tr 


7 


(-A - £W-)xrXn ] XR 


+ 2 dry >2 


(g>2 (2) (g>2 

WeXR 7 nXr 


1 


TV 


N 


J2 Tr X k ( ^(- A - ew -)i + J2 W ^ Xi ~ x i) I G N,k 
k= 1 yi=l i<j J 

We apply the inequality 

A + tB = (1 — t)A + t(A + B) > (1 — t) inf cr(A) + t inf a{A + B ) 


(6.23) 


4 = E<- A 


— A — EW _ 


A + B = H £)k , t=(k-l)/(N-l). 


1=1 


with 


(6.24) 
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We have 

and since we assume that a < 0, for e small enough 


lim inf cr(—A — EW-) = inf a{— A) = 0 
£—^0 


Thus 

and we may then write 
Tr* 


inf er(—A — ew-) > a> a £ > k 1 inf a(H e ,&). 
inf a (A) > inf a (A + B ) 

N k Tr G^j k E %k) 


1 1 


(-A - ew-)xrT n Xr\ + 2 Tr 7 2 


<S>2 (2) 02 

w £ Xr 7 n Xr \ >2^ 
k =1 


N k 


But since 


N k Tr (l x 

^2 K 11 


k =0 


N 


= Tr 


2 ( 1 ) 
XrT n 


. TV—>>00 


Tr 


Xi7 (1) 


and lim = a £ 

k —^oo k 


and k H > —4r^ is increasing we conclude 


lim inf Tr* 

N^oo ' 


(—A — ew^xr^nXR 


+ g Tr 7 2 


02 ( 2 ) ®2 
WcXr 7 n Xr 


> a £ Tr 


Xr7 (1) 


by monotone convergence. 

The T^-localized term. Using the results of Section 15.21 as above we have 


(6.25) 


Tr* 


Tvr^nVR 


+ 2 Tr 7 2 
N 


02 ( 2 ) 02 
w sVr In Vr 




k =1 


E-*< + vbX>(*< 

. *=1 i<j 


— Xi 


3) I ^N,k 


. (6.26) 


Here we use 


—A >0, w £ = w + 2 ew > (1 — 2 e)w and E°(k ) < ak < 0 


to obtain 
k 


J2~ A ' + wiri^2 w ^- x 3) ^ 


N — 1 

i =1 i<j 


> 


(1 - 2e)(fc - 1) 
aTTT Hk 

(l-2E)(k-l) 


N — l 


E°(k) > K\k) - 2sak 


for all k > 1. Combining with (16.261) we obtain 


Tr fl -A •nRljNVR 
We finally deduce that 


+ 2 Tr » 1 


02 ( 2 ) 02 
w sVr Tn Vr 


> 


N 

E 

k =1 


^TrG^ fc /£*>(*) 


A V k 


— 2 ECL 


lim inf Tr* 

TV—»oo ' 


TrjR^Nm 


H— Tr <-,2 
2 ’ 


02 ( 2 ) 02 

7Ar ; % 


> (l-2e)a(l—Tf[ X | 7 (1) ]) (6.27) 
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upon using 

N 


E 

k =0 


^ T 'G« = T '' 


2 ( 1 ) 
VrTn 


. TV—>-oo 


1 - Tr 


xh {1) 


and 


E°(k) 

Inn —-— = a 

k— »oo k 


E°(k ) 


as well as the fact that k i->- —is increasing. 
Conclusion. Inserting (16.2511 and (16.271) in (16.221) we find 


Op- > lim inf ( a F Tr 


xW l) 


+ (1 — 2e)a(l — Tr 


R —>oo 

= o £ TV[ 7 «] + (1 - 2e)a(l - TV[ 7 «]) 
Since we assumed a e < a < 0, we obtain 

TV[ 7 (1) ] = 1. 


2 ^(i) 


XrI 


(6.28) 


There is thus no loss of mass for the auxiliary problem defined in Lemma [6.41 This suffices 
to conclude that the reduced density matrices 7 jy converge strongly. Indeed, one may 
apply the weak de Finetti theorem to the sequence of weak limits 7 ^ to obtain a measure 
fi living on the unit ball of Sj. But ()6.28l) combined with (13.241) implies that the measure 
must in fact live on the unit sphere. We thus have, for all k > 0, 

Tr[ 7 «] = 1, 


which implies that the convergence of 7 jy to 7 ( fc ) is actually strong in trace-class norm. 
We can then go back to (16.221) to obtain 


lim inf 

TV—>■ oo 


Tr, 


(-A - eu;_) 7 £ ) 


+ 2^ 2 


( 2 ) 

Wel N 


> Try, 


(—A — ew -) 7 ^^ 




Wel 


( 2 ) 


We then apply the strong de Finetti theorem to the limits of the reduced density matrices 
to conclude that the right side is necessarily larger than efj(l). This gives (16.211) and 
concludes the proof. □ 


6.2. Concluding the proof in the general case. 

We have almost all the ingredients of the proof of Theorem 16.11 at our disposal. We 
only need a little bit more information on the translation-invariant problem, as we now 
explain. 

For k > 2 consider the energy 


hW ■= ^ inf Ufik 


A 


k 

2=1 


iE 

Kj 


W[Xj — Xj 


(6.29) 


i.e. an energy for k particles with an interaction strength proportional to t~[- In the last 
section we have already shown that the limit k —>• oo of such an energy is given by Ae^(A) 
when A is a fixed parameter. Using Fock-space localization methods, the energy of particles 
lost at infinity in the minimization of the energy for N > k particles will be naturally 
described as a superposition of energies for fc-particles systems with an interaction of 
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strength 1 /{N — 1) inherited from the original problem. In other words, we will have to 
evaluate a superposition of the energies b k (\) with 

\ fe " 1 

IV- 1’ 

a bit as in (16.231) and (16.261) . Since this A depends on k it will be useful to know that b k { A) 
is equicontinuous as a function of A: 

Lemma 6.5 (Equi-continuity of the energy as a function of the interaction). 

We set the convention 

6o(A) = &i(A) = 0. 

Then, for all A € [0,1] 

lim A6fc(A) = ejj(A). (6.30) 

k—¥ OO 

Moreover, for all 0 < A < A' < 1 

0 < b k { A) - b k ( A') < C\X - A'| (6.31) 

where C does not depend on k. 

Proof. We start by vindicating our claim that (16.301) is a direct consequence of the analysis 
of the previous section. For A = 0 there is nothing to prove. For A > 0 we use Theorem l6.3l 
to obtain 

lim A b k (X) = A inf ((u, Ku) + ^ ff w(x - y)\u(x)\ 2 \u{y)\ 

IMI =i V 2 J J 

= || |[^1 A ((“> Ku ^ + \ JJ w ( x ~ y)l M ( x )l 2 l ?/ (y)l 2 ) = 4(A)- 

The fact that 

b k ( A) > b k {\') for all 0 < A < \' < 1 

is a consequence of (16.241) . Indeed, one can see that either b k ( A) = 0 for all A or the same 
kind of argument as in (16.241) applies. 

For the equicontinuity (16.311) we fix some 0 < a < 1 and remark that, with 

S:= (A' — A)(a -1 — A) -1 , 

we have 


\*=1 i<j ) 


\-b 


A 




i =1 
k 


k - i ^ Wij 

Kj 


+ i^\ a Y, K ' + —Y. 


Wi 


i=l 


k~ 1 ^ ~ lJ 

i<j 


> (1 - S)b k ( A) - — 

a 

using the fact the spectrum of the operator appearing in the second line is bounded below. 
We deduce 

0 < b k ( A) - b k ( A') < 5{b k { A) + CoT 1 ) < C\\' - A| 
since b k { A) is uniformly bounded and |<5| < C |A — A'|. □ 
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We may know conclude the 

Proof of Theorem \6.1[ Let ^>n be a sequence of iV-body wave-functions such that 

(* N ,H%* N ) = E v (N) + o(N), 

(k) 

with 7 k n ' the associated reduced density matrices. After a diagonal extraction we have 


(fc) 

Tn 


7 


(*) 


weakly-*. Theorem 13.71 ensures the existence of a probability measure /r such that 

d/r(u)|u® fc )(u® A; |. (6.32) 


7 ( fc ) = 


' uGBS) 
.2 


We pick a smooth partition of unity x# + r) 2 R = 1 as previously and define the localized 
state and G V N constructed from |vJLv)(\lLv|. Using Lemma f3. 12 1 again we obtain 

E V (N ) 


lim 

TV—>-oo 


N 


= lim ( TV 

TV—>-oo 






+ 2 Tr » s 


( 2 ) 

wTn 


> lim inf lim inf < Tr^-, 

R-^OO TV—>-oo 


Txr^nXR 


+ g 


ig)2 (2) (g>2 

w Xr InXr 


+ TV 


7 


n ] vr 


T — Tr^ 


Igi2 (2) K>2 

7jv 7 r 


(6.33) 


First we use the strong local compactness and (16.321) for the xiR-localized term: 


lim inf < Try 

N —»oo I 


Txr^nXR 


+ j T W 


(g)2 (2) <g>2 

7]v ar 


> Try-, 


7\'R7 (1) XR 


+ 2 T r 7 2 


02 ( 2 ) 02 


£n[XRu]dfj,(u). (6.34) 


' BSj 


Our main task is to control the second term in the right side of (16.331) . We claim that 


lim inf ( Tr 

TV—>-oo 


Tvrx'nVR 


+ 2 T V 


(g>2 (2) ig)2 

Tn or 


> 


[ e^(l - (6.35) 

J RSj 


Indeed, using the ^-localized state GJj we have 


TVf 


-A Vr^nVR +2 Tr ^ 2 


(g)2 (2) ig)2 

w Vr Tn Vr 


N 




fc=i 


E- a ‘ + wtE 




i=l 


TV - 1 ^ ~ 13 

Kj 


N f k 

^E^ TrG ^ inf ^ E- a * + f3t 

fc=l \i =1 



AT 


- E Tr 

fc =0 


fc- 1 


_IV- 1, 

where is the function defined in (16.291) . On the other hand 


N 


lim V TrGV, -b 

oo N,k \ N 


k =0 


k 


N 


k- 1 
IV - 1 


-4|4))=o, 


WVA 


(6.36) 
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because, using the equicontinuity of {6fc}^ =1 and the convergence lim*,-^ A6&(A) = e^(A) 
given by Lemma 16.51 we get 


lim sup 

N^i-oc k=l,2,...,N 

We just have to combine this with 


-h 

N 


N 


k- 1 
N- 1 


k 


= 0. 


E TrG E = 1 


k =0 


to obtain (16.361) . At this stage we thus have 


lim inf Tr 

TV—>-oo 


io 


-A Vr^nVR + 2 Tr fl 2 


682 (2) 682' 
'TiV Vr 


N 

k =0 


We now use the fundamental relation (15.121) and Corollary 15.41 as indicated in Section T5.31 
to deduce 


N 


lim VTrG^ 

AT—>00 2-^ 


k =0 


k 


** 6H [n =n^o 


N 


^12 TrG N,N-k e H 


k =0 


( X ~ 7 ^) = J m e n ( 1 - ||xfl«l| 2 )rf/*(«)» 


which concludes the proof of (|6.35l) . 

There remains to insert (|6.34l) and (16.351) in (16.331) and use Fatou’s lemma. This gives 

e„(U> n m 


> lim inf 

R-*-°o \JbSj l 


> 


L 


> 


IBS D L 


£r[xru] + eui 1 ~ \\xru\\ 2 ) dfi(u] 
[xru] + en(l - IIxr^II 2 ) dfx(u) 

[u] + en(l — |M| 2 ) dfi(u) 
e u(\\ u \\ 2 ) + ®h(! - IM | 2 )1 d n(u) > e H (l), 


lim inf 

BSj R ~^°° L 

f cV r„.i , „o /i ii,.n 2 

Jbs i L 


(6.37) 


using the continuity of A i —> e^(A) and the binding inequality (16.91) . This concludes the 
proof of (16.101) . and the other results of the theorem follow by inspecting the cases of 
equality in (16.371) . □ 
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7. Derivation of Gross-Pitaevskii functionals 


We now turn to the derivation of non-linear Schrodinger (NLS) functionals with local 
non-linearities: 

£*#]==/ |vV’| 2 + P|V’l 2 + fH 4 - 

One may obtain such objects starting from a Hartree functional such as (13.811 and taking 
an interaction potential converging (in the sense of distributions) to a a Dirac mass 


w L (x) = L d w —>■ (/ v?J when L 


0 . 


Since we already obtained (13.81) as the limit of a IV-body problem, one might be tempted 
to see the derivation of the NLS functional as a simple passage to the limit in a one-body 
problem. The issue with such an approach is the total lack of control on the relationship be¬ 
tween the physical parameters N and L. One can see this as a problem of non-commuting 
limits: it is not clear at all that one can interchange the order of the limits N —> oo anco 
L —» 0. 

From a physical point of view, the good question is “What relation should N and L 
satisfy if one is to obtain the NLS energy by taking simultaneously the limits N —>• oo and 
L —>• 0 of the IV-body problem ?” The process leading to the NLS theory is thus more 
subtle than that leading to the Hartree theory. We shall first elaborate a little bit more 
on this point in the following subsection. 


7.1. Preliminary remarks. 

Until now we have worked with only two physical parameters: the particle number N 
and the interactions strength A. In order to have a well-defined limit problem we have 
been lead to considering the mean-held limit where A oc IV . In this case, the range of the 
interaction potential is fixed and each particle interacts with all the others. The interaction 
strength felt by a typical particle is thus of order XN oc 1, comparable with its self-energy 
(kinetic + potential). We saw that this equilibration of forces acting on each particle, 
combined with structure results a la de Finetti, naturally leads to the conclusion that 
particles approximately behave as if they were independent. It follows that Hartree-type 
descriptions are valid in the limit N —>• oo. 

There are other ways to justify such models: the equilibration of forces which allows the 
limit problem to emerge may result from a more subtle mechanism. For example, in the 
ultra-cold alkali gases where BEC has been observed, it has more to do with the dilute¬ 
ness of the system than with the weakness of the interaction strength. For a theoretical 
description of this situation, we can introduce in our model a length scale L characterizing 
the range of the interactions. Taking the total system size as a reference, a dilute system 
is materialized by taking the limit L —>• 0. The interaction strength for a typical particle 
is then of order \NL d (interaction strength x number of particles in a ball of radius L 
around a given particle). This parameter is the one we should fix when N —>• oo to obtain 
a limit problem. Different regimes are then possible, depending on the ratio between A 
and L. 


21 The L —¥ 0 limit of the IV-body problem is by the way very hard to define properly. 
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One may discuss the different possibilities starting from the /V-body HamiltoniarF^l 
N 1 

H N = Y / -^ + V(x J ) + ^-j N^wiN^Xi-Xj)) (7.1) 

j =1 1 <i<j<N 

which corresponds to choosing 

L = N~P, A oc N d P~ l . 

The fixed parameter 0 < f3 is used to set the ratio between L and A. We consider the 
reference interaction potential w as fixed, and we shall denote 

wn{x) := N d ^w(N^x). (7-2) 

For f3 > 0, wn converges in the sense of distributions to a Dirac mass 


w N ->■ 


iL w ) io ’ 


(7.3) 


materializing the short range of the interactions/diluteness of the system. Reasoning 
formally, one may want to directly replace wn by ( f Rd w) <5o in (17. lj) . In this case we are 
back to a mean-field limit, with a Dirac mass as interaction potential. This is of course 
purely formal (except in ID where the Sobolev injection H 1 C° allows to properly 

define the contact interaction). Accepting that this manipulation has a meaning and that 
one may approximate the ground state of (17.11) under the form 

V n = i/P n (7.4) 

we obtain the Hartree functional 

£#]==/ \V^\ 2 + V\^\ 2 + h^\\w*\^\ 2 ). (7.5) 

JM. d ^ 

Replacing the interaction potential by a Dirac mass at the origin leads to the Gross- 
Pit aevskii functional 

£ nls [ip]:= f \VM 2 + V\M 2 + ^\M A - (7.6) 

In view of (17.31) . the logical choice seems to imagine that when f3 > 0, we obtain, starting 
from 4m the above functional with 

a = f w. (7.7) 

J R d 


In fact one may prove the following results (described in the case d = 3; d < 2 leads to 
subtleties and d > 4 does not have much physical meaning): 

• If /3 = 0 we have the previously studied mean-field (MF) regime. The range of 
the interaction potential is fixed and its strength decreases proportionally to N _1 . 
The limit problem is then (17.51) . as previously shown. 


2 ^Once again, it is possible to add fractional Laplacians and/or magnetic fields, cf Remark 13.21 For 
simplicity we do not pursue this here. 
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• If 0 < P < 1, the limit problem is 117.61) with the parameter choice (17.71) . as can be 
expected. We will call this case the non-linear Schrodinger (NLS) limit. This case 
does not seem to have been considered in the literature prior to [ 115 : but, at least 
when w > 0 and V is confining, one may adapt the analysis of the more difficult 
case P = 1 . 

• If P = 1, the limit functional is now (17.61) with 

a = 4-/T x scattering length of w 

(see i125 . Appendix C] for a definition). In this case, the ground state of (17.11) 
includes a correction to the ansatz (EU), in the form of short-range correlations. 
In fact, one should expect to have 

N 

^ N (xi,...,x N ) « / [N p {xi- Xj)') (7.8) 

j =1 l<i<j<N 

where / is linked to the two-body problem defined by w (zero-energy scattering 
solution). It so happens that when P = 1, the correction has a leading order effect 
on the energy, that of replacing f Rd w by the corresponding scattering length, as 
noted first in [59]. We will call this case the Gross-Pitaevskii (GP) limit. It is 
studied in a long and remarkable series of papers by Lieb, Seiringer and Yngvason 
(see for example [ESI Q351 HZ] E26] US] E23] H25] h 
As already mentioned, the corresponding evolution problems have also been thoroughly 
studied. Here too one has to distinguish the MF limit IT). |5] [76] [1651 1157] . the NLS 
limit [65] 1156] and the GP limit [69: .1551 43 ] . 

There is a fundamental physical difference between the MF and GP regimes: In both 
cases the effective interaction parameter XNL 3 is of order 1, but one goes (still in 3D) 
from a case with numerous weak collistions when A = A r_1 , L = 1 to a case with rare but 
strong collisions when A = N 2 , L = N . The different NLS scaling sort of interpolate 
between these two extremes, the transition from “frequent weak collistions” to “rare strong 
collisions” happens at A = 1,L = A r_1 / 3 , i.e. P = 1/3. 

The difficulty for obtaining (17.61) by taking the limit N —>• 00 is that there are in fact two 
distinct limits N —>• 00 and wn —>• a5 0 to control at the same time. A simple compactness 
argument will not suffice in this case and one has to work with quantitative estimates. The 
goal of this chapter is to explain how one may proceed starting from the quantitative de 
Finetti theorem stated in Chapter [U Since this theorem is only valid in finite dimension, 
we will need a natural way to project the problem onto finite dimensional spaces. We thus 
deal with the case of trapped bosons, assuming that for some constants c, C > 0 

c\x\ s -C <V(x). (7.9) 

In this case, the one-body Hamiltonian — A + P has a discrete spectrum on which we have 
a good control thanks to Lieb-Thirring-like inequalities. 

The results we are going to obtain are valid for 0 < P < Pq where Pq = Po(d, s ) depends 
only on the dimension of the configuration space and the potential V. We shall give 
explicit estimates of Pq, but the method, introduced in [ 115| . is for the moment limited to 
relatively small p. In particular, we will always obtain (17.71) as interaction parameter. The 
main advantage compared to the Lieb-Seiringer-Yngvason method (1341 il26t 11281 il23j is 
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that we can in some cases avoid the assumption w > 0, always made in these papers (see 
also [125J ). In particular we present the first derivation of attractive NLS functional^ in 
ID and 2D. 

Recently, a blend of the following arguments, the tools of |[13lt (126111281 il23j and new 
a priori estimates for the many-body ground state has allowed to extend the analysis to 
the GP regime |149 j. One can also adapt the method to deal with particles with statistics 
slightly deviating from the bosonic one (“almost bosonic anyons”), cf m- 


7.2. Statements and discussion. 

We take comfortable assumptions on w: 

w € L°°(R d ,R) and w(x) = w(—x). (7-10) 

Without loss of generality we assume 

sup |u>| = 1 

m d 

to simplify some expressions. We also assume that 

x eA (1 + \x\)w(x) € L 1 (M <i ), (7-11) 

which simplifies the replacement of wjy by a Dirac mass. As usual we use the same notation 
wn for the interaction potential 117.211 and the multiplication operator by wn(x — y) acting 
on L 2 ( R m ). 


For /3 


0 we have shown previously that 


lim 

N—> oo 


E(N) 

N 


= e H - 


(7.12) 


We now deal with the case 0 < (3 < /3o(d, s) < 1 where we obtain the ground state energy 
of 1(7151) : 

e n i s := inf £ n ls W\ (7.13) 

II Z , 2 (K rf ) 1 

with a defined in (17.71) . Because of the local non-linearity, the NLS theory is more delicate 
than Hartree’s theory. We shall need some structure assumptions on the interactions 
potential (see [1151 for a more thorough discussion): 

• When d = 3, it is well-known that a ground state for (17.61) e xists if and only if a > 0. 
This is because the cubic non-linearity is super-critical 24 ! in this case. Moreover, 
it is easy to see that N~ 1 E(N ) —>• —oo if w is negative at the origin. The optimal 
assumption happens to be a classical stability notion for the interaction potential: 



y)p(y)dxdy > 0, for all p € L l (M d ), p 


> o. 


(7.14) 


This is satisfied as soon as w = w± + u> 2 , w\ > 0 and W 2 > 0 with W 2 the Fourier 
transform of W 2 - This assumption clearly implies f Rd w > 0, and one may easily see 
by changing scales that if it is violated for a certain p > 0, then E(N)/N —» —oo. 


09 

One often uses the non-linear optics vocabulary to distinguish the repulsive and attractive cases : 
repulsive = defocusing, attractive = focusing. 

24 One may for example consult ll02j for a classification of non-linearities in the NLS equation. 
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When d = 2, the cubic non-linearity is critical. A minimizer for (17.61) exists if and 
only if a > —a* with 

a* = \\Q\\ 2 L 2 (7.15) 

where Q € H 1 (M 2 ) is the unique [108] (modulo translations) solution to 

-AQ + Q-Q 3 = 0. (7.16) 

The critical interaction parameter a* is the best possible constant in the interpo¬ 
lation inequality 


[ \u\ 4 <c([ |Vu|' 
Jr 2 \J r 2 


(7.17) 


We refer to |88 l ll42j for the existence of a ground state and to [196| for the 
inequality (17.171) . A pedagogical discussion of this kind of subjects may be found 
in [Sj- 

In view of the above conditions, it is clear that in 2D we have to assume 
f w > —a*, but this is in fact not sufficient: as in 3D, if the interaction po¬ 
tential is sufficiently negative at the origin, one may see that N~ 1 E(N) —>• —oo. 
The appropriate assumption is now 

IMlia ||Vtt||| 2 + \ [[ \u(x)\ 2 \u(y)\ 2 w(x-y)dxdy > 0 (7.18) 

Z J J R 2 xR 2 

for all u € i7 1 (M 2 ). Replacing u by Au(Ax) and taking the limit A^Owe obtain 


Ml l 2 II v «IIl 2 + 2 


w 


f \u(x)\ 4 dx > 0, V«€i7 1 (M 2 ), 

Jr 2 


which implies that 


/ 

Jr* 


w(x) dx > —a* 


A scaling argument shows that if the strict inequality in (17.181) is reversed for a 
certain u, then E(N)/N —>• —oo. The case where equality may occur in (17.181) is 
left aside in these notes. It requires a more thorough analysis, see [88j where this 
is provided at the level of the NLS functional. 

• When d = 1, the cubic non-linearity is sub-critical and there is always a minimizer 
for the functional (17.61) . In this case we need no further assumptions. 

We may now state the 

Theorem 7.1 (Derivation of NLS ground states). 

Assume that either d = 1, or d = 2 and (17.181) holds, or d = 3 and (17.141) holds. Further 
suppose that 


0 < /? < Po(d,s) := 


2 ds + sd 2 + 2d 2 
where s is the exponent appearing in ClHD- We then have 
(1) Convergence of the energy: 

E(N) 


< 1 . 


(7.19) 


^nls* 


N 


(7.20) 
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(2) Convergence of states: Let 4 'n be a ground state of (17.11) and 


7^ :=TY n+1 _> iV [|\I'iv>(*iv|] 


its reduced density matrices. Along a subsequence we have, for all n € N 


lim = f dfi{u)\u &n ) (u® n 

N^-OO J\A , 


wGA^nls 


(7.21) 


strongly in 6 1 (L 2 (W ln )). Here fa is a probability measure supported on 


-Mnls — |w £ i (K ), ||u||^2 — l>£nlsM — 6nlsj 1 • 


(7.22) 


In particular, when (ESJ) has a unique minimizer u n \ s (up to a constant phase), 
we have for the whole sequence 




(7.23) 


Uniqueness of u n \ s is ensured if a > 0 or \a\ is small. If these conditions are not 
satisfied, one can show that there are several minimizers in certain trapping potentials 
having degenerate minima [9] [88] . 

Remark 7.2 (On the derivation of NLS functionals.). 


(1) The assumption /3 < /3o{d , s ) is dictated by the method of proof but it is certainly 


not optimal, see [1231U2B1H331H3511139] . One may relax a bit the condition on /?, 
at the price of heavier computations, something we prefer to avoid in these notes, 
see [115j . In ID, one may obtain the result for any /3 > 0. 

(2) It is in fact likely that one could replace /3o(d,s) by f3(d, oo) for any s > 0, using 
localization methods from Chapter [5] Indeed, ground states always essentially live 
in a bounded region. We refrain from pursuing this in order not to add another 
layer of localization (in space) to the proof: we will already rely a lot on finite 
dimensional localization. In any case, /3 < (3(d,oo) is not optimal either. 

(3) Let us inspect in more details the conditions on /3o(d, s ) we obtain. For the case of 
a quadratic trapping potential V(x) = |x | 2 for example, we can afford /3 < 1/24 in 
3D, (5 < 1/12 in 2D and /3 < 1/4 in ID. The method adapts with no difficulty to 
the case of particles in a bounded domain which corresponds to setting formally 
s = oo. We then otain /3o(d, s) = 1/15 in 3D, 1/8 in 2D and 1/3 in ID. Improving 
these thresholds in the case of (partially) attractive potentials remains an open 
problem. 

(4) When /3 is smaller than the critical (3o(d,s ) by a given amount, the method gives 

quantitative estimates for the convergence (17.201) . see below. We refer to [ 115 . 
Remark 4.2] for a discussion of the cases where a convergence rate for the minimizer 
can be deduced, based on tools from EDO and assumptions on the behavior of 
the NLS functional. □ 

The proof of this result occupies the rest of the chapter. We proceed in two steps. 
The bulk of the analysis consists in obtaining a quantitative estimate of the discrepancy 
between the IV-body energy per particle N~ 1 E(N ) and the Hartree energy 



(7.24) 
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given by minimizing the functional 

£hN := / [\S7u\ 2 + V\uf] dx+ ]- [[ \u{x)\ 2 w N (x - y)\u(y)\ 2 dxdy. (7.25) 

JR d v J l J j]Rd X Rd 

The latter objects still depend on IV when /3 > 0, whence the necessity to avoid compact¬ 
ness arguments and obtain precise estimates. Once the link between N~ 1 E(N) and (|7.24l) 
is established, there remains to estimate the difference |e n i s — en|, which is a much easier 
task. Most of the restrictive assumptions we have made on w are used only in this second 
step. The estimates on the difference |en — N~ 1 E(N)\ are valid without assuming (17.1111 
and (|7.14l) or (17.1811 . They thus give some information on the divergence of N~ 1 E(N) in 
the case where en does not converge to e n i s : 


Theorem 7.3 (Quantitative derivation of Hartree’s theory). 
Assume (17.91) and (17.1011 . Let 

1 + 2 d/3 


t : = 


2 + d/2 + d/s 


If 

t > 2d/3, 

then, for all d> 1 there exists a constant C d > 0 such that 

E(N ) 


eH > 


N 


> e n - C d N 


-t+2d/3 


(7.26) 

(7.27) 

(7.28) 


Remark 7.4 (Explicit estimates in the mean-field limit). 

(1) Condition (17.271) is satisfied if 0 < j3 < f3o(d,s). For the proof of Theorem 17. II we 
are only interested in cases where |en| is bounded independently of N, and (17.2811 
then gives non-trivial information only if (17.2711 holds. 

(2) The result is valid in the mean-field case where (3 = 0 and thus en does not depend 
on N. We then obtain explicit estimates improving on Theorem l3.6l These present 
a novelty in the case where Hartee’s functional has several minimizers, or a unique 
degenerate minimizer. In other caset-Fl. better estimates are known, with an error 
of order IV -1 given by Bogoliubov’s theory [1171 176 . 1861 56!. See ( 150] for exten¬ 
sions of Bogoliubov’s theory to cases of mutiple and/or degenerate minimizers. 

□ 


The proof of Theorem 17.31 occupies Section 17.31 We then complete the proof of Theo¬ 
rem o in Section PI 


7.3. Quantitative estimates for Hartree’s theory. 

The main idea of the proof is to apply Theorem 14.11 on a low-energy eigenspace of the 
one-body operator 

T = —A + V 

acting on Sj = L 2 (R d ). Assumption (17.91) ensures that the resolvent of this operator is 
compact and thus that its specturm is made of a sequence of eigenvalues tending to infinity. 

^The simplest example ensuring uniqueness and non-degeneracy is that where w > 0 with w the Fourier 
transform of w. 
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We denote P_ and P + the spectral projectors corresponding to energies respectively below 
and above a given truncation A: 

P_ = 1 ( _oo i A)(T), P+ = l fl -P_=Pi. (7.29) 

We also denote 

N\ := dim (P_.fi) = number of eigenvalues of T smaller than A. (7.30) 

Since the precision of the quantitative de Finetti theorem depends on the dimension of 
the space on which it applies, it is clearly necessary to have at our disposal a convenient 
control of N\. The tools to achieve this are well-known under the name of Lieb-Thirring 
inequalities, or rather Cwikel-Lieb-Rosenblum in this case. We shall use the following 
lemma: 


Lemma 7.5 (Number of bound states of a Schrodinger operator). 

Let V satisfy m- For all d> 1, there exists a constant C d > 0 such that, for all A large 
enough 

N\ < C d A d/s+d/2 . (7.31) 

Proof. When d > 3, this is an application of [1241 Theorem 4.1]. For d < 2, the result 
follows easily by applying [43J Theorem 2.1] or [1851 Theorem 15.8], see [ 115] for details. 
The familiarized reader can convince herself that the right side of (17.311) is proportional 
to the expected number of energy levels in the semi-classical approximation. We refer 
to [1241 Chapter 4] for a more thorough discussion of this kind of inequalities. □ 


In the sequel we shall argue as follows: 

(1) The eigenvectors of T form an orthogonal basis of L 2 (M. d ) on which the N particles 
should be distributed. The methods of Chapter [5] provide the right tools to analyze 
the repartition of the particles between P_f) and P + fj. 

(2) If the truncation A is chosen large enough, particles living on excited energy levels 
will have a much larger energy per unit mass than the Hartree energy we are 
aiming at. There can thus only be a small number of particles living on excited 
energy levels. 

(3) Particles living on P—Sj form a state on J-f N (P^S)) (truncated bosonic Fock space 
built on P-Sj). Since P_fj has finite dimension, one may use Theorem 14. II to de¬ 
scribe these particles. These will give the Hartree energy, up to an error depending 
on A and the expected number of PC-localized particles. More precisely, in view 
of (14. 2[) . we should expect an error of the form 


(A + N*P) x N a 

n1 


(7.32) 


i.e. dimension of the localized space x operator norm of the projected Hamiltonian 
/ number of localized particles. 

(4) We next have to optimize over A, keeping the following heuristic in mind: if A 
is large, there will be many P_-localized particles, which favors the denominator 
of (17.321) . On the other hand, taking A small favors the numerator of (17.321) . 
Picking an optimal value to balance these two effects leads to the error terms of 
Theorem 17.31 
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Proof of Theorem \7.3\ The upper bound in (17.281) is as usual proven by taking a trial state 
of the form u® N . Only the lower bound is non-trivial. We proceed in several steps. 

Step 1, truncated Hamiltonian. We first have to convince ourselves that it is legitimate 
to think only in terms of P + and P_-localized particles, as we did above. This is the oject 
of the following lemma, which bounds from below the two-body Hamiltonian 

H 2 = T®1 + 1®T + w N (7.33) 

in terms of its P_-localization, P_ <® P_ H 2 P_ <8 P- and a crude bound on the energy of 
the P+-localized particles. 


Lemma 7.6 (Truncated Hamiltonian). 

Assume that A > CN d & for a large enough constant C > 0. Then 

H 2 > p® 2 H 2 P® 2 + ^-P ® 2 - 

Proof. We denote 

H$ = T®1 + 1®T 

the two-body Hamiltonian with no interactions. We may then write 

h 2 = (P- + P+f 2 H°2 (P- + P+f 2 = E H ° 2 P fe ® p z 

i,j,k, ££{—,+} 

= ^ Pi ® Pj ^2 Pi ® Pj ■ 

i,3 6{-,+} 


Indeed, 

Pi <8 Pj H 2 Pk <8 Pc = 0 

if i ^ k or j ^ £, since T commutes with P± and P_P + = 0. We then note that 


P+TP+ > A P + and P_TP_ > -CP-, 

which gives 

P + <8 P + P + <8 P + > 2A P + <8 P + 

P + ®> P_ H 2 P + ® P_ > (A - C)P + <8 P_ 

P_ (8 P + P_ <8 P + > (A — C)P + (8 P_. 

We conclude 

H 2 > P® 2 P® 2 + (A - c)n 

where 

n = P+ (8 P+ + P_ <8 P+ + P+ <8 P- = lf-,2 - P® 2 . 
We turn to the interactions: 

w N = (p® 2 + n) w N (p® 2 + n). 

We have to bound from below the difference 


w N - P® 2 w n P® 2 


by controling the off-diagonal terms 

Uw N P® 2 + P® 2 w n n 


(7.34) 

(7.35) 


(7.36) 


(7.37) 


(7.38) 
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in (17.3811 . To this end we write 

wn = w~x — with > 0 

and we use the well-known fact that the diagonal elements of a self-adjoint operator 
control the off-diagonal elements. 

Since are positive as multiplication operators on L 2 (M? d ) we have for all b > 0 

(V/2 p ®2 ± ft- 1 / 2 if) w% (lj x / 2 P® 2 ± b-^uj > 0. 

Combining these inequalities appropriately we obtain 

n iv n p® 2 + p® 2 w N n > -bP® 2 \w N \p® 2 - 6 _1 n|i«jv|n 
for all b > 0. We then recall that, as an operator, 

I| < ||whv||L°° < N d/3 
and we choose b = 2N d @/A to obtain 

o ]\[2d/3 A 

Uw n p ® 2 + p® 2 w n n > ——p® 2 - f n 

Inserting this bound in (17.381) we get 

w N > P® 2 w n P® 2 - ^-J—P- 2 - + N d P\ n (7.39) 

by simply bounding from below 

nuiTvii > -n|u;7v|n > -N df} n. 

Combining with (17.371) and (17.391) we obtain for all A > 1 the lower bound 

h 2 > p® 2 h 2 p® 2 - ^L-^p® 2 + ^ - N dfi - c^j n 

Since we assume A > CN d @ for a large constant C > 0, we may use P_ <g) P_ < 1^2 and 
P_ <g> P + , P + (g) P_ > 0 to deduce 

9 /u2 d/3 a 

P 2 > P® 2 H 2 P® 2 - ^Pf 2 , 

which concludes the proof. □ 

Step 2, estimating the localized energy. Let Tat be a minimizer for the IV-body 
energy, T N = I^at) <^at| and 

7jV^ = Tr n+ i^.Ar[r7v] 

the corresponding reduced density matrices. We can now think only in terms of the P_ 
and P+-localized states defined as in Lemma 15.21 by the relations 

(G^) W = Pf^PfL (7.40) 

2 ®Cf for a positive hermitian matrix (rriij)i<i,j< n the inequality 2| m,ij\ < 
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We recall that are states on the truncated Fock space, i.e. 

N 

E Tr «‘[ G w] = 1 - < 7 - 41 > 

k =0 

We now compare the Hartree energy en and the localized energy of Tat defined by the 
truncated Hamiltonian of Lemma 17.61 


Lemma 7.7 (Lower bound for the localized energy). 

If A > CN d d for a large enough constant C > 0 we have, 


Tr 


p® 2 h 2 p® 2 2n 


I A rp 

+ I Tl ' 


p®2 (2) 

r \~ In 


>en~C 


AN± 

N 



(7.42) 


This lemma is proved by combining Theorem 14.11 and the methods of Chapter [5j We 
define an approximate de Finetti measure starting from Gl. The idea is related to that 
we used for the proof of the weak de Finetti theorem in Section f5.3l 


N 


dgL N {u) = ^ 


k =2 


-1 


dim ((P-Sj)ty (u® k ,G 


<S>k /'-y— 


T N,k L 


du 


(7.43) 


where du is the uniform measure on the finite-dimensional sphere SP-S). The choice of the 
weights in the above sum comes from the fact that we want to approximate the localized 
two-body density matrix P® 2 ^^ P®' 2 , which is the purpose of 


Lemma 7.8 (Quantitative de Finetti for a localized state.). 

For all A > 0 , we have 


Tr jT 


P 


®2 ( 2 ) p < 8>2 
- In - 


L 


SP-S) 


u® 2 )(u® 2 \d^ N ( u) 


< 81Va 

- ~JT' 


Proof. Up to normalization, G Nk is a state on (P-S))f k . Applying Theorem 14.II with the 
explicit construction (14.611 we thus have 


Tr^2 


Tr 3 ^ fc 



u &2 )(u® 2 \dfi N:k (u ) 


ISP-Si 


N A 

< 8—f- Tr- 



with 

dnN,k(u) = dim (P-S)) k (u® k ,Gj fk u® k 'jdu. 
In view of (I7.4UD and (17.431) we deduce 


Tic 


p®2^(2)p(g)2 _ 


I SP-S) 


\u® 2 ){u® 2 \dn N (u) 


N IN'” 1 

<E 


k=2 


k\ 8 N A 
2) — ^ 


G 


N,k 


There remains to use the normalization (17.411) and 


-l 


k(k — 1) < k 


N(N -1) - N 


to conclude the proof. 


□ 


We proceed to the 
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Proof of Lemma | 7. 7\ We start with the TL-localized term. By cyclicity of the trace we 
have 


Tr 


p® 2 h 2 p^^ 


<g>2 (2) 


= Tr 


We then apply Lemma 17.81 which gives 


Tr ft 2 


p®2 p^p®2 ^p®2^(2)p®2\ 

8N\ 


P® 2 ^ 2) P® 2 - [ \u® 2 )(u® 2 \dp N (u) 

JsP-Si 


< 


N 


On the other hand we of course have 

\\P® 2 H 2 P® 2 \\ < 2A + 11uijv11 l°° < 3A 
in operator norm, and thus 


-Tr 

2 


P®2 H2 P® 2 ^ 2) 




'SP_7 


Tr# [H 2 \u® 2 )(u® 2 \] dfr N - 


CAN a 

N 


I SP-Sj 


£n[u\dp N - 


CAN a 
N ' 


By the variational principle £u[u] > en, we deduce 


-Tr 

2 


P® 2 H 2 P® 2 7 (2) > eH ^ 


N 'N N 


k =2 


Tr flfc (G 


T N,k 


CANa 

N 


where the computation of f dp n is straightforward using Schur’s formula (14.511 . 
For the P + -localized term we use (|7.40l) . (15.1011 and (15.12|) to obtain 


A 


jTt[P® 2 7 <?*P? 2 ] = TE 


A 


N 


k=2 


a AT-2 

Ay 

4 2- 

k =0 

Gathering (17.4611 , (17.4711 and recalling that 


- 1 /,a k 2 


N 2 


+ 0(N 


-i A 


N 

N 

2 

-l 


-i 


-l 


Tr 


N — k 


G~N,k 


Tr 


G 


N,k 


N-k 

2 


(N - kf 
N 2 


+ 0(N 


-i\ 


we find 


-Tr 

2 

N 

>£ t.- 

k =0 


p®2 F2 p®2 7 (2)j + ^ Tr [p®2 7 (2)j 


G 


N,k 


P_ 

N 2 


en + 


A 

Tr[j 
(N-k) 2 A 


N 2 


C(|e H | + A) CANa 


N 2 


N 


(7.44) 


(7.45) 


(7.46) 


(7.47) 


(7.48) 


The first error term comes from the fact that the sums in (17.4611 and (17.4711 do not run 
exactly from 0 to N, and we have used the normalization of the localized states (17.41 II to 
control the missing terms. 

It is easy to see that for all p, q, 0 < A < 1 

2 

p\ 2 + q(l — A) 2 >p-y. 












































102 


NICOLAS ROUGERIE 


We then take p = en, Q = A/4, A = k/N and use (17.411) again to deduce 
TV 


p®2^ 2 p®2 7 (2)j + ^ TV(p ®2 7 (2) ) 


> e H - -r - 


C(|e H | + A) CAN\ 


A IV 2 

from (17.481) . There remains to insert the simple estimate 

M < C + ||waHU~ < C + < C + A 

to obtain the sought-after result. 


N 


(7.49) 

□ 


Step 3, final optimization.. There only remains to optimize over A. Indeed, we recall 
that by definition 

E(N) 1 ( 2 ) 1 

=2 Tr fl 2 [^27V] 

with the two-body Hamiltonian (17.331) . Combining Lemmas 17.61 and 17.71 we have the lower 
bound 

E(N) ^ _ CN 2d P „AN a „ A 


> e H - 


-C- 


N ° N 2 


N ~ A 

for all A > CN d ^ with (7 > 0 large enough. Using Lemma 17.51 this reduces to 

E(N) CN 2d>3 „ ^i+d/s+rf/2 


IV 


> e H - 


Optimizing over A we get 
where 


t = 2s 


A ~ C “- 
A = N t 
1 + 2 d/3 


IV 


(7.50) 


4s + ds + 2 d 

and the condition t > 2 d/3 in (17.271) ensures that A N 2d P for large N. We thus conclude 

E(N) 


en > 


N 


> e H - C d N 


-t+2d/3 


which is the desired result. 


□ 


Remark 7.9 (Note for later use.). 

Following the steps of the proof more precisely we obtain information on the asymptotic 
behavior of minimizers. More specifically, going back to ([7.451) . using Lemma 17.61 and 
dropping the positive P+ localized terms we have 

e H > ^ J Tr[P? 2 H 2 P® 2 ^} + o(l) > [ S H [u]d» N (u) + o(l) 

1 JSSj 

and thus 

o(l) > [ (£n[u] - en) dp N (u) (7.51) 

Jssj 

when N —>• oo. We do not specify here (cf Item (3) in Remark 17.21) the exact order of 
magnitude of the o(l) obtained by optimizing in Step 3 above. Estimate (17.511) morally 
says that pjy must be concentrated close to the minimizers of £h- This will be of use in 
the proof of Theorem 17.11 □ 




























DE FINETTI THEOREMS AND BOSE-EINSTEIN CONDENSATION 


103 


7.4. From Hartree to NLS. 

There remains to deduce Theorem o as a corollary of the above analysis. We start 
with the following lemma: 

Lemma 7.10 (Stability of one-body functionals). 

Consider the functionals (17.2511 and (17.611 . Under the assumptions of Theorem \7.1\ there 
exists a minimizer for £ n \ s . Moreover, for all normalized function u G L 2 (R. d ) we have 

IIMIIhi <C(Sh[u] + C) (7.52) 

and 

2 


\Sn[u] -£ nl sMI < CN-P( 1+ / |V«| 


Consequently 


|^H ^nls | — C I\ 


(7.53) 


(7.54) 


Proof. The stability assumptions we have made guarantee that minimizing sequences for 
fnis are bounded in H 1 . Assumption (17.311 allows to easily estimate the difference between 
the Hartree and NLS functionals for H 1 functions. Details are omitted and may be found 
in 1 1151 . □ 

We are now equiped to complete the derivation of NLS functional. 

Proof of Theorem 1 7 ,1\ Combining (17.5411 with (I7.28|l concludes the proof of (17.2011 . There 
thus remains to prove convergence of states, the second item of the theorem. We proceed 
in four steps: 

Step 1, strong compactness of reduced density matrices. We extract a diagonal 
subsequence along which 


(n) 

Tn 


7 


(n) 


when N ^ oo, for all n € N. On the other hand we have 


TV 


Of 


= TV 


(i) 


(~A + V)tN 


< C. 


(7.55) 

(7.56) 


independently of N. To see this, pick some a > 0, define 
N 


HN,a ~ ( — Aj + V(Xj )) + 


1 “h Q. 

N-l 




j =1 1 

and apply Theorem 17.31 to this Hamiltonian. We find in particular H^, a > —CN and 


deduce 


6nls T o(l) ^ 


(Tjy, H n ^ n ) 
N 


7 —0(1 + a) + 


a 


1 -\- OL 


Tr [T'y 


»'] 


which gives (17.56(1 . Since T = —A + V has compact resolvent, (17.55(1 and ((7.5611 imply 
that, up to a subsequence yjy strongly converges in trace-class norm. As noted previously, 
Theorem 13.51 then implies that also "yff strongly converges, for all n > 1. 

Step 2, defining the limit measure. We simplify notation by calling r/v the best 
bound on | E(N)/N — e n i s | obtained previously. Let d/i n be defined as in Lemma 17.81 It 
satisfies 

H N {SP-fj) = Tr [p® 2 7 ^ } P® 2 . 
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We have 


Tr 


J3®2^(2)pi§)2 _ 


ISP-sj 


\u® 2 )(u® 2 \dfj,N(u) 


8 N a Ai+d/«+d/2 
< —^ < C --- > 0 . 


N 


N 


One may on the other hand deduce from the energy estimates of Section 17.31 a control on 
the number of excited particles: 


1 - MSP-ft) = Tr [(1 - P? 2 hN 


< 


r/v 

A ' 


By the triangle and Cauchy-Schwarz inequalities we deduce 


Ti¬ 


ls) 

IN 


ISP-Sj 


\u® 2 ){u® 2 \d[lN{u) 


A \-\-df s-\-d /2 



r_N_ 

A 


(7.57) 


(7.58) 


We now denote Pj< the spectral projector of T on energies below a truncation K, defined 
as in (17.291) . Since 7 ^ and Pk —> 1 

lim lim hn(SPkSj) = 1 . 

K—> 00 TV—>-oo 

This condition allows us to use Prokhorov’s Theorem and [186 > Lemma 1] to ensure that, 
after a possible further extracion, hn converges to a measure ji on the ball BS j. Passing 
to the limit, we hnd 


7 ^ = f \u® 2 )(u® 2 \dn(u) 
JBSi 


and it follows that fi has its support in the sphere SS) since Tr^ 2 )] = 1 by strong conver¬ 
gence of the subsequence. 

Step 3, the limit measure only charges NLS minimizers. Using (17.511) and 


m (SP_i5) = l + 0 (^), 


we deduce that 


'SP-Si 


(£hM - eu)dn N (u) < o(l) 


in the limit N —>• 00 . By the estimates of Lemma 1 7.1(11 it follows that, for a large enough 
constant B > 0 independent of N, 

B 2 


C 


L 


Vu\\ L 2 >B 


djJLN{u) < 


||Vu|| L a >B 


[£h[u] - en)d/j, N (u) < o(l), 


and 


L 


(fnisM - e n\s)d/JN < C{ 1 + B A )N ,3 + (£ n [u] - eu)d/j N (u) < o(l)- 

..^ u \\l2<B J \[\7u\\ l 2<B 

Passing to the limit N —>• oo, we now see that fi is supported in A/l n i s . 

At this stage, using (|7.58|) and the convergence of /j ,n we have, strongly in trace-class 
norm, 
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where /x is a probability measure supported in -Ad n is- Taking a partial trace we obtain 

|xx)(xx|d/x(ix) 






and there only remains to obtain the convergence of reduced density matrices of order 
n > 2. 

Step 4, higher order density matrices. We want to obtain 

f \ u ® n )(u® n \dv(u), 

^Mnls 

in trace-class norm when N —>■ oo. In view of the definition of /x, it suffices to show that 


Tr 


[ \u® n )(u® n \d m (u) 

■) SPSS 


(7.59) 


where /xat is the measure defined by applying Lemma [T751 to T N ■ To this end we start by 
approximating using a new measure, a priori different from hn = ■ 

N /atN -l /) X 

dim {P_fy k s (u® k , G~ Nk u® k ^du. (7.60) 


d^ n N {u) = 

k=n 


Proceeding as in the proof of Lemma 17.81 we obtain 


Tr, 


p®n^(n) p®n _ 


'SPSS 


\u® n )(u® n \d^ N (u) 


An estimate similar to (I7.57|) next shows that 


Tr, 


(rc) 

Tn 


ISP-Sj 


\u® n )(u® n \dv n N (u) 


< C 


0. 


nN A 

N ' 


(7.61) 


Using again the bound 


-i 


n 


= (£] +°( JV 




as well as the triangle inequality and Schur’s formula ()4.5f) we deduce from (17.6111 that 


Tr 


(n) 

Tn 


'SPSS 


\u® n )(u® n \dfi N (u) 


N 


< 


k 


IT n 

N 


E((v) -(v) 

k =0 
n—1 

/ k- \ " 

G~ f 


G 


N,k 


k =0 
2 

+ E 

k =0 


k 




N 


T N,k 


Tr, 


G 


N,k 


+ o(l). 


(7.62) 


Finally, combining the various bounds we have obtained 


N ' 2 


k =2 


E(v) 


G 


N,k 


1. 
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But, because of (17.411) it follows that 

N 

E 

k =0 


f k\ 2 

r 

\n) Tr ^ 

G N,k 


i. 


Recalling the normalization 

N 

^ G Nk 

k =0 

we may thus apply Jensen’s inequality to obtain 


= 1, 


N 




k =0 


lfe\" 

n) 


G 


IV,A: 


> 


< N 

E 

\k=0 


N 


n/1 


Tr 




G 


IV, fc 


1. 


There only remains to insert this and (17.411) in (|7.62l) to deduce ([7.59D and thus conclude 
the proof of the theorem. □ 
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Appendix A. A quantum use of the classical theorem 

This appendix is devoted to an alternative proof of a weak version of Theorem 13.61 The 
method, introduced in [100] is less general than those described previously. This cannot 
be helped since it consists in an application of the Hewitt-Savage (classical de Finetti) 
theorem to a quantum problem. Here we follow an unpublished note of Mathieu Lewin 
and Nicolas Rougerie jllSj . 

In some cases (absence of magnetic fields essentially), the wave-function minimizing 
a A-body energy can be chosen strictly positive. The ground state of the quantum prob¬ 
lem may then be entirely analyzed in terms of the A-body density p^ N = |\kjv| 2 , which 
is a purely classical object (a symmetric probability measure) whose limit can be de¬ 
scribed using Theorem 12. 11 This approach works only under assumptions on the one-body 
Hamiltonian that are much more restrictive than those discussed in Remark 13.21 


A.l. Classical formulation of the quantum problem. 

We here consider a quantum A-body Hamiltonian acting on L 2 (W lN ) 

N 1 

h n = E ( T i + v ( x j)) + w ( Xi ~ (A- 1 ) 

j =1 1 <i<j<N 

where the trapping potential V and the interaction potential w are chosen as in Section 
In particular we assume that V is confining. We shall need rather strong assumptions on 
the kinetic energy operator T. The approach we shall discuss in this appendix is based on 
the following notion: 


Definition A.l (Kinetic energy with positive kernel). 

We say that T has a positive kernel if there exists T(x, y ) : W l xM^-) M + such that 

(ip, Tip) = I[ T(x,y) \ip(x) — ip(y)\ 2 dxdy (A.2) 

J JR d xR d 

for all functions ip € L 2 (M. d ). □ 

It is well-known that the pseudo-relativitic kinetic energy is of this form. Indeed 


(ip,V^A ^ = 


T(^±i) ff l'ip(x) — ip(y)\ 2 


dxdy , 


(A-3) 


2vr (n+l)/2 JJ RdxRd \x - y 

see mi Theorem 7.12]. More generally one may consider T = |p| s , 0 < s < 2: 

w.x)-m? ixdv 

/ Rd x «d \x-y\ d+s 

recalling the correspondance (11.2411 . 

The non-relavistic kinetic energy does not fit in this framework, but one may neverthe¬ 
less apply the considerations of this appendix to it because 


{fp, bl s > = (b, (-A) s/2 b) = c d , s Jj r 


(ip,-A ip) = / | Vg/>| 2 = C d lim (1 — s) 

J W d S H 


\ip(x) - ip(y)\ 4 
\x — y\ d + 2s 


■dxdy 


(A.4) 


with 


C d = 


cos 6 da 


-l 


S d-i 
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Here 5 rf_1 is the euclidean sphere equiped with its Lebesgue measure da and 6 represents 
the angle with respect to the vertical axis. One may thus see —A as a limiting case of 
Definition IA. 11 Formula ()A.4D is proved in [2H Corollary 2] and |143| . see also [25, TTl . 

The cases with magnetic field T = (p + A) 2 and T = |p + A| are not covered by this 
framework. One may deal with them with the methods of the main body of the course, 
as already mentioned, but not with those of this appendix. 


An important consequence of the above choice of kinetic energy is that, by the triangle 
inequality 

(^ry,)>(h/,|,T|y,|) 

and thus the total IV-body energy 




satisfies 

<?at[Tat] > £n [|^n|] • 

The ground state energy can thus be calculated using only positive test functions 

E(N) = inf {S N [V N }, Tat G ^(M dJV )} = inf {s N [^ N }, € L 2 s (R dN ), > o} . (A.5) 

This remark actually allows to prove a fact mentioned previously: the bosonic ground state 
is identical to the absolute ground state in the case of a kinetic energy of the form (IA.2I) 
or ([A.41) . see 1371 Chapter 3]. 

We recall the definition of Hartree’s functional: 


£ h [u] = (u,Tu)+f V\u \ 2 + i 

J R d 1 


u(x)\ 2 w(x - y)\u{y)\ 2 dxdy, 


its infimum being denoted en- In the sequel of this appendix we prove the following 
statement, which is a weakened version of Theorem 13.61 


Theorem A.2 (Derivation of Hartree’s theory, alternative statement). 

Let V and w satisfy the assumptions of Section 14.21 in particular ([3.121) . We moreover 
assume that T has a positive kernel in the sense of Definition IA.il or is a limit case of 
this definition, as in (IA.4D . We then have 


lim 

TV—>• oo 


E(N) 

A 


= e H - 


Let Tv > 0 be a ground state of H^, achieving the infimum in (IA.5D and 
Pn\x i,...,x n ):= / |Tat {xi, ... ,x N )\ 2 dx n+ i ... dx N 

J^d(N-n) 

be its n-body reduced density. There exists a probability measure y on M.u, the set of 
minimizers of £h (modulo a phase), such that, along a subsequence 

lim p$ = [ |tt lg>n | 2 dy(u) for all n € N, (A.6) 

N ^-°° Jmh 

strongly in L 1 (M dn ) . In particular, if en has a unique minimizer (modulo a constant 
phase), then for the whole sequence 


lim = |u § n | 2 for all n G N. (A.7) 
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Remark A.3 (Uniqueness for Hartree’s theory). 

In the case of a kinetic energy with positive kernel, uniqueness of the minimizer of £h is 
immediate if w > 0. Indeed, the kinetic energy is then a strictly convex functional 

of p = y/\W> see H22J Chapter 7]. If in is positive the functional <?h is thus itself strictly 
convex as a function of p. □ 

The particular case where T = — A has been dealt with by Kiessling in |100j . and we 
shall follow his method in the general case. It consists in treating the problem as a purely 
classical one, which explains that we only obtain the convergence of reduced densities (1A.6I) 
instead of the convergence of the full reduced density matrices (13.191) . We pursue in the 
direction of (1A.5D by writing 

E(N) = inf [S N [y^] , /jl n € V s (R dN )} (A.8) 

where fi N plays the role of | 'I/at | 2 and we used the fact that we may assume > 0. 
The object we have to study is a symmetric probability measure of N variables, and our 
strategy shall be similar to that used to prove Theorem (12.61) : 

• Since the system is confined, one may easily pass to the limit and obtain a problem 
in terms of a classical state with infinitely many particles /z <E V s (R m ). We 
then use Theorem m to describe the limit /z( n ) of for all n, using a unique 
probability measure € V(V(R d )). 

• The subtle point is to prove that the limit energy is indeed an affine function of /z. 
This uses in an essential way the fact that the kinetic energy has a positive kernel 
(or is a limit case of such energies), as well as the Hewitt-Savage theorem. 

These two steps are contained in the two following sections. We then briefly conclude 
the proof of Theorem IA.2I in a third section. 


A.2. Passing to the limit. 

The limit problem we aim at deriving is described by the functional (compare 
with (12.501) 1 


E\n\ := limsup — T ( J /z( n ) J 
n —kxd Tl \ J 

+ [ V{x)d^ [l \x)+ \ j j w(x - y)dn {2 \x,y), (A.9) 

J R d ^ J JR d xR d 

where /z 6 V s (W m ) and we set 


T (y/K) ■ = 


V k 1 ni 




(A.10) 


for all probability measure /z n £ T 15 (M.^”"). In fact we shall prove in Lemma (1A.5I) below 
that the limsup in (1A.9D is actually a limit. 


Lemma A.4 (Passing to the limit). 

Let /Zjy € V s {R dN ) achieve the infimum in (1A.8D . Along a subsequence we have 

/zj? /*(") € P s (M dri ) 
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for all n € N, in the sense of measures. The sequence (//”)) N defines a probability 
measure t 1 £ P s (M rfN ) and we have 

inf >£[/*]• ( A - 11 ) 

N^rOO 1\ 

Proof. Extracting convergent subsequences is done as in Section 12.31 The existence of the 
measure /i € V s (Mh m ) follows, using Kolmogorov’s theorem. 

Passing to the liminf in the terms 

1 n 

tfE/ V{xj)dn N {x i,...,x N )= V(x)dp$(x) 

A ■, jRdN J R d 

J = 1 

and 

11 ]\r 1 

T7T7—T E / w{x i -x j )dn N {x 1 ,...,x N ) = - j j w(a;-y)dn%\x,y) 

N ^ ~ 1 dR dN IJjRdxRd. 

uses the same ideas as in Chapters [2] and [3j We shall not elaborate on this point. 

The new point is to deal with the kinetic energy in order to obtain 

lim inf -j-T (VaTJv) > lim SU P f \fpW) • (A. 12) 

N —>oo 1\ n—>-oo fl \ / 

To this end we denote Tat = ^/T/v a minimizer in (IA.5I) and we then have 

" N 

_J =1 
n 

Erf 

j=i 

(n) 

where 7 ]^' is the n-th reduced density matrix of Tat. It is a trace-class operator, that we 
decompose under the form 

+OO 

TN ] =J2 X n\ U n)( u n\ 
k =1 

with normalised in and = 1. Inserting this decomposition in the 

previous equation and recalling (1A.10I) we obtain by linearity of the trace 


= 1 tv 


= -Tr 
n 



1 

n 


+00 

X>1 


k=\ 


T 



> 



n 



n 
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where the inequality in the second line uses the convexity of the kinetic energy as a 
functions of the density p, already recalled in Remark I A. 31 cf [ 122L Chapter 7]. In the last 
equality 


CO 


T/v 


/ j ''n I “n I 
k =1 


_ _ \ \k \ k |2 

/h,( n ) / , Ai I A 


is the densit }0 of 7 ^ and it is easy to see that 


_ ( n ) 

P 7 <>) — Pn 


L 


which yields 


Jjd(iV-n) 


±T(^)>±-T 


... ,x N )\ 2 dx n+ i... dx N , 


A 


To obtain (IA.12D . we first pass to the liminf in N (using Fatou’s lemma), then to the 
limsup in n. □ 

The notion of kinetic energy with positive kernel is already crucial at this level. It 
provides the convexity property that we just used. It will play an even greater role in the 
next section. 


A.3. The limit problem. We now have to show that the functional (IA.9I) is affine on 
V s (M. df3 '). The last two terms obviously are, which is not suprising since they are classical 
in nature. It thus suffices to show that the first term, which encodes the quantum nature 
of the problem, is also linear in the density: 

Lemma A. 5 (Linearity of the limit kinetic energy). 

The functional 


is affine on V s 


T {y/JT) := limsup — T 

n—>■ 00 Tl 

pdN\ 


) = lim — T 

n—>-oo n 




(n) 


Kiessling [100] gave a very elegant proof of this lemma in the case of the non-relativistic 
kinetic energy. He notes that in this case 

— 2 


1 . 


n 


/xv 


(n) = 


[ 

n J 


Vj-a/m 


t( n ) 


— V [ 

4n J R‘ 


Vj log n 


(n) 




(n) 


and that the last expression is identical to the Fisher information of the probability mea¬ 
sure /x ^ n \ The quantity we study can thus be interpreted as a “mean Fisher information” 
of the measure /x G V S (W^), in analogy with the mean entropy introduced in (12.501) . 

This quantity has an interesting connection to the classical entropy of a probability 
measure. Letting /x ( n > evolve following the heat flow, one may show that at each time 
along the flow, the Fisher information is the derivative of the entropy. Since the heat 
flow is linear and the mean entropy is affine (cf Lemma 12.71) , Kiessling deduces that the 
mean Fisher information is affine. Another point of view on this question is given in m 
Section 5]. 


27 


Formally, the diagonal part of the kernel. 
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Here we shall follow a more pedestrian approach whose advantage is to adapt to the 
general kinetic energies described in Definition I A. II among other to the pseudo-relativistic 
kinetic energy. 


Proof. Theorem 12. II implies that V s (W m ) is the convex envelop of the symmetric probabil¬ 
ity measures characterized by = p® n , p £ To prove the lemma it thus suffices 

to take 

Hi = P?\ H 2 = pT » H = | Hi + \P2 
with pi,p 2 €E pi 7 ^ p 2 and to prove that 


«<”>) - \T ( JPP ) - It ( JPi 


A n ) 


< o(n ) —>• 0 


(A.13) 


when n —>• oo. For of the above form (i = 1,2) we clearly have 


T 


(v^F)=nT(v^). 


Thus the above estimate implies both that the functional under consideration is affine, 
and that the limsup is a limit, given for a general /i £ V s (W m ) by 


lim -T (y/Jw) = [ T (VP) dP(p) 
n^-oo n \ v / JpeV(R d ) 

where P is the classical de Finetti measure of /i defined by Theorem 12.II 

We now prove our claim (IA.13D . By symmetry of and in view of Definition IA.ll 
we have to compute 

2 


=n f dX ff T(x,y) J(X) - J (y) 

V v / j iH"- 1 ) JJ R d xR d V V 

where we denote 

X = (xi,...,x n ), Y = (yi,x 2 ■ ■ ■ ,x n ), X = (x 2 , . • • ,x n ). 
We start by claiming that for all X, Y 

\I(X,Y)\:= 

^^\X)-^P(Y) 


(A.14) 




^^\X)-^^\Y) 


1/2 


< C [ Y\p 1 (x j )p2{x j ) 
U= 2 


vViOd) - Vpi(yi) + Vp2(xi) - Vp2(vi) • (a.is) 


We will prove (IA.15D in the case 

p 2 (yi) < P 2 ^i) and pi(yi) < pi(xi). (A.16) 

The adaptation to the other cases is left to the reader. We simplify notation by setting 


ui = yfp u Ui= Jpf n . 
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After expanding the squares we obtain 
2 1(X,Y) = Ui(X)Ui(Y) + U 2 (X)U 2 (Y) 

- jU 2 (X)U 2 (Y) + U 2 (X)U 2 (Y) + U?(X)U*(Y) + U*(X)U?(Y). 

Then 

Ul(X)Ul(Y) + U 2 {X)U 2 {Y) + U 2 (X)U 2 (Y) + U 2 {X)U 2 {Y) 

= (C/i(A)?7i(y) + U 2 (X)U 2 {Y)) 2 +Ul{X)U%{Y)+U%{X)Ul(Y)-2U l {X)U 2 {Y)U 1 {Y)U 2 (X) 
and thus, 

2\I(X,Y)\ = I (Ui(X)Ui(Y) + U 2 {X)U 2 (Y)) 


x 1-, 1 + 


U*(X)U*(Y) + U$(X)U?(Y) - 2U 1 (X)U$(Y)U 1 (Y)U 2 {X) 


{U 1 (X)U 1 {Y) + U 2 (X)U 2 (Y)Y 


< 


Uf(X)U 2 (Y) + U*(X)U?(Y) - 2U 1 (X)U*{Y)U 1 (Y)U 2 (X) 
U!(X)Ui(Y) + U 2 (X)U 2 {Y) 

n ! x , s \ |«i(xi)w 2 (yi) - Uiiy^u^xx )] 2 


u= 2 

n 

~[ui(xj)u 2 (a 
4= 2 

( U 

V = 2 

«i(yi) 2 


I^G/i) (ui(xi) - ui(yi)) + ni(^i) {u 2 {yi) - u 2 (zi))|' 
ui(xi)ui(yi) + u 2 (xi)u 2 (yi) 

u 2 (yi ) 2 


+ 


ui{xi)ui(yi) + u 2 (xi)u 2 (yi) 
(u 2 {yi) - u 2 {x 1)) 2 


(ui(xi) - ui(yi)Y 


+ u 2 (xi)u 2 (yi) 

Estimate (IA.15I) immediately follows in the case (1A. 161) and by similar considerations in 
the other cases. 

Inserting ()A.15I) in (1 A. 141) and recalling (1A. 101) we obtain 

* \ n— 1 


/T 


(n) ) _ 


m!” 1 


— -T 
2 ± 


/4 n) 


<Cn( yYply/p 2 (T (VpT)+T (v^pi)) 


(A.17) 

by Fubini. We may assume that T (^/pi) an d ^ (v^) are finite otherwise all the quan¬ 
tities we estimate are equal to +00 and there is nothing to prove. There remains to note 
that 


S:= f yfpis/p2 < r ( [ Pl+ [ P2 
Jm d * \J R d 

since pi 7 ^ p 2 by assumption. We conclude that 

V^) - y (v^F) - Y (y, 


T 


and 5 n 1 —>• 0 when n —» + 00 , as desired. 


4”’ 


< 1 


< Cn5 n ~ l 
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To deal with the case of the usual, non-relativistic, kinetic energy, we first apply the 
preceding arguments to the kinetic energy T s defined by the positive kernel 

T s (x,y) = C d \x - y\~ [d+2s) 

with 0 < s < 1 fixed, to obtain the analogue of (1A.17I) with T = T s . We then multiply 
this inequality by (1 — s) and take the limit s —>• 1. Using (1A.4I) this gives (IA.17I) with T 
the non-relativistic kinetic energy. We may then pass to the limit n —>• oo to obtain the 
desired result. □ 

A.4. Conclusion. 

The upper bound follows by the usual trial state argument. Combining Lemmas IA.4I 
and IA.5I as well as the representation of p, given by Theorem 12.11 we deduce 


lim inf 

iV—>• oo 


E(N) 

N 


= £ 


I 

Jv{ «. d ) 


V(R d ) 

£[p 1 


P®°°dP{p) 


dP{ P ) 


= £h [yfp\ d p(p) > e H dP(p) = e H , 

Jv{ R d ) Mr*) 

where P is the de Finetti measure. This gives the energy convergence. The convergence 
of reduced densities follows as usual by noting that equality must hold in all the previous 
inequalities. 
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Appendix B. Finite-dimensional bosons at large temperature 

Until now we have only considered mean-field quantum systems at zero temperature, 
and obtained in the limit N —>• oo de Finetti measures concentrated on the minimizers of 
the limit energy functional. It is possible, taking a large temperature limit at the same 
time as the mean-field limit, to obtain a Gibbs measure. In this appendix we explain this 
for the case of bosons with a finite dimensional state-space, following [SHI 119] and [1041 
Chapter 3]. 

In infinite dimension, important problems arise, in particular for the definition of the 
limit problem. The non-linear Gibbs measures one obtains play an important role in 
quantum held theory [HU 1181111931 ITT-, 8l] and in the construction of rough solutions to 
the NLS equation, see for example [1101122112311195112911301128111941IH2] . We refer to the 
paper [Il4| for results on the “mean-held/large temperature” limit in infinite dimension 
and a more thorough discussion of these subjects. 


B.l. Setting and results. 

In this appendix, the one-body state space will be a complex Hilbert space Sj with finite 
dimension 

dimf) = d. 

We consider the mean-held type Hamiltonian 

N 1 

h n = ^2 h j + —— y wij (B.i) 

j =1 l<i<J<AT 

where h is a self-adjoint operator on S) and w a self-adjoint operator on fj (g>f), symmetric 
in the sense that 

w(u (8) v) = w{y G> u), Vu, v € Sj. 

The energy functional is as usual dehned by 

= i'&N, Hn^n) 

for Ttv € and extended to mixed states Ttv over fj A = 0^15 by the formula 

£7v[riv] = Tr^jv [H N F N ] . 


The equilibrium state of the system at temperature T is obtained by minimizing the 
free-energy functional 

•FjVpMv] := ^Ar[Tjv] + TTV [Ttv logTiv] (B.2) 


amongst mixed states. The minimizer is the Gibbs state 

= exp (-T-'Hn) 

N Tr [exp (-T-'H/v)] ‘ 

The associated minimum free-energy is obtained from the partition function (normalisation 
factor in (IB.31) 1 as follows: 


Ttv = inf {Fn[Tn], Ftv € 5^)} 


—T log Tr 



(B.4) 


We shall be interested in the limit of these objects in the limit 

N oo, T = tN, t hxed 


(B.5) 
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which happens to be the good regime to obtain an interesting limit problem. We will in 
fact obtain a classical free-energy functional that we now define. 


Since Sj is finite-dimensional, one may define du, the normalized Lebesgue measure on 
its unit sphere SSj. The limiting objects will be de Finetti measures, hence probability 
measures n on SSj, and more precisely functions of L^STj, du). We introduce for these 
objects a classical free-energy functional 


Fal/i] = / £h [u]fi(u)du + t / n(u) log (n(u)) du 

JSSj J SSj 


(B. 6 ) 


and denote F c \ its infimum amongst positive normalized L 1 functions, 
the Gibbs measure 

_ exp (-G^hM) 
fssj exp (—du 


It is attained by 


(B.7) 


and we have 


^1 




j£h[«] 



Here £n[ u ] is the Hartree functional 

£n[n\ = (u® N , H n u® n )^ n = (u , hu)^ -\-\{u®u,wu® u) f)2 ■ (B. 8 ) 

The theorem we shall prove, due to Gottlieb [85j (see also [8T1I95] and, independently j!04 . 
Chapter 3]) is of a semi-classical nature since it provides a link between the quantum and 
classical Gibbs states, (IB.311 and (IB.71) : 


Theorem B.l (Mean-field/large temperature limit in finite dimension). 

In the limit (USD. we have 

F n = -T log dim (ft?) + NF C i + O(d). (B.9) 

Moreover, denoting 7 ^ the n-th reduced density matrix of the Gibbs state (IB.31) . 

-> [ | u® n )(u® n \ii cl (u)du (B.10) 

J SSj 

strongly in the trace-class norm of6 1 (Ij n ). 


Remark B.2 (Mean-field/large temperature limit). 

A few comments: 

(1) One should understand this theorem as saying that essentially, in the limit under 
consideration, 

£n ~ [ \u 0N )(u (Z ' N \n c i(u)du. 

JSS) 

The Gibbs state is thus close to a superposition of Hartree states. The notions of 
reduced density matrices and de Finetti measures provide the appropriate manner 
to make this rigorous. We will see that the de Finetti measure (lower symbol) 
associated to Tat by the methods of Chapter [4] converges to fj, c \ (u)du. 
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(2) Note that the first term in the energy expansion (IB.91) diverges very rapidly, 
see (14.171) . The classical free-energy only appears as a correction. In view of the 
dependence on d of this first term, it is clear that the approach in this appendix 
cannot be adapted easily to an infinite dimensional setting. 

(3) Our method of proof differs from that of [83] • We shall exploit more fully the semi- 
classical nature of the problem by using the Berezin-Lieb inequalities introduced 
in [161 12(')[ lt'8'4] . The method we present [I IT) owes a lot to the seminal paper [120] 
and is reminiscent of some aspects of m- 

(4) It will be crucial for the proof that the lower symbol of Fn is an approximate de 
Finetti measure for Tjy. This will allow us to apply the first Berezin-Lieb inequality 
to obtain a lower bound to the entropy term. A new interest of the constructions 
of Chapter |4] is thus apparent here. We use not only the estimate of Theorem 14.11 
but also the particular form of the constructed measure. 

□ 


B.2. Berezin-Lieb inequalities. 

We recall the resolution of the identity (14.51) over S) N given by Schur’s lemma. We thus 
have for each state Tat € S(Sj^) a lower symbol defined as 

fi N {u) = dim (fjf) Tr [r A r|u 0jV )(u® Ar |] . 

The first Berezin-Lieb inequality is the following statement: 

Lemma B.3 (First Berezin-Lieb inequality). 

Let Tat € S(Sj^) have lower symbol and f : M + —> M be a convex function. We have 

TV [/(I»] > dim (fi?) (B.ll) 

The second Berezin-Lieb inequality applies to states having a positive upper symbol 
(see Section fOli . One may in fact show that every state has an upper symbol, but it is 
in general not a positive measure. 


Lemma B.4 (Second Berezin-Lieb inequality). 

Let Tat € S(S)N) have upper symbol /r/v > 0, 

Tat= f \u® N )(u® N \p.]\r(u)du 
Ju&SS) 

and f : R + —> R be a convex function. We have 

Tr[f(T N )]<6mi{f>T) J f (- 


Vn 


dim (fj 


ad 


du. 


(B.12) 


(B.13) 


Proof of Lemmas \B.<A and B.4 We follow [ 184] , Since Tat is a state we decompose it under 
the form 




k=1 


with normalized and = 1. We denote 
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and by (14.51) we have 


dim (ijf) I n k N (u)du = (v$, V#\ = 1. 
JSS) N ' 


On the other hand, since {V^)k is a basis of fj;/, for all u € SSj 


£*&(«) = £ \(vl^ N ) 


= 1. 


(B.14) 


(B.15) 


First inequality. Here hn is the lower symbol of Tjv and we have 

Hn(u) = dim (ijf) X nFn(u) 


and thus 


dim (ijf) f , 
Jss .j 


FN 


du < dim (fjf) / ^ / (\ k N ) fi k N (u)du 

k 


dim (SjN 

by Jensen’s inequality and (IB.15I) . Next 

dim (ijf) f ^Tf (a&) fJ, k N (u)du = E / (a = Tr [f(T N )} 
Jss i k k 

using (IB. 141) . 

Second inequality. Here Tjv and are related via (IB.121) . We write 
Tr [/(IV)] = E / (A k N ) = E / 


Fn{u) 


= E/f/ HN{v)n k N (u)d\ 
<Sdim W)/ M / VdimW) 

=dk(,f) l'(di 


ia k N (u)du 


Vn(u) 
dim (Sj?) 


using Jensen’s inequality and (IB.14D to prove the inequality in the third line and 
then (IB. 151) to conclude. □ 


We have here presented a specific version of these general inequalities. It is clear that 
the proof applies more generally to any self-adjoint operator on a separable Hilbert space 
having a coherent state decomposition of the form (14.51) . In the next section these inequali¬ 
ties will be used to deal with the entropy term by taking f(x) = x log x. This will complete 
the treatment of the energy using Theorem 14.11 and make the link with the discussion of 
Chapter 0] 
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B.3. Proof of Theorem IB.lL 
Upper bound. We take as a trial state 

rt st := [ \u® N )(u® N \ii c i(u)du. 

JSS 3 

The energy being linear in the density matrix we have 

Sn[tT}= f £ N [\n® N )(u® N \] f i cl (u)du = N [ £ H [nMu)d 
JSSj JSS) 

For the entropy term we use the second Berezin-Lieb inequality, Lemma IB. 41 with /(x) 
xlogx. This gives 


u. 


Tr [r^ st log r^ st ] < dim (tff) 


<SS) 


dd (u) 
dim (fj ^' 


log 


dd (u) 

dim 


du 


log dim (fjf) + / /i c iO) log (Mci(u)) du. 
Jss i 


ISS) 

Summing these estimates we obtain 

F n < Fn [r^ st ] < -T log dim (iof) + NF C \ 
since fj, c \ minimizes F c \. 

Lower bound. For the energy we use the reduced density matrices as usual to write 


e N \p N \=NTig> h^ ] 

Denoting 

/LvO) = dim (f)f) (u® N ,T n u® n ) 

the lower symbol of Tjv, we recall that it has been proved in Chapter [I] that 

d 

\u){u\HN{u)du < C'i 
'SSj 

\u® 2 )(u® 2 \nN(u)du 


N m 

+ T7 Tr fi 2 


( 2 ) 

w 1n 


Tit 


Tic 


( 2 ) 

IN 


<SS) 


< C\ 

< c 2 


N 

d_ 

N' 


Since we work in finite dimension, h and w are bounded operators and it follows that 

£jv[riv] > N f Ti% [h\u){u\]/j,N(u)du + — f Tt^ [w\u® 2 )(u® 2 \] [ii\r(u)du — Cd 
Jss) 2 Jss) 

= N S}i[u\nN(u)du — Cd. 

Jss) 

To estimate the entropy we use the first Berezin-Lieb inequality, Lemma IB. 31 with /(x) 
xlogx. This yields 


Tr [r N log T n ] > dim ) 


>SS) 


dim ' 


log 


m(u) 

dim 


du 


= -logdim (fjf) + / /j, N (u)log(fj, N (u))du. 

JSS) 
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There only remains to sum these estimates to deduce 

Fn = Fn\T n ] > -T log dim (iof ) + NT cl [p N ] - Cd 
> -T log dim (ijf) + NF C \ - Cd 


since 


by definition. 



Convergence of reduced density matrices. The lower symbol fj,j\r(u)du is a probability 
measure on the compact space SSj. We extract a converging subsequence 


HN{u)du —> n{du) € V(SSj) 

and it follows from the results of Chapter 0] that, for all n > 0, along a subsequence, 

—»■ [ \u® n ){u® n \dff{u). (B.16) 

JSSj 

Combining the previous estimates gives 

F c i>F cl [fi N ]-C^. (B.17) 

To pass to the liminf N —> oo in this inequality, the energy term is dealt with as in the 
preceding chapters. For the entropy term we note that since du is normalized 



Hn log fj,Ndu > 0 


since this quantity may be interpreted as the relative entropy of rn with respect to the 
constant function 1. Using Fatou’s lemma we thus deduce from (IB.17I) that 


F c \ > 


and thus dfj,(u) = fi c \(u)du by uniqueness of the minimizer of J- c \. Uniqueness of the 
limit also guarantees that the whole sequence converges and there only remains to go back 
to (IB. 161) to conclude the proof. □ 
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